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the  nonunifonn  flow  field,  the  store  forces  and  moments,  and  the  store 
trajectory  are  presented.  Comparisons  between  theory  and  experiment  for  flow 
fields,  store  loading  distributions,  store  forces  and  moments,  and  store 
trajectories  are  shown  and  discussed. 
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This  report  "Prediction  of  Supersonic  Store  Separation 
Characteristics  Including  Fuselage  and  Stores  of  Noncircular 
Cross  Section,"  describes  a  combined  theoretical-experimental 
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data.  Volume  II,  "Users  Manual  for  the  Computer  Program," 
presents  detailed  instructions  on  the  use  of  the  computer 
program  with  emphasis  on  preparation  of  input  data  and  inter¬ 
pretation  of  output.  Volume  III,  "Appendices  A  and  B, 

Details  of  Program  I,"  provides  additional  descriptions  of 
the  individual  subroutines  and  program  variables  passed  between 
modules  in  the  first  of  two  programs.  Volume  IV,  "Appendices 
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of  the  individual  subroutines  and  program  variables  passed 
between  modules  in  the  second  program. 
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PREDICTION  OF  SUPERSONIC  STORE  SEPARATION  CHARACTERISTICS 
INCLUDING  FUSELAGE  AND  STORES  OF  NONCIRCULAR 

CROSS  SECTION 


Volume  I.-  Theoretical  Methods  and  Comparisons 

with  Experiment 


1 .  INTRODUCTION 


This  report  presents  the  results  of  a  combined  theoretical/ 
experimental  research  program,  the  purpose  of  which  is  to 
extend  and  improve  the  supersonic  store  separation  prediction 
method  described  in  References  1  and  2.  This  volume  presents 
the  theoretical  methods  used  and  compares  predictions  with 
experiment.  The  computer  programs  developed  during  the  present 
investigation  are  described  in  Volumes  II  through  IV  of  this 
report.  References  3  through  5. 

The  aircraft  configuration  treated  in  References  1  and 
2  was  limited  to  an  axisymmetric  fuselage,  a  wing  without 
dihedral  or  incidence  relative  to  the  fuselage,  a  single  pylon 
under  the  wing  or  fuselage,  and  single  axisymmetric  store 
which  could  have  planar  or  cruciform  fins.  The  forces  and 
moments  acting  on  the  store  were  calculated  using  slender-body 
theory . 

In  the  present  work  the  axisymmetric,  or  circular,  fuselage 
has  been  retained  as  an  option.  A  second  option  has  been 
added  in  order  that  a  fuselage  which  is  noncircular  in  cross 
section  can  be  modeled  with  or  without  engine  inlets.  The 
wing  model  has  been  generalized  to  include  dihedral,  which  may 
vary  across  the  span,  and  wing  incidence  relative  to  the 
fuselage.  Incidence  is  only  included  in  the  wing  boundary 
condition  for  reasons  to  be  discussed.  The  pylon  model  has 
been  extended  to  allow  for  breaks  in  leading-edge  and  trailing- 
edge  sweep.  The  option  of  adding  an  ejector  rack  to  the  pylon 


has  been  added.  Up  to  seven  stores  may  now  be  used  on  the 
parent  aircraft.  They  may  be  circular  or  elliptic  in  cross 
section.  In  the  trajectory  calculation,  a  store  modeled  using 
the  circular  option  may  have  one  set  of  planar  or  cruciform 
fins.  If  the  elliptic  option  is  used,  two  sets  of  fins  may 
be  included  with  one,  two,  three,  or  four  fins  ,in  each  set.  For 
certain  fin  arrangements  tracking  of  the  trailing-edge  vortices 
from  the  forward  set  of  fins  to  the  aft  set  may  be  performed 
so  their  influence  can  be  included  in  the  loads  on  the  aft 
set  of  fins. 

The  experimental  program  carried  out  in  conjunction  with 
the  present  work  will  not  be  described  in  this  report.  It  is 
fully  documented  in  Reference  6.  The  purpose  of  the  program 
was  to  provide  systematic  data  at  supersonic  speeds  for  use 
in  evaluating  and  improving  analytical  techniques  for  predict¬ 
ing  store  separation  trajectories.  The  program  was  conducted 
under  Air  Force  Contract  F33615-75-C-3053  and  the  present  con¬ 
tract.  Five  wind-tunnel  entries  were  made  during  the  period 
from  June  1975  through  February  1978.  The  test  program  was 
conducted  in  Tunnel  A  of  the  von  K£rm£n  Gas  Dynamics  Facility 
(VKF)  at  Arnold  Engineering  Development  Center  (AEDC)  at 
nominal  Mach  numbers  of  1.5,  1.6,  2.0,  and  2.5  with  most  testing 
done  at  1.5  and  2.0.  Four  types  of  data  were  obtained;  flow- 
field  survey  data,  store  pressure-distribution  data,  store 
force-and-moment  data,  and  store  separation  trajectories. 
Numerous  parent  aircraft  configurations  were  used.  They  varied 
from  a  simple  circular  fuselage  with  or  without  a  wing  to  a 
complex  noncircular  fuselage  with  engine  inlets,  pylon,  rack, 
and  stores.  The  parent  aircraft  configurations  were  built  up 
component  by  component  in  order  to  isolate  interference  effects. 

The  next  two  sections  of  this  report  describe  the  general 
approach  to  the  problem  and  the  flow  models  based  on  linear, 
potential  flow  theory  used  for  the  aircraft  components  and  the 
stores.  The  work  of  Reference  1  will  be  heavily  referenced  with 
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only  unifications  and  additions  described  in  detail.  This 
will  be  followed  by  a  section  which  describes  nonlinear 
corrections  made  to  the  flow  models  in  order  to  place  shock 
waves  more  nearly  in  their  correct  positions. 

The  last  sections  of  the  report  describe  the  parent  air¬ 
craft  flow  field  calculation,  the  store  force  and.  moment 
calculation,  including  store-aircraft  interference  effects, 
and  the  store  trajectory  calculation.  Many  comparisons  between 
experiment  and  theory  are  also  presented. 
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2 .  GENERAL  APPROACH 


It  is  possible  to  determine  the  forces  and  moments  acting 
on  an  external  store  in  the  presence  of  a  parent  aircraft  in 
the  following  manner.  Consider  first  the  entire  airframe 
excluding  the  store  in  question.  The  nonuniform  flow  field, 
induced  by  the  flow  model  representing  the  airframe,  is  com¬ 
puted  in  the  region  occupied  by  the  store.  If  the  store  is  now 
placed  in  thxs  nonuniform  flow  field,  then  the  calculated 
loading  acting  on  it  includes  primary  interference.  The 
primary  interference  can  be  viewed  as  the  first  term  in  an 
iterative  procedure.  Higher-order  solutions  would  be  generated 
as  follows.  After  modeling  the  store  subjected  to  the  primary 
interference  flow  field,  its  effects  on  the  airframe  can  be 
determined.  These  effects,  in  terms  of  velocities  normal  to  the 
wing,  pylon,  and  fuselage,  must  then  be  cancelled  by  an  additional 
set  of  singularities  distributed  over  these  components.  A  new 
nonuniform  flow  field  can  now  be  computed  in  the  region  occupied 
by  the  store,  and  the  loading  acting  on  it  can  be  recalculated. 

At  this  stage,  the  result  is  associated  with  the  second  term  in 
the  iteration  procedure.  The  difference  between  the  results  of 
the  second  and  the  first  loading  calculation  is  due  to  additional 
interference.  This  process  can  be  continued  until  the  effects  of 
additional  interference  becomes  negligible. 

In  the  work  of  Reference  1,  a  modification  to  the  first 
iteration  was  used  to  account  for  aircraft-store  mutual  inter¬ 
ference.  This  consisted  of  including  in  the  first  iteration 
the  effects  of  the  store  in  the  boundary  condition  applied 
to  the  aircraft  flow  model.  This  approach  has  two  difficulties. 
First,  it  is  impractical  to  include  enough  panels  in  modeling  the 
aircraft  to  describe  accurately  the  multiple  reflections  of  the 
store  shock  wave  between  the  store  and  the  airplane.  Second, 
the  separating  store  moves  relative  to  the  aircraft  and  the 
point  at  which  its  shock  wave  strikes  the  aircraft  moves. 
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changing  the  boundary  conditions.  This  requires  that  the 
strengths  of  the  singularities  modeling  the  wing-pylon- 
fuselage  combination  be  recalculated  at  each  point  in  the 
trajectory,  a  lengthy  process.  It  was  recommended  that 
store  effects  be  removed  from  the  aircraft  boundary  condition 
and  store-aircraft  interference  be  accounted  for  by  imaging. 

The  wing  and/or  fuselage  are  treated  as  reflection  planes; 
that  is,  an  image  store  is  placed  on  the  other  side  of  the 
wing  and/or  inside  the  fuselage.  This  approach  has  been 
adopted  in  the  present  work. 

The  approach  taken  in  the  present  work  is  to  model 
all  bodies  (fuselage,  ejector  rack,  and  stores)  as  isolated 
bodies.  Singularities  are  used  to  model  both  volume  and 
angle-of-attack  effects.  Wing  and  pylon  thickness  are  modeled 
by  source  panels  and  the  lifting  surfaces  are  modeled  by 
u-velocity  panels.  Mutual  interference  between  wing  and  pylon 
is  accounted  for  and  their  mutual  interference  with  the 
fuselage  is  included  by  placing  an  interference  shell  of 
u-velocity  panels  on  the  fuselage. 

With  the  singularity  distributions  which  represent  the 
various  aircraft  components  determined,  the  trajectory  of  the 
separated  store  can  be  calculated.  This  requires  that  the 
forces  and  moments  acting  on  the  store,  including  damping, 
be  calculated  at  each  point  in  the  trajectory.  This  is  done 
by  calculating  the  three-dimensional  flow  field  in  which  the 
store  is  located.  To  do  this  the  separated  store  is  removed 
from  the  flow  field  and  the  perturbation  velocity  field  due 
to  the  aircraft  components  is  calculated  at  specified  points 
in  the  region  the  store  had  occupied.  These  points  are  the 
control  points  to  be  used  in  modeling  the  store  in  the  non- 
uniform  flow  field.  The  store  is  then  put  back  in  the  flow 
field  and  the  image  store (s)  is  located  relative  to  the  store. 
The  free-stream  velocity  components  that  the  store  sees  and  the 
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velocities  due  to  store  angular  motion  are  added  to  the 
perturbation  velocities. 

The  solution  for  the  store  singularity  strengths 
begxns  at  the  store  nose  and  works  its  way  aft  one  axial 
station  at  a  time.  As  soon  as  an  image  store  influence  is 
felt  at  a  control  point,  the  image  store  induced  velocities 
are  added  to  the  previously  calculated  velocities.  After 
all  of  the  singularity  strengths  have  been  calculated,  the 
store  surface  pressure  distribution  is  calculated  using  the 
full  Bernoulli  equation.  These  pressures  are  integrated 
radially  and  axially  to  determine  the  forces  and  moments 
acting  on  the  store.  A  similar  procedure  is  used  to  cal¬ 
culate  the  empennage  forces  and  moments. 
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3.  FLOW  MODELS 


In  general,  the  flow  models  for  the  components  to  be 
described  are  based  on  linear  theory  for  steady,  inviscid 
flow.  Because  the  flow  regime  considered  here  is  supersonic, 
the  associated  potential  flow  is  governed  by  the  wave 
equation.  In  a  rectangular  coordinate  system  (x,y,z)  with 
the  x-axis  in  the  direction  of  the  free-stream  vector,  the 
potential  <p  is  given  by 


n  ^4  -  2%  -  *4 

3xz  3y  3z 
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(1) 


This  equation  is  valid  if  the  perturbation  velocities  u,v,w 
given  below  in  terms  of  the  potential  are  small. 
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Because  of  the  linearity  of  the  potential  $  in  Equation  (1) , 
it  is  possible  to  make  use  of  superposition.  This  principle 
allows  for  the  representation  of  the  aircraft  components 
by  distributions  of  potentials  or  singularities.  Flow  con¬ 
ditions  at  a  given  field  point  are  given  by  the  summed 
influences  of  all  the  singularities.  Before  describing  the 
flow  models,  the  aircraft  components,  the  store,  and  flight 
conditions  will  be  summarized. 


3.1  Description  of  Aircraft  Components  and 
Flight  Conditions 

The  configuration  of  interest  consists  of  a  wing  attached 
to  a  fuselage.  A  pylon  can  be  positioned  under  the  wing  or 


under  the  fuselage  centerline.  An  ejector  rack  can  be 
attached  to  the  pylon.  Up  to  seven  stores  can  be  positioned 
relative  to  the  aircraft. 

The  geometric  parameters  of  the  aircraft  components 
described  below  are  accounted  for  in  the  computer  program. 
Flight  conditions  cover  the  supersonic  speed  regime  and  are 
specified  below  together  with  the  expected  range  of  validity 
of  the  present  method. 

3.1.1  Fuselage 

Two  options  exist  for  modeling  the  fuselage.  These 
options  will  be  referred  to  throughout  this  report  as  the 
"circular"  fuselage  option  and  the  "noncircular"  fuselage 
option. 

In  the  circular  fuselage  option,  the  body  is  axisymmetric 
and  its  axial  radius  distribution  is  specified  by  a  set 
of  polynomials. 

The  noncircular  fuselage  may  have  an  arbitrary  cross 
section  which  may  vary  over  the  length  of  the  fuselage.  A 
number  of  options  exist  in  the  computer  program  for  specifying 
the  cross-sectional  shape  as  a  function  of  axial  location. 

All  of  the  options  provide  information  such  that  quadrilateral 
panels  can  be  laid  out  on  the  fuselage  surface.  This  option 
also  allows  for  the  inclusion  of  external-compression-ramp 
inlets  on  the  sides  of  the  fuselage.  The  mass-flow  ratio 
through  the  inlets  is  varied  by  specifying  panels  on  the  face 
of  the  inlet  as  "open"  or  "closed". 

3.1.2  Wing 

The  wing  may  have  twist  and  camber  and  the  leading  and 
trailing  edges  can  have  breaks  in  sweep.  Dihedral  angle  is 
also  variable  across  the  wing  span  and  incidence  relative  to 
tne  fuselage  centerline  may  be  specified  and  is  included  in  the 
wing  boundary  condition.  Vertical  location  of  the  wing 
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attachment  on  the  fuselage  is  not  limited  to  the  midplane. 
Effects  of  wing  thickness  are  included  in  the  program.  The 
wing  leading  edge  can  be  blunt  or  sharp. 

3.1.3  Pylon 

Leading  and  trailing  edges  may  be  swept  with  breaks  in 
sweep  across  the  span.  The  length  of  *the  pylon  root  chord 
need  not  equal  the  local  wing  chord.  Thickness  is  accounted 
for  and  the  leading  edge  can  be  blunt  or  sharp.  The  pylon 
can  be  located  under  the  wing  or  fuselage  but  is  aligned  with 
the  free  stream. 

3.1.4  Ejector  Rack 

The  ejector  rack  is  assumed  to  be  attached  to  the  pylon 
and  can  only  be  used  with  a  pylon.  Only  the  rack  body  is 
modeled  and  this  is  modeled  as  an  axisymmetric  body.  Its 
axial  radius  distribution  is  specified  by  a  set  of  poly¬ 
nomials. 

3.1.5  Stores 

Two  options  exist  for  modeling  stores.  They  will  be 
referred  to  throughout  this  report  as  the  "circular"  store 
option  and  the  "elliptical"  store  option.  Up  to  seven  stores 
can  be  included  in  the  parent  aircraft  model ,  including  the 
store  to  be  separated. 

The  circular  store  body  shape  is  described  by  a  set  of 
polynomials  which  gives  its  radius  as  a  function  of  axial 
location.  Up  to  seven  different  shapes  may  be  used.  A 
particular  snape  may  be  used  for  more  than  one  store. 

The  elliptical  store  option  is  restricted  to  two  dif¬ 
ferent  shapes;  again,  more  than  one  store  can  be  modeled  with 
a  given  shape.  Several  options  exist  for  describing  the 
shape  of  the  elliptic  body.  They  all  provide  information  which 
allow  computation  of  the  body  cross-sectional  shape  as  a 
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function  of  axial  location.  This  information  allows  quadri¬ 
lateral  panels  to  be  laid  out  on  the  body  surface. 

The  store  to  be  separated  may  have  fins.  If  the 
circular  store  option  is  used,  it  can  have  one  set  of  fins. 

They  may  be  planar  or  cruciform  and  may  be  located  axially 
anywhere  on  the  store  body.  Use  of  the  elliptic  store  option 
allows  up  to  two  sets  of  fins  to  be  placed  on  the  store  at 
arbitrary  locations.  Each  set  can  consist  of  one,  two,  three, 
or  four  fins.  The  fins  can  be  located  anywhere  on  the  cir¬ 
cumference  of  the  body. 

All  of  the  stores  may  be  at  arbitrary  pitch  angle  rela¬ 
tive  to  the  aircraft  but  must  be  at  zero  yaw  angle.  Stores 
using  the  circular  body  option  must  be  at  zero  roll  angle 
while  those  using  the  elliptic  option  may  be  at  arbitrary 
roll  angle. 

At  the  beginning  of  the  trajectory  the  store  may  have 
arbitrary  translational  and  rotational  velocities  relative  to 
the  aircraft.  The  store  may  be  ejected  with  an  ejector  force 
versus  time  or  displacement  history.  There  can  be  one  or  two 
ejector  feet.  The  store  can  be  powered  by  specifying  a  thrust 
force  versus  time  history.  Store  mass  and  moments  and  products 
of  inertia  are  constant  and  the  center  of  gravity  may  be 
located  off  the  store  longitudinal  axis. 

3.1.6  Flight  Conditions 

The  flight  path  is  straight  but  not  necessarily  horizontal 
angle  of  attack  is  constant,  and  angle  of  sideslip  is  zero. 

The  Mach  number  should  lie  between  1.2  and  3.0.  The  angle  of 
attack  should  be  in  the  0°-10°  range. 

3.2  Flow  Models  for  the  Fuselage,  Ejector 
Rack,  and  Store  Bodies 

This  section  will  describe  the  flow  models  of  the  fuselage 
rack,  and  stores  as  they  comprise  part  of  the  parent  aircraft. 
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In  the  present  work  they  are  modeled  as  isolated  bodies. 

The  modeling  of  the  separated  store  during  its  trajectory 
when  forces  and  moments  must  be  calculated  will  be  presented 
in  Section  6. 

3.2.1  Circular  Fuselage 

The  potential  flow  method  used  to  model  the  circular 
fuselage  is  the  one  which  was  used  in  Reference  1  with  one 
modification  which  will  be  discussed  later.  The  method  is 
described  in  Section  3.2  and  Appendix  I  of  that  reference. 
Fundamentally,  an  axisymmetric  body  in  supersonic  flow  can  be 
represented  by  a  distribution  of  line  sources  or  sinks  and 
doublets  on  the  body  centerline  to  account  for  volume  and 
angle  of  attack  effects,  respectively.  The  strengths  of  these 
singularities  are  determined  from  the  flow  tangency  conditions 
applied  at  points  on  the  body  surface.  The  source  strength 
constants  are  determined  using  Equations  (1-14)  and  (1-17)  and 
the  doublet  strength  constants  using  Equations  (1-27)  and 
(1-28)  of  Appendix  I  of  Reference  1.  The  strengths  of  the 
singularities  vary  linearly  with  length.  The  perturbation 
velocities  produced  by  the  line  singularities  at  a  field  point 
are  given  by  Equation  (1-30)  in  the  appendix. 

During  the  course  of  the  present  investigation  it  was 
found  that  this  flow  model  worked  well  as  long  as  the  field 
point  at  which  the  velocities  were  to  be  calculated  was  not 
downstream  of  the  point  at  which  the  last  singularities  model¬ 
ing  the  body  began  inf luencing  the  flow  field.  As  the  field 
point  moved  downstream  of  the  point  at  which  the  last  singu¬ 
larities  were  "felt",  the  induced  velocities  increased  without 
bound  due  to  the  linearly  varying  strengths  of  the  singularities. 
To  eliminate  this  problem  additional  singularities  are  intro¬ 
duced  at  the  body  base  to  cancel  the  effects  caused  by  the 
singularities  modeling  the  body.  The  potential  and  velocity 
equations  will  now  be  derived,  first  for  the  sources  and  sinks 
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and  then  for  the  doublets.  The  derivations  parallel  those  in 
Appendix  I  of  Reference  1. 

The  following  sketch  shows  a  body  of  revolution  of 
length  xL. 


s<r° 


*L  V  ? 


Line  sources  and  sinks  of  linearly  increasing  strength  modeling 

the  body  volume  start  at  CQ  =  0,  C-^,  ^2'***^n  known 

source  constants  of  Kq,  K^,  At  the  base  of  the  body, 

x  =  xL,  the  source  strengths  are  KQxL,  (xL  -  » 

K„(xt  -  C_)  ,  ...  K  (xT  -  £  )  .  Therefore,  the  function 
z  z  n  Li  n 

representing  the  singularities  to  be  introduced  at  the  body 
base  to  cancel  the  singularities  behind  the  base  is 


fLU)  -  -koC  -  kx(c  -  q)  -  k2(C  -  Z2) 


■-  Kn<t  -  V 


The  potential  for  these  sources  and  sinks  starting  at 
x_  =  C  =  xT  is 


*a(VrB> 


--f 


XB-6rB 


fL(C)dC 


- t 


Since  the  integrand  is  singular  at  the  upper  limit  of  inte¬ 
gration,  the  procedures  described  on  page  108  of  Reference  1 
can  be  followed  and  the  axial  potential  rewritten  as 


cosh 
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XB  '  XL 
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f L (XB  ~  BrBcosh  a)da  (5) 


For  the  ntla  term  in  Equation  (3)  where  f^  n  =  -Kn(£  "  5n)  » 
Equation  (5)  becomes 


cosh 
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and  after  integrating 


(xB'rB)  =  Kn(xB  ~  Vcosh 


-1 


XB  -  XL 
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The  total  potential  for  all  the  source-sink  singularities 
originating  at  the  body  base  is  obtained  by  summing  the  right 
side  of  Equation  (7)  over  the  N  singularities. 

The  axial  and  radial  perturbation  velocities  can  be 
determined  by  differentiation  of  the  potential 
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Similar  expressions  can  be  derived  for  the  line  doublets 

introduced  at  the  body  base.  Referring  to  the  preceeding 

sketch,  line  doublets  modeling  the  body  angle  of  attack  effects 

start  at  =  0,  ^2'*”^n  ^nown  doublet  constants 

K,  K,  .  At  the  base  of  the  body,  x  *  xT , 

c,o  a,i  a,/  a,n  I* 

the  doublet  strengths  are  Qx L,  ^(xL  -  2  (xL  -  , 

. . . K.  (xT  -  £  ).  The  function  representing  the  doublet 
a ,  n  l  n 

singularities  to  be  introduced  at  the  body  base  is 


dL<«>  "  -Kd,o«  '  Kd,l<5  -  <1»  -  Kd,2<5  -  «2>-  •••  -  Kd,n<5  -  V 
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The  potential  for  these  doublets  starting  at  x0  =  5  =  is 
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This  potential  is  integrated  following  the  same  procedure  used 
for  the  sources  to  give 


<j>c  (xB,rB, 0)  =  -cos  6 
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The  total  potential  due  to  all  the  doublets  introduced  at  the 
body  base  is  obtained  by  summing  the  right  side  of  Equation 
(12)  over  the  N  singularities. 

The  axial,  radial,  and  tangential  perturbation  velocities 
are  obtained  by  differentiating  the  potential, 
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From  Equation  (12) 
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The  velocities  induced  by  all  of  the  line  sources  and 
doublets  starting  at  the  body  base  are  obtained  by  summing 
Equations  (9)  and  (14)  over  the  N  singularities. 

An  examination  of  Equations  (9)  and  (14)  shows  that  the 
last  terms  in  the  equations  for  \  B  and  vfi  are  singular  on  the 
Mach  cone  originating  at  the  body  base,  (x0  -  xL)/$rfi  =  1. 

The  singular  terms  in  Equation  (9)  for  the  source  induced 
velocities  are  multiplied  by  the  source  strength  Kn(xL  -  ?n) • 
For  a  closed  body  the  sum  of  the  source  strengths  is  zero, 
that  is 


N-l 

I  Kn(ltL  -  V  -  ° 


Making  use  of  this  fact,  the  singular  terms  in  Equation  (9) 
have  been  set  to  zero. 

For  the  doublets  the  requirement  that  the  sum  of  their 
strengths  be  zero  does  not  apply.  Thus,  the  singular  terms  in 
Equation  (14)  cannot  be  eliminated  for  this  reason.  A  numerical 
study  of  the  relative  magnitudes  of  the  singular  and  nonsingular 
terms  was  performed.  It  was  found  that  in  a  short  distance 
downstream  of  the  Mach  cone  originating  at  the  body  base  the 
singular  terms  were  small  compared  to  the  nonsingular  terms. 
Thus,  to  avoid  numerical  difficulties  on  or  very  near  the 
base  Mach  cone,  the  singular  terms  in  the  uB  d  and  vfi  d 
expressions  in  Equation  (14)  have  been  set  to  zero.  The  last 
term  in  the  wn  .  equation,  which  is  proportional  to  the  sum  of 

D  f  U 

the  doublet  strengths  £  n(*L  ~  £n)  has  been  set  to  zero 
since  it  was  also  found  to  be  small. 

The  total  axial,  radial,  and  tangential  body  induced 
perturbation  velocities  at  field  point  (xfi,rB,0)  are  obtained 
by  adding  the  velocities  given  by  Equation  (1-30)  of  Appendix  I 
of  Reference  1  to  those  given  by  Equations  (9)  and  (14)  after 
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summing  these  latter  two  equations  over  the  number,  NXBODY, 
line  sources  and  doublets.  The  resulting  equations  are 


nwidf 


and 


xBd  'rBd  “  coordinates  of  the  body  definition  point 
n  •  n 

uB  =  axial  perturbation  velocity  component 
vB  =  radial  perturbation  velocity  component 
Wg  =  tangential  perturbation  velocity  component 

The  positive  directions  are  shown  again  for  convenience  in  the 
following  sketch. 


To  show  how  the  introduction  of  the  source  and  doublet 
singularities  at  the  body  base  improves  the  flow-field  predic¬ 
tions  behind  the  body  base,  calculations  were  performed  for  a 
circular  arc  body  of  length  to  diameter  ratio  (L/D)  of  5.  For 
zero  degrees  angle  of  attack,  sources  only,  the  results  of 
including  these  singularities  are  shown  in  Figure  1.  The  dashed 
curves  are  the  velocities  without  including  the  sources  starting 
at  the  body  base  and  the  solid  curves  include  these  singular¬ 
ities.  As  can  be  seen,  when  the  singularities  are  included  the 
velocities  return  to  the  free-stream  values. 


Similar  results  for  an  angle  of  attack  of  5°  are  shown  in 
Figure  2.  Here,  the  dashed  curves  have  only  source  singularities 
starting  at  the  body  base.  The  solid  curves  have  doublet  singu- 
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larities  also  starting  there.  With  the  doublets,  the  velocities 
again  return  to  the  free-stream  values  downstream  of  the  body  base 

3,2.2  Noncircular  Fuselage 

The  method  used  in  the  present  work  to  model  a  noncircular 
fuselage  is  an  adaptation  of  program  WDYBDY  used  in  the  work 
of  Reference  7.  This  program  was  extracted  from  the  advanced 
paneling  program  described  in  References  8  and  9.  The  method 
is  used  to  model  both  the  noncircular  fuselage  and  the 
elliptic  store  body  which  will  be  described  in  Section  3.2.6. 

The  theoretical  approach  is  based  on  representing  the  non¬ 
circular  fuselage  body  alone  by  supersonic,  linear  flow  theory. 
The  body  is  modeled  using  quadrilateral  supersonic  body  source 
panels  which  can  be  inclined  to  the  flow.  In  this  way  the 
panels  account  for  both  volume  and  angle  of  attack  effects. 

A  typical  panel  layout  on  an  elliptic  body  is  shown  in  Figure  3. 
The  inclination  of  a  panel  with  respect  to  the  body  axis  is 
limited  to  the  semiapex  angle  of  the  Mach  cone  associated  with 
the  free-stream  Mach  number.  This  limits  the  local  body  surface 
slope  which  can  be  modeled.  An  approximate  method  for  modeling 
engine  inlets  where  the  inclination  angle  is  greater  than  the 
Mach  cone  semiapex  angle  is  described  in  the  next  section. 

In  the  noncircular  fuselage  model  a  vertical  plane  of  symmetry 
is  assumed,  the  cross-sectional  shape  is  arbitrary  and  may 
vary  over  the  fuselage  length,  and  the  sideslip  angle  is  zero. 

The  panel  strengths  are  determined  by  solving  the  following 
set  of  simultaneous  algebraic  equations  for  N  panels  for  the 
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A ,  •  y  = 
13  T3 
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i  =  1, 2, .  .  .N 


(17) 


The  quantity  is  the  aerodynamic  influence  coefficient  of 
panel  i  on  panel  j,  Yj  is  the  unknown  strength  of  panel  j,  and 
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<WN/VJ  .  is  the  component  of  the  free-stream  velocity  normal 
of  the  panel  control  point,  its  centroid  of  area. 

The  calculation  of  the  A^j  matrix  is  described  in 
Reference  8.  All  the  fuselage  panels  have  unswept  leading  and 
trailing  edges  so  that  Equations  (88)  through  (95)  of  that 
reference  give  the  basic  velocity  expressions  for  one  corner 
of  the  panel  used  in  calculating  the  A^j's.  The  superposition 
of  the  four  corners  of  the  panel  is  given  by  Equations  (21) , 

(22),  and  (23)  of  Reference  8.  Equation  (13)  of  Reference  7 
is  the  expression  used  for  (Wfj/Vro)^. 

3.2.3  Fuselage  Mounted  Engine  Inlets 

This  section  is  concerned  with  the  flow  modeling  of  super¬ 
sonic,  ramp  type  engine  air  inlets  mounted  to  the  sides  of  the 
parent  aircraft  fuselage.  The  primary  objective  is  to  determine 
the  effects  of  the  inlets  on  the  flow  fields  in  which  the 
separated  store  may  be  immersed.  Effects  of  the  inlets  are 
also  included  in  accounting  for  fuselage/wing  interference. 

The  inlet  is  assumed  to  operate  at  on-design  condition  with  or 
without  flow  blockage.  In  either  case,  the  contribution  to 
the  flow  field  from  the  inlet  is  characterized  by  the  presence 
of  a  shock  system.  The  method  described  below  is  an  approxi¬ 
mate  one  and  is  based  on  the  supersonic  source  panels  used  for 
modeling  the  noncircular  fuselage  of  the  parent  aircraft  and 
the  store  body  if  the  elliptical  store  option  is  used. 

First,  a  description  is  given  of  the  inlet  of  interest 
and  the  flow  conditions.  Then,  the  technical  approach  is 
discussed.  The  method  used  to  generate  the  shock  system 
external  to  the  supersonic  inlet  will  be  discussed  in  Section 
4.2. 

The  configuration  treated  here  is  a  fuselage  mounted, 
cwo-dimensional,  external  compression  type  engine  air  inlet 
as  shown  in  the  following  sketch.  For  this  type  of  inlet,  the 
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leading  edges  of  the  sideplates  are  swept  back  at  angle 
90°  -  ^SE*  The  sharp  unswept  leading  edge  of  the  compression 
ramp  is  at  angle  0R  which  extends  into  the  inlet  duct,  and 
the  opposite  edge  forms  the  cowl  with  a  rounded  or  sharp  edge 
or  lip. 

The  oncoming  flow  is  assumed  at  free-stream  Mach  number, 

M^,  and  can  be  inclined  with  respect  to  the  duct  center  line. 

The  inlet  is  shown  "on-design"  and  operates  at  "critical  flow" 
condition.  The  first  designation  indicates  that  the  oblique 
shock  emanating  from  the  ramp  leading  edge  falls  on  the  intake 
opening.  Angle  0S  is  the  oblique  shock  angle  corresponding  to 
ramp  angle  0R  and  free-stream  Mach  number  M^.  For  the  on- 
design  case,  0_  =  For  critical  flow,  a  normal  shock 

stands  in  the  throat.  If  the  oncoming  Mach  number  is  below 
design,  the  oblique  shock  still  originates  from  the  ramp  edge 
but  stands  ahead  of  the  cowl  lip.  When  the  inlet  sees  higher 
than  design  Mach  number,  the  oblique  shock  may  be  turned  inside 
the  duct.  In  this  effort,  only  the  on-design  condition  will  be 
considered.  An  on-design  inlet  can  operate  subcritically , 
critically,  or  supercritically .  When  the  mass  flow  rate  in  the 
duct  is  less  than  that  contained  in  the  stream  tube  ahead  of  the 
inlet,  some  flow  spillage  occurs.  Under  this  partial  blockage 
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condition,  the  normal  shock  moves  ahead  of  the  throat  and  the 
inlet  operates  at  subcritical  condition.  This  situation  is 
depicted  in  the  following  sketch.  The  prediction  method  i3 
aimed  at  being  able  to  handle  the  zero  and  partial  blockage 
case. 

throat 


Length  A  is  the  height  of  the  stream  tube  containing  the 
c 


Length  AQ  is 


maximum  mass  flow  which  can  enter  the  inlet, 
the  height  of  the  stream  tube  actually  entering  the  inlet 
under  a  partial  blockage  condition.  The  mass  flow  ratio  MFR 
is  defined  as 

MFR  =  (18) 


The  norma]  and  oblique  shocks  combine  to  form  the  indicated 
coalesced  shock,  and  the  stream  tube  containing  the  spilled 
flow  is  deflected  from  the  intake. 
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The  descriptions  of  the  inlet  and  its  flow  features 
given  above  are  by  no  means  complete.  They  are  intended  only 
to  provide  the  nomenclature  and  parameters  used  in  the  predic¬ 
tion  method  discussed  below.  More  detailed  description  of 
supersonic,  external  compression  inlet  Characteristics  are 
available,  for  example  refer  to  Reference  10. 

The  method  described  here  for  modeling  the  inlets  makes 
use  of  the  supersonic  source  panels  already  employed  to  model 
the  noncircular  fuselage.  The  slant  area  of  the  intake  opening 
is  covered  by  the  source  panels.  These  inlet  panels  have  a 
special  boundary  condition  treatment  to  account  for  flow  block¬ 
age.  The  steps  in  the  procedure  are  as  follows.  The  first 
step  concerns  the  solution  for  the  panel  strengths.  Then,  the 
shock  shape  is  determined.  Finally,  the  contributions  to  the 
flow  field  are  calculated.  These  latter  two  are  discussed  in 
Sections  4.2  and  5.4. 

The  fuselage  source  panels  and  the  inlet  panels  are 
treated  integrally  by  one  set  of  simultaneous  equations, 
accounting  for  a  vertical  plane  of  symmetry,  to  solve  for  their 
strengths  subject  to  the  boundary  conditions  to  be  subsequently 
discussed.  In  this  process,  the  Mach  cones  determining  the 
region  of  influence  of  the  inlet  panels  require  special  treat¬ 
ment.  This  is  necessary  because  the  angles  of  incidence  of 
the  inlet  panels  are  such  that  the  panel  planes  usually  lie 
outside  the  Mach  cones  from  their  corners  based  on  the  free- 
stream  Mach  number.  In  fact,  they  are  superinclined  panels. 

The  present  source  panel  solution  cannot  handle  such  a  case. 

To  alleviate  this  problem  the  oblique  shock-wave  angle  is 
employed  to  specify  an  effective  Mach  number  for  use  in  the 
calculation  of  the  contributions  from  the  inlet  panels  such 
that  the  panels  lie  behind  their  Mach  cones.  The  Prandtl- 
Glauert  factor  0  for  the  source  panels  on  the  inlet  face  is 
calculated  as  follows.  For  given  free-stream  Mach  number,  Mm, 
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and  specified  ramp  angle  0_,  the  oblique  shock  angle  0 „  is 

K  S 

obtained  from  Reference  11.  For  the  Mach  wave  triangle  we  have 


6^  =  tan  (90°  -  v)  =  -  1  (19) 

In  an  analogous  manner,  replacing  the  Mach  wave  with  the 
oblique  shock  leads  to 

=  tan  (90°  -  6s)  (20) 

In  effect,  the  on-design  oblique  shock  is  used  to  position  a 
"Mach  wave"  slightly  ahead  of  the  inlet  panels.  It  is  this 
modified  "Mach  wave"  that  determines  the  region  of  influence  of 
the  inlet  panels.  Above  all,  this  scheme  allows  the  source 
panels  placed  on  the  inlet  panels  to  have  a  solution;  that  is, 
generate  perturbation  velocities  at  other  panel  control  points. 

The  original  flow  tangency  boundary  condition  associated 
with  the  source  panels  used  for  modeling  a  body  in  supersonic 
flow  is  described  on  page  52  of  Reference  8.  A  modified 
version  allowing  for  pitch  and  roll  is  discussed  in  Reference  7, 
Equation  (13) .  However,  the  boundary  condition  applied  at  con 
trol  points  of  the  source  panels  on  the  inlet  requires  special 
treatment  to  reflect  specified  mass  flow  ratio,  MFR,  Equation 
(18)  . 
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In  the  source  panel  solution,  the  flow  tangency  condition 
is  expressed  in  terms  of  velocity,  not  mass  flow,  normal  to 
the  plane  of  the  source  panels.  It  is  clear,  however,  that 
fully  blocked  panels  have  zero  mass  flow  and  zero  normal 
velocity  at  their  control  points.  Open  panels  allow  mass  flow 
at  free-stream  velocity  to  pass  through  at  their  control  points. 
Thus,  in  these  two  limits,  the  boundary  condition  can  be 
expressed  in  either  velocity  or  mass  flow  normal  to  the  panel. 

Refer  now  to  the  first  two  sketches  in  this  section.  In 
the  first,  all  of  the  mass  flow  contained  in  the  stream  tube 
of  height  A  enters  the  inlet  and  the  mass  flow  ratio,  MFR, 
equals  1.  The  boundary  condition  assigned  to  the  inlet  source 
panels  for  this  case  states  that  the  perturbation  velocities 
induced  by  all  source  panels  must  cancel  one  another.  To 
accomplish  this,  the  right-hand  side  of  the  inlet  source 
panel  flow  tangency  condition  is  set  equal  to  zero.  Equation 
(17),  Equation  (138)  in  Reference  8,  is  modified  for  the  full 
flow  situation  to 


0,  i  =  index  of  source  panels  on  inlet 
N  =  total  number  of  source  panels 
on  fuselage  and  inlet 


Here,  represents  the  aerodynamic  influence  coefficients 

and  Yj  are  the  source  panel  strengths.  Consequently,  all  the 
source  panels  modeling  the  fuselage  and  inlets  cause  the 
resultant  velocity  at  the  control  point  of  the  inlet  panels  to 
be  equal  to  that  of  the  free  stream.  Note  that  with  the 
vertical  plane  of  symmetry  assumption,  the  effects  of  the 
"image**  inlet  is  included  in  all  of  the  above.  Any  effects 
of  the  wing  thickness  or  lifting  panels  on  the  inlet  velocity 
are  neglected. 


The  second  sketch  in  this  section  depicts  the  partial 
blockage  case.  Now,  only  the  stream  tube  with  height  AQ 
enters  the  inlet.  The  remaining  stream  tube  with  height 
(A  -  A  )  is  spilled  from  the  intake.  For  this  case,  the 

L*  v 

region  designated  "OPEN"  on  the  intake  opening  is  made  to 
act  as  if  there  is  no  blockage.  The  panels  in  this  region  have 
the  boundary  condition  given  by  Equation  (21)  with  changed 
limits  on  the  index  i. 


A  *1^3  ■  °' 


i  =  index  of  OPEN  source  panels  on 
inlet 


(22) 


Now,  the  source  panels  in  the  region  designated  "BLOCKED"  are 
made  to  block  the  flow  and  actually  become  a  continuation  of 
the  upper  stream  surface  forming  the  stream  tube  containing 
the  spilled  flow.  The  flow  tangency  condition  is  applied  to 
these  panels. 


N 

l 

j=l 


AijYi  = 


W. 


N 


i  *  index  of  BLOCKED  source 

panels  on  inlet 

Q 


(23) 


Quantity  is  given  by  Equation  (13)  in  Reference  7  and 

represents  the  contribution  from  the  free  stream.  Regions 
OPEN  and  BLOCKED  are  determined  by  the  intersections  of  the 
stream  tubes  and  the  plane  of  the  inlet.  As  such,  the  lengths 
of  the  two  regions  are  proportional  to  the  heights  of  the  stream 
tubes,  Aq  and  (Ac  -  AQ) ,  respectively.  Thus,  the  *s»tio  of  the 
inlet  area  covered  by  the  OPEN  source  panels  to  the  total 
inlet  area  reflects  the  mass  flow  ratio  MFR,  Equation  (18) . 

To  solve  for  the  source  panel  strengths  modeling  a  non¬ 
circular  fuselage  with  engine  inlets,  the  set  of  equations 
gi^en  by  Equation  (17)  is  modified  in  accordance  with  Equation 
(21)  or  Equations  (22)  and  (23) .  In  calculating  the  aerodynamic 
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influence  coefficient  for  an  inlet  panel  acting  on  another 
panel,  given  by  Equation  (20)  is  used. 

3.2.4  Ejector  Rack  Body 

A  sketch  of  a  triple  ejection  rack,  TER,  wind  tunnel 
model  is  shown  in  Figure  4.  It  consists  of  the  rack  body  with 
two  stubby  pylons  at  the  shoulder  locations  and  a  long  pylon 
at  the  bottom  location.  This  same  configuration  exists  at  the 
two  longitudinal  stations  on  a  multiple  ejection  rack,  MER. 

In  the  present  work  the  pylons  are  not  modeled  and  the  rack 
body  is  modeled  as  a  body  of  revolution. 

The  rack  body  is  modeled  as  if  no  parent  aircraft  were 
present,  that  is  as  a  body  in  free-stream  flow  at  angle  of 
attack  only.  Volume  effects  are  modeled  by  a  distribution  of 
line  sources  along  the  body  axis  and  angle  of  attack  effects 
by  a  distribution  of  line  doublets  along  the  axis.  The  method 
of  calculating  the  singularity  strengths  and  the  perturbation 
velocity  field  is  identical  to  that  described  in  Section  3.2.1 
for  the  circular  fuselage. 

3.2.5  Circular  Store  Body 

The  circular  store  body  option  is  restricted  to  a  body 
of  revolution.  As  in  the  case  of  the  circular  fuselage  and  the 
ejector  rack  bodies,  the  circular  store  is  modeled  in  the  free- 
stream  at  angle  of  attack.  A  distribution  of  line  sources 
along  the  store  axis  models  the  volume  effects  and  a  distribution 
of  line  doublets  models  the  angle  of  attack  effects.  The  method 
of  calculating  the  singularity  strengths  and  the  perturbation 
velocity  field  is  identical  to  that  described  in  Section  3.2.1 
for  the  circular  fuselage. 

3.2.6  Elliptic  Store  Body 

The  elliptic  store  body  option  models  a  body  of  elliptic 
cross  section.  The  ratio  of  major  axis  length  to  minor  axis 
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can  vary  over  the  store  length.  This  option  can  also  be  used 
to  model  a  circular  store  which  is  to  be  separated.  In  the 
circular  store  option,  the  separated  store  can  have  only  one 
set  of  fins  and  they  must  be  planar  or  cruciform.  The 
elliptic  store  option  allows  the  separated  store  to  have  one 
or  two  sets  of  fins  with  1,  2,  3,  or  4  fins  in  each  set.  The 
fins  may  be  oriented  circumferentially  anywhere  on  the  body. 
Also,  they  do  not  have  to  be  perpendicular  to  the  tangent  to  the 
body  surface. 

The  elliptic  store  body  is  modeled  using  the  same  method 
as  was  described  in  Section  3.2.2  for  the  noncircular  fuselage. 
The  body  is  modeled  by  supersonic  source  panels  and  their 
strengths  are  calculated  using  Equation  (17). 

The  elliptic  store  model  does  not  assume  a  vertical  plane 
of  symmetry  as  was  done  for  the  noncircular  fuselage.  Source 
panels  are  laid  out  on  the  entire  body.  This  is  done  since, 
even  though  the  store  is  being  modeled  in  the  free- stream, 
it  can  be  rolled  as  well  as  pitched  relative  to  the  free-stream 
’• "^locity  vector.  An  example  of  this  is  an  elliptic  store  on 
one  of  the  shoulder  locations  on  a  TER. 

3.3  Flow  Model  for  Wing-Pylon  Combination  and 
Interference  Shell  on  Fuselage 

This  section  describes  the  representation  of  the  wing- 
pylon  in  terms  of  elements  from  linear,  potential  flow  t'  eory. 
'"so,  methods  used  to  account  for  mutual  interference 

between  the  wing-pylon  and  the  fuselage  are  discussed. 
Specifically,  constant  u-velocity  panels  are  distributed  over 
the  wing,  pylon,  and  part  of  the  fuselage.  In  addition,  con¬ 
stant  sc  ;e  panels  distributed  over  the  wing  and  pylon 
account  for  thickness  effects.  The  constant  u-velocity  panels 
laid  out  on  the  fuselage  form  the  interference  shell  and  serve 
primarily  to  cancel  interference  effects  from  the  wing-pylon 
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combination.  The  flow  model  is  an  extension  to  that  described 
in  Section  3.3  of  Reference  1  in  that  wing  dihedral  and  wing 
incidence  are  now  included. 

3.3.1  Geometrical  Layout 

The  left  wing  and  pylon  and  the  left  half  of  the  fuselage 
are  divided  up  into  area  panels  of  trapezoidal  shape.  The 
right  half  of  the  airframe  is  accounted  for  by  employing  sym¬ 
metry  properties  discussed  later.  Figure  5  shows  a  swept- 
wing,  swept-pylon  combination  attached  to  a  circular  fuselage 
above  its  midplane.  The  rectangular  coordinate  systems 
associated  with  the  fuselage  (xp,yp,zp)  and  the  wing  (xw,yw,zw) 
are  shown.  The  origin  of  the  former  is  at  the  fuselage  nose. 
The  wing  coordinate  system  has  its  origin  in  the  yp  *  0  plane 
at  the  point  where  it  is  intersected  by  the  line  connecting  the 
leading  edges  of  the  root  chords  of  the  left  and  right  wing 
panels.  The  coordinate  system  (x,y,z)  with  its  origin  at 
corner  1  of  the  cross-hatched  panel  is  associated  with  a  semi¬ 
infinite  triangle  which  is  part  of  the  solution  for  the  panel 
to  be  described  later.  The  panel  lies  in  the  x,y  plane. 

A  distribution  of  constant  u-velocity  type  panels  to  model 
loading  and  a  distribution  of  constant  source  type  panels  to 
model  thickness  will  be  laid  out  on  the  wing  and  pylon.  The 
distributions  of  the  two  types  of  panels  need  not  be  the  same 
in  the  chordwise  direction  or  in  the  spanwise  direction.  The 
sides  of  both  types  of  panels  are  taken  parallel  to  the 
centerline  of  the  fuselage.  Their  leading  and  trailing  edges 
are  swept  along  constant  percent  chordlines  of  the  wing  or 
pylon.  A  leading  or  trailing  edge  is  called  subsonic  or  super¬ 
sonic  depending  on  whether  or  not  the  component  of  the  free- 
stream  velocity  perpendicular  to  the  edge  in  question  is 
subsonic  or  supersonic,  respectively.  Figure  5  shows  a 
simplified  layout  of  2  chordwise  by  6  spanwise  panels  on  the 
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left  wing  half.  A  distribution  of  2  chordwise  by  2  spanwise 
panels  covers  the  pylon.  Generally,  more  panels  are  required 
to  obtain  accurate  results. 

Constant  u-velocity  panels  only  are  laid  out  over  part 
of  the  fuselage  and  designated  body  interference  panels.  The 
leading  and  trailing  edges  of  these  panels  are  unswept.  Only 
the  left  half  of  the  aircraft  configuration  needs  to  be  covered 
because  of  symmetry  about  the  yF  =  0  plane  provided  that  the 
effects  of  the  other  half  are  accounted  for  as  discussed  later. 
The  length  of  fuselage  covered  by  the  body  interference  panels 
is  determined  by  the  range  of  influence  of  the  wing  and  pylon. 
This  range  is  taken  as  the  length  along  the  body  measured  from 
the  leading  edge  of  the  wing-fuselage  junction,  or  root  chord, 
to  the  trailing  edge  of  the  wing  for  a  supersonic  wing 
trailing  edge.  For  the  case  of  a  subsonic  trailing  edge,  the 
length  should  extend  back  further.  In  this  case,  the  inter¬ 
ference  shell  is  sufficient  to  account  for  most  of  the  wing- 
pylon  to  fuselage  interference.  In  the  simplified  layout  of 
Figure  5,  the  body  interference  shell  consists  of  5  circum¬ 
ferential  rings  with  4  panels  per  half  ring.  The  interference 
shell  shape  is  constant  over  its  length.  For  a  noncircular 
fuselage  it  must  be  on  or  outside  the  source  panel  representa¬ 
tion  of  the  body,  see  Section  3.2.2,  and  is  not  circular  as 
shown  in  Figure  5. 

Each  constant  u-velocity  panel  contains  a  control  point 
placed  at  the  95-percent  location  on  the  chord  containing  the 
panel  centroid  as  shown  in  Figure  5.  The  flow  tangency  con¬ 
dition  is  applied  at  these  points. 

3.3.2  Constant  u-Velocity  and  Constant  Source  Panel  Solutions 

An  explanation  of  the  basic  solutions  for  the  two  types 
of  panels  will  not  be  presented  in  this  report.  They  are 
discussed  in  detail  in  Section  3.3.2  and  Appendix  II  of 
Reference  1. 
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With  the  induced  velocity  expressions,  or  solutions, 
for  the  two  types  of  panels  known,  it  is  possible  to  formulate 
the  influence  of  a  constant  u-velocity  panel  of  unknown 
strength  and  the  influence  of  a  constant  source  panel  with 
known  strength  at  any  control  point  accounting  for  the  region 
of  influence  inherent  in  supersonic  flight.  By  considering 
interactions  between  all  constant  u-velocity  panels  laid  out 
on  the  wing,  pylon,  and  body  interference  shell,  there  results 
a  set  of  simultaneous  equations  from  which  the  unknown  constant 
t  u-velocity  panel  strengths  can  be  determined.  This  formula¬ 
tion  will  now  be  discussed. 

3.3.3  Flow  Tangency  Boundary  Condition 

The  flow  tangency  boundary  condition  states  that  there  is 
no  flow  through  the  constant  u-velocity  panels  on  the  wing, 
pylon,  and  interference  shell  at  the  control  points.  A  few 
control  points  are  shown  in  Figure  5.  The  tangency  condition 
for  a  wing  panel  with  dihedral,  4>v,  is  shown  in  Figure  6  and 
can  be  written  in  the  wing  coordinate  system  shown  in  Figure 
5  as 


NPTOT  WWM  VW 

I  Y"  cos  ♦  -  sin  4>v  =  (ctp  +  iw  +  )cos  ^ 

n=l  oo  oo  v  \w 


+  JSLy.  sin  - 
00 


^i,  V  +  wWi,v 


COS  <f) 


V  =  1,2,... NPANLS 


In  the  above  equation  NPTOT  is  the  total  number  of  u-velocity 
panels  on  the  wing,  pylon,  and  fuselage  interference  shell  and 
NPANLS  is  the  number  on  the  wing. 


The  free-stream  components  are  the  fuselage  angle  of 
attack,  otp,  and  the  wing  incidence  angle,  iw*  Wing  incidence 
is  only  included  in  the  boundary  condition  in  the  computer 
program.  The  wing  is  not  placed  at  incidence  geometrically 
relative  to  the  fuselage  interference  shell  since  Mach  cones 
on  the  shell  and  the  wing  must  be  parallel  to  get  a  valid 
solution.  This  restricts  the  incidence  to  small  angles. 

The  angle  is  the  local  angle  of  attack  due  to  wing  twist 
and/or  camber  at  the  Vth  control  point  as  shown  in  the  following 
sketch.  The  velocities  u^  v,  vwi  y ,  ww^  v  are  interference 


velocities  induced  at  the  control  points  by  the  circular 

fuselage  line  sources  and  doublets  or  the  noncircular  fuselage 

source  panels,  wing  thickness  source  panels  when  there  is 

wing  dihed  al,  and  pylon  thickness  source  panels.  The 

velocities  vw  ,  ww  on  the  left-hand  side  of  the  equation 
v  f  rv  v  *  n 

are  perturbation  velocities  induced  at  the  control  point  by 
the  constant  u-velocity  panels  distributed  over  the  wing, 
pylon,  and  fuselage  interference  shell. 

The  pylon  boundary  condition,  also  shown  in  Figure  6,  is 
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(25) 


V  =  NPANLS+1,  NPANLS+2, . . .NPANLS+MP 
in  the  wing  coordinate  system.  The  velocities  Vw  are 

v  f  XI 

perturbation  velocities  induced  at  the  control  point  by  the 
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u-velocity  panels.  The  velocity  vw^  is  the  sum  of  the  inter¬ 
ference  velocities  due  to  the  circular  fuselage  Line  sources 
and  doublets  or  the  noncircular  fuselage  source  panels  and 
wing  thickness  source  panels.  The  index  MP  is  the  number  of 
u-velocity  panels  on  the  pylon. 


Similarly,  letting  vm  be  the  perturbation  velocity  in 

V  p  ii 

the  direction  normal  to  the  body  interference  panel,  the 
flow  tangency  condition  states  that  the  net  velocity  normal 
to  the  body  interference  panel  is  zero. 


NPTOT  VN  „ 
V  v  >  n 

L  V 

n=l  00 


(26) 


V  *  NPANLS+MP+1 , NPANLS+MP+2 , . . .NPTOT 
U  =  V-NPANLS-MP 

The  right-hand  side  of  the  above  equation  represents  the 
external 1  induced  perturbation  velocity  normal  to  the  body 
interference  panel  under  consideration.  It  contains  inter¬ 
ference  velocities  induced  by  source  panels  laid  out  on  the 
wing  and  the  pylon  to  model  thickness.  In  this  way,  inter¬ 
ference  on  the  fuselage  due  to  wing-pylon  thickness  is 
accounted  for.  Angle  $  is  called  the  body  panel  orientation 
angle  measured  relative  to  a  plane  parallel  to  the  (x_,,y_,)  or 

r  r 

(xw,yw)  plane  as  indicated  in  the  following  sketch  looking  in 
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the  positive  x„  direction.  The  angle  <fc  shown  in  the  sketch 

*  y 

is  a  negative  angle.  Finally,  the  summation  on  the  left-hand 
side  of  Equation  (26)  represents  the  perturbation  velocity 
normal  to  the  interference  panel  under  consideration  induced 
by  all  constant  u-velocity  panels.  So  far,  the  strengths  of 
these  panels  are  still  unknown. 

In  the  boundary  conditions  formulated  above,  advantage 
is  taken  of  the  fact  that  the  yF  =  0  or  yw  =  0  plane  is  a 
plane  of  symmetry.  Panels  are  laid  out  over  the  left  wing  and 
pylon  and  only  the  left  half  of  the  fuselage  is  covered  with 
interference  panels.  Perturbation  velocities  on  both  sides 
of  Equations  (24),  (25),  and  (26)  are  induced  not  only  by  the 
panels  to  the  left  of  the  symmetry  plane  but  also  must  con¬ 
tain  contributions  from  the  right  half.  The  required  procedure 
is  discussed  below  in  connection  with  the  aerodynamic  influence 
coefficients. 

Next,  it  will  be  shown  how  the  perturbation  velocities 
are  expressed  in  terms  of  aerodynamic  influence  coefficients. 
These  coefficients  relate  the  perturbation  velocity  components 
induced  at  some  point  by  a  constant  u-velocity  panel  to  its 
strength  u^/V^  and  the  coordinates  of  the  point  relative  to 
the  panel  corners.  A  similar  procedure  relates  the  perturba¬ 
tion  velocities  at  a  field  point  to  the  known  strength  of  a 
source  panel. 

3.3.4  Aerodynamic  Influence  Coefficients 

Appendix  II  of  Reference  1  contains  the  results  of  the 
theory  for  determining  the  perturbation  velocities  induced  by 
a  semi-infinite  triangle  subject  to  a  constant  jump  condition 
in  axial  velocity.  The  appendix  also  specifies  the  perturba¬ 
tion  velocities  induced  by  a  semi-infinite  triangle  with  a 
constant  strength  source  distribution.  The  results  are  given 
below  in  symbolic  form.  In  terms  of  the  coordinate  system 
associated  with  the  trangle  shown  in  Figure  7,  the  perturbation 


velocities  due  to  a  triangle  with  sweptback  leading  edge  and 
with  constant  jump  in  u-velocity  are  expressed  as  follows. 
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The  quantity  u+/Voo  is  the  constant  strength  of  the  singularity 
distributed  over  the  triangle.  In  fact,  it  is  proportional  to 
the  axial  velocity  u+  in  the  plane  of  the  triangle  (z  =  0) , 

The  computer  program  has  been  arranged  to  predict  axial 
velocity  u+  in  the  zero  plane  and  requires  a  very  small  nega¬ 
tive  value  for  z  to  produce  axial  velocity  u_  shown  in  Figure  7. 
The  complete  expressions  for  the  perturbation  velocities  are 

given  in  Appendix  II  of  Reference  1.  Functions  F  ,F  ,F  are 

U  V  w 

called  influence  functions  and  depend  on  coordinates  x,y,z  of 
the  point  at  which  the  velocities  are  computed  relative  to  the 
apex  of  the  semi-infinite  triangle  shown  in  Figure  7.  The 
influence  functions  can  be  determined  from  Appendix  II  of 
Reference  1,  Equations  (II-4)  and  ( 11-12) ,  and  depend  on  the 
leading-edge  sweep  \|/  and  the  factor  6  given  in  terms  of  the 
free- stream  Mach  number  by 

0  =Vm2  -  1  (28) 

00 

Perturbation  velocities  induced  by  the  semi-infinite 
triangle  with  constant  source  strength  are  related  to  a  surface 
tangent,  tan  0,  to  be  discussed  later,  and  the  coordinates  of 
the  field  point  relative  to  the  triangle  apex.  Employing  the 
subscript  t  to  indicate  association  with  thickness,  the  per¬ 
turbation  velocities  are  expressed  as 
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Functions  FUt'Fvt'Fwt  are  the  influence  functions  associated 
with  the  thickness  and  depend  on  coordinates  x,y,z  of  the  field 
point  relative  to  the  apex  of  the  semi-infinite  triangle  such 
as  the  one  with  sweptback  leading  edge  shown  in  Figure  7. 

These  functions  can  be  determined  from  Appendix  II  of 
Reference  1,  Equations  (11-15)  and  (11-16),  and  depend  on  the 
leading-edge  sweep  \pT„  and  the  factor  &  given  by  Equation  (28). 
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In  order  to  formulate  the  perturbation  velocities  induced 
by  a  trapezoidal  constant  u-velocity  panel,  the  superposition 
principle  described  in  Section  3.3.2  of  Reference  1  and 
depicted  in  Figure  4  of  that  reference  must  be  employed.  For 
example,  for  a  panel  with  sweptback  leading  and  trailing  edges 
on  the  left  wing,  the  upwash  at  the  Vth  control  point  induced 
by  the  n  constant  u-velocity  panel  in  the  coordinate  system 
of  the  n  panel  is  expressed  as 
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The  subscripted  coordinates  represent  the  coordinates  of  the 

control  point  relative  to  one  of  the  four  panel  corner  points. 

The  numbering  sequence  of  the  corner  points  is  kept  consistent 

in  its  relation  to  the  directions  of  the  coordinate  systems 

shown  in  Figure  4  of  Reference  1.  The  leading-edge,  and 

trailing-edge  sweep  angles  of  the  panel  are  designated  \pr„ 

and  respectively.  Here  function  Fw  is  defined  in  Equation 

(27) .  Similar  expressions  can  be  written  for  u  and  v 

v ,  n  v ,  n 

using  Fu  and  F  ,  respectively.  The  left-hand  sides  of 
Equations  (24),  (25),  and  (26)  require  the  velocities  in  the 
wing  coordinate  system.  The  transformation  is 
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In  a  similar  manner,  the  perturbation  velocities  induced 
by  a  trapezoidal  constant  source  panel  are  obtained  by  means  of 
a  superposition  scheme  using  the  solutions  for  four  semi¬ 
infinite  triangles  with  their  apexes  at  each  of  the  panel 
corners.  The  source  strength  for  one  panel  is  given  by  the 
surface  slope  of  the  wing  or  pylon  thickness  distribution  at  the 
panel  centroid.  It  is  expressed  as  the  tangent  of  the  angle  9 
associated  with  the  thickness  envelope  as  shown  in  the  following 
sketch.  Therefore,  the  upwash  due  to  thickness  at  the  v 
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thickness  envelope 


source  panel  centroid 


mean  camber  surface 
made  to  be  flat 


tangent  to  thickness  envelope 
at  axial  location  of  the 
centroid 


control  point  induced  by  the  constant  source  panel  with  index 
nt  on  the  left  wing  is  expressed  in  the  coordinate  system  of 

j.  t 

the  nt  panel  as 
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(32) 


Here  pertains  to  the  leading-edge  or  trailing-edge  sweep  of 
the  constant  source  panel  accounting  for  thickness.  Similar 
expressions  can  be  written  for  ut4l  .  and  vt ,  ..in  Equation 

v  f nt  v r nt 

(29)  using  FUt  and  FVt«  The  thickness  panel  induced  velocities 

make  up  part  of  the  u^  ,  vwi  ,  and  wwi  v  terms  on  the  right- 

hand  sides  of  Equations  (24),  (25),  and  (26)  in  the  wing 

coordinate  system.  The  Bame  transformation  as  was  used  in 

Equation  (31)  can  be  used  to  put  ut,  .,  vt,,  ^  and  wt„  . 

uv,nt  v,nt  ^v,nt 

into  that  system. 

The  bracketed  terms  in  Equations  (30)  and  (32)  are  the 
upwash  aerodynamic  influence  coefficients  associated  with  a 
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constant  u-velocity  panel  and  a  constant  source  panel, 
respectively. 


3.3.5  Symmetry  Considerations 

In  order  to  account  for  the  effects  of  the  constant 
u-velocity  and  wing  thickness  source  panels  on  the  right  half 
of  the  aircraft  configuration,  the  perturbation  velocity  com¬ 
ponents  induced  by  a  panel  on  the  left  are  augmented  by  a 
contribution  from  its  image  panel  on  the  right.  However, 
the  following  simplification  is  used.  The  effect  of  the  image 
panel  at  a  given  point  can  be  obtained  using  the  panel  on  the 
left  with  a  change  in  sign  of  the  yw  or  yp  coordinate  of  the 
point.  The  sidewash  calculated  this  way  at  (xw,-yw, zw)  is 
then  transformed  back  by  a  change  in  sign.  The  procedure  is 
shown  in  the  following  sketch  for  a  panel  on  the  left  wing  half 
of  a  wing-alone  configuration.  The  coordinate  system  (x,y,z) 
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associated  with  the  semi- infinite  triangle  with  its  apex  at 
corner  1  will  be  used.  The  panel  lies  in  the  x,y  plane.  It 
is  the  same  as  the  coordinate  system  of  Figure  7.  The  coordi¬ 
nates  of  point  P  relative  to  corner  1  (xw1,yw1,z^1)  are  then 
given  in  the  wing  system  by 


A*w 
-  4*w 


*w  -  XW1 


(33) 


Let  the  panel  be  at  dihedral  angle,  <|u  ,  with  respect  to  the 
yw  axis.  Then  the  location  of  point  P  relative  to  corner  1 
in  the  panel  (x,y,z)  system  is 


x  =  -  Axw 

y  =  Aywcos  4>i  +  Azwsin  ) 

z  =  Aywsin  <J»i  -  Azwcos  4>i 


(34) 


With  these  local  coordinates,  the  direct  influence  of  corner  1 
at  point  P  can  be  calculated.  The  influences  of  the  other 
corners  of  the  panel  are  added  in  accordance  with  the  super¬ 
position  scheme  as  indicated  for  the  w  component  in  Equations 
(30)  and  (32). 

The  location  of  the  image  point,  P^,  in  the  wing  coordi¬ 
nate  system  is 
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and  Equation  (34)  can  again  be  used  to  locate  this  point  in 
the  panel  system.  The  influence  of  corner  1  is  now  calculated 
at  point  P^.  The  influences  of  the  other  corners  of  the  panel 
are  added  in  accordance  with  the  superposition  scheme.  The 
resulting  velocity  components  u,v,w  in  the  panel  system  repre¬ 
sent  the  influence  of  the  panel  on  the  left  wing-half  on  the 
point  P^.  With  a  change  in  the  t;:?gn  of  the  sidewash,  v,  they 
also  represent  the  influence  of  the  u'age  panel  on  the  point 
P.  Thus,  at  a  given  point  P  there  will  be  a  direct  influence 
induced  by  a  constant  u-velocity  or  constant  source  panel  and 
an  indirect  effect  induced  by  the  image  of  the  panel  under 
consideration. 

This  method  of  accounting  for  symmetry  is  also  applied 
to  constant  u-velocity  and  constant  source  panels  on  a  wing- 
mounted  pylon  and  the  constant  u-velocity  panels  on  the  fuselage 
interference  shell.  The  influence  functions  required  in  Equa¬ 
tions  (30)  and  (32)  and  defined  by  Equations  (27)  and  (29)  for 
these  panels  are  obtained  from  Appendix  II  of  Reference  1, 
Equations  (II-4) ,  (11-12),  and  Equations  (11-15),  (11-16), 
using  suitable  transformations.  For  the  pylon,  a  90°  counter¬ 
clockwise  rotation  about  the  x-axis  (viewing  upstream)  places 
the  semi-infinite  triangles  in  the  vertical  plane  and  the 
influence  functions  can  be  determined.  The  superposition 
scheme  gives  the  solution  for  a  trapezoidal  panel  on  the  pylon. 
Effects  of  the  image  wing-mounted  pylon  panel  are  obtained  in 
the  manner  described  above  for  a  panel  on  the  wing.  The  same 
procedure  is  applied  to  the  pylon  source  panels.  The  influence 
functions  associated  with  a  body  interference  panel  are  also 
obtained  after  performing  a  rotation  about  the  x-axis.  The 
angle  of  rotation  is  different  for  each  panel  and  is  related 
to  the  body  orientation  angle  <f>^  mentioned  earlier  in  connec¬ 
tion  with  the  boundary  conditions.  Image  panel  effects  are 
accounted  for  using  the  same  procedure  as  used  for  a  panel  on 
the  wing. 
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The  terms  on  both  sides  of  Equations  (24) ,  (25)  ,  and  (26) 
can  now  be  specified.  Equation  (30)  for  w/V^  and  a  similar 
expression  for  v/Vtt  are  used  on  the  left-hand  sides.  These 
expressions  give  velocities  in  the  coordinate  system  of  the 
nth  panel  while  Equations  (24) ,  (25)  ,  and  (26)  require 
velocities  in  tho  wing  coordinate  system.  The  velocity  trans¬ 
formation  '  .  c;.-r  panel  system  to  the  wing  system  is 
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(36) 


The  u^  ,  vwi  ,  and  wwi  v  terms  in  Equations  (24), 

(25) ,  and  (26)  are  externally  induced  perturbation  velocities. 
In  Equation  (24)  for  the  wing  they  are  the  sums  oi  the 
velocities  induced  at  the  v*"*1  control  point  by  the  circular 
fuselage  line  sources  and  doublets  or  the  noncircular  fuselage 
source  panels,  the  wing  thickness  source  panels,  and  the  pylon 
thickness  source  panels. 

The  circular  fuselage  induced  velocities  in  the  axial, 
radial,  and  tangential  directions  are  given  by  Equation  (16) . 
The  positive  directions  are  shown  in  the  sketch  following  the 
equation.  The  velocities  are  transformed  into  the  wing  coordi¬ 
nate  system  shown  in  Figure  5  using 
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The  noncircular  fuselage  induced  velocities  due  to  the 
source  panels  are  calculated  using  the  superposition  scheme 
given  by  Equations  (21)  through  (23)  of  Reference  8  and  the 
velocity  expressions  given  by  Equations  (88)  through  (95) 
which  are  multiplied  by  the  panel  strength.  These  velocities 
are  in  the  panel  coordinate  system.  The  program  extracted  from 
Reference  7  rotates  them  into  the  body  coordinate  system  shown 
in  Figure  3  of  this  report.  The  signs  on  the  u  and  w  velocity 
components  are  changed  to  put  them  in  the  wing  or  fuselage 
coordinate  system  shown  in  Figure  5. 

Wing  and  pylon  source  panel  effects  are  given  by  Equation 
(32)  for  wty  nt«  Similar  expressions  can  be  written  for  u^v  nt 
and  vtv  nt  using  FUt  and  FVt»  These  are  velocities  induced 
at  the  control  point  by  the  nt*"*1  panel  in  the  coordinate 
system  of  the  nt  panel.  They  can  be  resolved  into  the  wing 
system  by 
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The  total  perturbation  velocity  induced  at  the  v  u-velocity 
panel  control  point  by  all  of  th*  wing  and  pylon  thickness  source 
panels  are  obtained  by  summing  Equation  (38)  over  all  of  the 
source  panels. 

The  vw^  ^  velocities  appearing  on  the  right-hand  side 
of  Equation  (25) ,  the  pylon  u-velocity  boundary  condition, 
are  the  sums  of  those  induced  by  the  circular  fuselage  line 
sources  and  doublets  or  noncircular  fuselage  source  panels  and 
the  wing  thickness  source  panels.  The  vM.  and  w  . 
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velocities  appearing  in  the  fuselage  interference  shell  boundary 
condition.  Equation  (26) ,  are  the  sums  of  the  wing  and  pylon 
thickness  source  panel  velocities.  These  velocities  are  cal¬ 
culated  in  the  manner  just  described  for  the  wing. 

3.3.6  Solution  for  Constant  u- Velocity  Panel  Strengths 

After  recasting  the  perturbation  velocities  on  the  left 
side  of  Equations  (24),  (25),  and  (26)  in  terms  of  the  influence 
functions  specified  by  Equation  (27)  and  accounting  for  symmetry, 
there  results  a  set  of  simultaneous  equations  in  which  the 
unknowns  are  the  NPTOT  values  of  panel  strengths  u+/Vco.  The 
values  can  be  obtained  through  a  matrix  solution  for  given 
angle  of  attack,  ctp,  wing  incidence  angle,  iw,  and  Mach 

number ,  M  . 

9  00 

3.3.7  Superposition  Scheme  for  Flow-Field  Calculations 

In  the  preceeding  sections  the  solution  for  a  constant 
u-velocity  panel  or  thickness  constant  source  panel,  such  as 
the  cross-hatched  one  shown  in  Figure  5,  is  obtained  through 
a  superposition  scheme  using  four  semi-infinite  triangles  with 
their  apexes  at  each  of  the  four  corners  of  the  panel.  At  a 
given  field  point  relative  to  the  panel,  the  influence  func¬ 
tions  due  to  the  triangles  at  the  corners  of  a  panel  are 
calculated,  summed,  and  multiplied  by  the  panel  strength.  This 
is  repeated  sequentially  for  all  panels.  This  scheme  must  be 
used  when  determining  the  strengths  of  the  u-velocity  panels. 
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Once  the  constant  u-velocity  panel  strengths  are  known, 
a  different  superposition  scheme  can  be  used  which  eliminates 
many  redundant  influence  function  calculations.  Consider  the 
wing  shown  in  the  following  sketch.  Corner  5  is  associated 


with  four  panels  so  that  using  the  four  corner  superposition 
scheme  requires  calculating  the  influence  functions  for  that 
corner  four  times,  once  for  each  of  the  surrounding  panels. 

This  redundant  calculation  can  be  eliminated  if  the  wing  is 
represented  by  a  series  of  corners  and  a  net  singularity 
strength  calculated  for  that  corner.  In  calculating  velocities 
at  a  field  point,  the  influence  functions  for  the  corner  are 
calculated  once  and  multiplied  by  the  net  strength. 

The  following  sketch  shows  a  corner  k  surrounded  by  four 
panels  of  known  strength.  The  net  strength  of  the  u-velocity 


panel  singularity  associated  with  this  corner  is 
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and  for  corner  5  in  the  previous  sketch 
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For  corners  which  do  not  have  panels  in  some  of  the  quadrants 
surrounding  the  point,  the  strengths  in  those  quadrants  are 
treated  as  zero.  Thus,  for  corner  1,  for  example. 
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and  for  corner  6 


<v6  ■  <V0R 


-  (VUL-  <u+>2il  -  '“+>2i2 


The  wing  shown  in  the  first  sketch  in  this  section  has 
a  break  in  leading  edge  sweep.  Since  the  influence  functions 
are  functions  of  this  angle,  two  rows  of  corner  points  are 
defined  at  this  location,  one  associated  with  the  panels  to 
the  left  and  the  other  associated  with  the  panels  to  the 
right.  In  calculating  the  net  strengths  of  corners  7,  8,  and 
9,  the  strengths  of  panels  1,3  and  2,3  are  assumed  to  be  zero 
For  example, 
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Similarly,  for  corners  10,  11,  and  12,  the  strengths  of 
panels  1,2  and  2,2  are  assumed  zero.  For  corner  12 


This  same  procedure  is  followed  at  points  where  breaks  in 
trailing-edge  sweep  angle  or  wing  dihedral  angle  occur. 

The  procedure  just  described  for  the  wing  is  followed 
for  the  pylon  and  fuselage  u-velocity  panel  interference 
shell.  The  latter  is  treated  like  a  wing  with  a  change  in 
wing  dihedral  angle  at  each  panel  side  edge. 

The  superposition  scheme  just  described  for  the  constant 
u-velocity  panels  is  also  used  for  the  wing  and  pylon  thick¬ 
ness  source  panels.  A  net  source  panel  strength  is  calculated 
for  each  corner. 


4.  NONLINEAR  CORRECTION  TO  FLOW  MODELS 


The  flow  models  described  in  Section  3  for  the  various 
parent  aircraft  components  are  based  on  supersonic,  linear 
flow  theory.  Linear  theory  causes  disturbances  to  travel 
along  Mach  waves  determined  only  by  the  free- stream  Mach 
number  which  was  used  in  the  solution  for  the  singularity 
strengths.  In  actuality,  the  disturbances  travel  along 
curved  Mach  waves  whose  shapes  are  determined  by  the  local 
Mach  number  variation  along  the  waves.  If  straight  linear 
theory  is  used  in  calculating  the  flow  field  under  the  air¬ 
craft  in  the  region  in  which  a  separated  store  is  immersed, 
parent  aircraft  disturbances  will  strike  the  store  in  the 
wrong  location.  This  will  cause  errors  in  the  calculated 
forces  and  moments. 

In  order  to  place  the  disturbances,  particularly  shock 
waves,  close  to  their  correct  locations  a  series  of  nonlinear 
corrections  are  made  to  the  linear  theory  solutions.  They  are 
engineering  corrections  and  have  been  selected  because  they 
require  little  computer  time  and  seem  to  be  quite  accurate. 

4 . 1  Body  Shock  Wave  Shape 


4.1.1  Circular  Body 

The  procedure  used  to  determine  the  shock  wave  shape  for 
a  circular  or  axisymmetric  body  will  now  be  described.  This 
method  is  used  for  the  circular  fuselage,  the  ejector  rack,  and 
the  circular  stores  which  are  modeled  by  line  sources,  sinks, 
and  doublets.  Experimental  data  in  Reference  6  show  that  a 
good  approximation  to  the  shock  wave  shape  at  angle  of  attack 
is  to  calculate  it  for  the  body  at  zero  degrees  angle  of 
attack.  This  aligns  the  longitudinal  axis  of  the  axisymmetric 
shock  wave  with  the  free  stream.  The  data  show  that  at  five 
degrees  angle  of  attack  there  is  little  movement  of  the  shock 
relative  to  the  free  stream. 
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The  nonlinear  axi symmetric  sh^ck  wove  shape  is  calculated 
in  the  manner  which  will  now  be  described.  Consider  the 
following  sketch  which  shows  an  axisymmetric  body  and  the 
linear  theory  Mach  cone  "riginating  at  the  body  nose.  The 
source-sink  distribution  modeling  the  body  is  calculated  as 


described  in  Section  3.2.1.  The  nonlinear  shock  shape  is  also 
shown  in  the  sketch.  It  is  curved.  The  slope  of  the  shock 
at  any  radial  point  r  is  assumed  to  be  determined  by  the  stream- 
wise  component,  of  the  Mach  number  immediately  ’  ehind  the 
shock.  If  the  slope,  is  known  as  a  function  of  r,  then 

the  shape  of  the  shock  wave  can  be  calculated  from  the  following 
integral 


x 


r 

6^ (r)dr 
o 


(40) 


where  =  vM£  -  1.  At  the  tip  of  the  nose,  x  =  r  =  0,  B^ 
is  determined  from  the  shock  wave  angle  for  a  cone  with  the 


same  angle  as  that  at  the  tip  of  the  nose  of  the  body.  This 
shock  angle,  0,  can  be  determined  from  Chart  5  of  Reference  11 
Then 


(Bp)  =  tan  (90°  -  6)  (41) 

r=0 

The  values  of  B^  at  other  values  of  r  are  determined 
from  M^(r).  At  a  given  value  of  r,  is  calculated  in  a 
manner  similar  to  that  used  in  Section  3.4.1  of  Reference  1 
for  the  wing  leading-edge  shock  wave.  At  a  radial  distance 
from  the  body  axis,  for  example  r^  in  the  preceeding  sketch, 
the  velocity  field  induced  by  the  source-sink  distribution 
is  calculated  beginning  at  the  linear  theory  Mach  cone  location 
and  extending  a  short  distance  downstream.  From  this  velocity 
field,  the  maximum  change  in  flow  angle,  Av,  is  determined. 

It  can  be  calculated  using  the  source  terms  in  Equation  (16) 
and  is 


where 


tan  Av 


iyvJ 

U/Vro 


1  + 


(42) 


(43) 


The  velocities  and  are  the  axial  and  radial  perturbation 
velocities  which  are  positve  aft  and  outward,  respectively. 

The  change  in  flow  angle,  Av,  corresponds  to  a  compression 
of  the  free-stream  flow.  Assuming  that  this  compression  is 
isentropic,  the  Prandtl-Meyer  angle,  v^,  at  the  point  where 
the  maximum  Av  occurs  is 


=  v  -  Av 


(44) 
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Prandtl-Meyer  angle  is  associated  with  the  free-stream  Mach 
number,  Mro,  and  is  given  by  Equation  (171c)  of  Reference  11. 

v  =  2.4495  tan-1  ro.40825(M2  -  1) 1/2 1  -  tan'^M2  -  1)1/2  (45) 

00  L  00  J  00 


In  accordance  with  Reference  12,  the  Mach  number  associated 
with  Prandtl-Meyer  angle  can  be  determined  most  conveniently 
as  folxows: 

2  3  ^ 

1  +  1.3604  p  +  0.0962  -  0.5127 

M1  1  -  0.6722  p  -  0.3278  p2 


P 


v 


max 


2/3 


}  (46) 


v 

max 


1.4 


J 


The  author  of  Reference  12  claims  0.05-percent  accuracy  over 
the  entire  range  of  Mach  number  (1  <  <  °°)  .  The  component  of 

the  Mach  number,  M^,  in  the  axial  free-stream  direction,  M^, 
is  then  given  by 


=  Micos  (47) 

The  above  calculation  is  repeated  for  a  series  of  radial 
locations,  such  as  r^r^r...^  shown  in  the  preceeding  sketch, 
and  then  the  nonlinear  shock  wave  shape  is  calculated  using 
Equation  (40) . 

Shock  wave  shapes  for  an  ogive-cylinder  body  calculated 
using  the  method  just  described  are  shown  by  the  dashed  curves 
in  Figures  8  and  9.  Also  shown  are  the  linear  theory  Mach  cone 
and  the  shock  wave  shape  calculated  using  an  Euler  equation 
computer  program.  Reference  13.  The  present  method  agrees 
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well  with  the  Euler  equation  solution.  In  Figure  8  the  hori¬ 
zontal  bars  show  the  shock  wave  location  as  estimated  from  the 
experimental  flow-field  data  of  Reference  6.  Flow  field 
measurements  were  made  every  0.25  inch.  The  bars  represent 
the  two  bracketing  data  points.  The  present  method  shows  good 
agreement  with  the  data.  The  crosses  in  Figures  8  and  9  are 
shock  wave  shapes  calculated  using  the  source  panel  method. 
This  will  be  described  in  the  next  section. 


It  was  mentioned  previously  that  with  a  body  at  angle  of 
attack  there  is  little  movement  of  the  shock  wave  relative  to 
the  free-stream  when  compared  with  the  shock  wave  produced 
by  the  body  at  zero  degrees  angle  of  attack.  Since  the  shock 
wave  shape  calculated  in  the  manner  just  described  is  for  the 
body  at  zero  degrees  angle  of  attack,  a  transformation  is 
required  to  locate  a  point  on  the  shock  wave  in  the  body 
coordinate  system  with  the  body  at  angle  of  attack.  Let  xQ 
and  rQ  be  the  coordinates  of  a  point  on  the  shock  at  a  =  0°. 
The  coordinates  x^  and  rfl  in  the  body  system  at  angle  of 
attack  are 


x^  =  xQcos  (otcos  9)  +  rQsin  (otcos  0)  j 
r^  =  -xQsin  (otcos  6)  +  r0cos  (otcos  6)  J 

where  9  is  measured  a<?  shown  in  the  following  sketch 


(48) 


shock 
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The  cos  8  variation  is  used  so  that  correct  motion  or  the 
shock  relative  to  the  body  is  obtained. 

4.1.2  Noncircular  Body 

The  shock  wave  shape  for  a  noncircular  body,  noncircular 
fuselage  or  elliptic  store,  is  calculated  in  the  same  manner 
as  was  just  described  for  a  circular  body.  The  circular  body 
shock  wave  is  axisymmetric  so  the  shape  is  only  calculated 
in  one  meridian  plane.  For  a  noncircular  body  the  shape  is 
not  axisymmetric  so  it  is  calculated  in  a  number  of  meridian 
planes. 

The  shock  wave  shape  for  the  noncircular  fuselage  is 
calculated  from  the  velocities  induced  by  the  body  source  panels 
whose  strengths  are  determined  with  the  fuselage  at  angle  of 
attack.  In us,  a  transformation  similar  to  that  give  by 
Equation  (48)  is  not  required.  The  velocities  used  in  the  shock 
shape  calculation  are  in  the  wind  axis  system.  Let  these  be 
u' ,  v',  and  w' .  They  are  determined  from  those  calculated  in 
the  body  axis  system  of  Figure  3  by 


u'  =  ucos  a  +  wsin  a 


v '  =  v 


w'  =  -usin  a  +  wcos  a 


Then  Equations  (42)  and  (43)  are  replaced  by 


tan  Av  = 


/v '  2  +  w '  2/V 


where 


JL  =  i  +  Hi 


The  shock  wave  shape  for  the  elliptic  store  is  calculated 
with  the  store  at  a  =  0°,  that  is,  the  source  panel  strengths 
were  calculated  for  a  =  0°.  Thus,  the  velocities  in  the  wind 
axis  system  are  equal  to  those  in  the  body  system  of  Figure  3. 


u'  =  u 


v '  =  v 
w1  =  w 


The  shock  shape  is  calculated  using  the  method  of  Section  4.1.1 
with  Equation  (50)  and  (51)  replacing  Equations  (42)  and  (43) . 

Since  the  elliptic  store  shock  wave  shape  is  calculated 
for  a  =  0°,  a  transformation  is  required  to  locate  a  point  on 
the  store  shock  wave  in  the  body  coordinate  system  when  the 
store  is  at  angle  of  attack.  The  transformation  used  is 


x  =  x  cos  a'  +  r  sin  a' 
a  o  o 


where 


-x^sina*  +  r„cos  o' 
o  o 


a'  =  a(l  -  e  )cos  6 
a 


and  is  similar  to  Equation  (48) .  The  additional  factor  (1  -  ea) 
appearing  in  the  expression  for  a'  is  taken  from  Reference  14. 
The  parameter  eQ  defines  the  amount  of  deviation  from  a  full 
rotation  of  the  shock  about  store  nose  as  it  goes  to  angle  of 
attack.  This  factor  is  strictly  applicable  only  to  a  cone.  A 
plot  of  e  as  a  function  of  cone  semi-apex  angle,  6  ,  and  Mach 
number,  M^,  is  presented  in  Figure  10. 

The  accuracy  of  the  shock  wave  shape  calculated  for  a 
body  modeled  by  source  panels  is  shown  in  Figures  8  and  9.  The 
ogive-cylinder  body  was  modeled  by  source  panels.  The  shock 
wave  shapes  are  well  predicted. 


4.2  Inlet  Shock  Wave  Shape 


The  noncircular  fuselage  and  engine  inlets  are  modeled 
by  source  panels  as  described  in  Sections  3.2.3  and  3.2.4. 

Once  the  panel  strengths  are  known,  an  approximate  shape  of  the 
shock  system  associated  with  the  inlets  can  be  determined  in 
a  way  which  will  now  be  described.  It  will  be  assumed  that 
the  flow  ahead  of  the  inlets  is  at  free  stream  conditions  and 
that  the  shock  shape  depends  on  the  local  inlet  configuration 
details  only.  Inlet  induced  effects  are  scanned  along 
strategically  located  axial  traverses  for  the  maximum  down- 
wash  velocity  component.  The  approximate  shape  of  the  shock 
system  is  then  determined  from  an  integration  scheme.  The 
objective  is  to  generate  a  table  containing  coordinates  of  the 
shock  shape  for  the  inlet  flow  conditions  at  hand. 


Consider  first  the  zero-blockage  case,  i.e.,  the  inlets 
are  "flowing  full"  and  operate  on  design.  The  flow  field  is 
calculated  at  points  along  traverses  laid  out  parallel  to  the 
centerline  of  the  parent  aircraft.  A  set  of  planes  containing 
axial  traverses  are  shown  in  the  cross  flow  plane  for  a 
simplified  fuselage  inlet  in  the  following  sketch.  Plane  (1) 
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contains  traverses  parallel  to  the  duct  outboard  wall,  in 
plane  (2)  they  are  outboard  along  the  extension  of  the  ramp 
leading  edge,  and  in  plane  (5)  the  traverses  run  parallel  to 
the  fuselage  centerline.  Traverses  are  located  in  plane  (3) 
passing  through  the  upper  outboard  corner  of  the  inlet  opening 
and  those  traverses  in  plane  (5)  lie  in  the  plane  of  the  inboard 
duct  wall.  The  length  of  the  traverses  are  made  to  cover  the 
range  of  influence  of  the  inlets.  In  these  flow  field  calcula¬ 
tions,  the  only  effects  included  are  those  of  the  source  panels 
laid  out  on  the  inlet  opening  and  the  panels  on  the  upper, 
lower,  and  outboard  duct  walls  sharing  corner  points  with  the 
panels  on  the  opening.  The  object  is  to  proceed  along  each 
traverse  and  determine  the  first  influence  from  inlet  source 
panels.  The  Prandtl-Glauert  factor,  $I,  associated  with  the 
inlet  panels  and  their  immediate  neighbors  on  the  duct  walls 
is  given  by  Equation  (20)  so  that  the  region  of  influence  of  a 
given  shared  corner  (A,B,C,  or  D)  is  the  same.  In  the  sketch 
below,  the  first  influence  comes  from  corner  A.  It  is  followed 
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in  rapid  succession  by  influences  from  corners  B,  C,  and  D. 

2 

Over  this  short  range,  the  maximum  m  the  quantity  vCmax  1S 
searched  for. 


„2  _  (v2  +  w2) 

vc  ^ 

max  V 

00 


(55) 


Here,  v  and  w  are  the  perturbation  velocities  induced  by  the  set 
of  source  panels  on  the  inlet  and  the  panels  on  the  side  walls 
nearest  the  inlet  face.  The  shock  "surface"  for  the  full  flow 
case  is  then  laid  as  follows.  In  the  region  below  the  fuselage 
and  inlets,  the  shock  is  assumed  to  be  attached  to  the  cowl 
lip  and  its  inboard  extension.  Actually,  it  is  the  continuation 
of  the  oblique  shock  which  is  taken  slightly  ahead  of  the  cowl 
lip.  In  planes  2  and  4  ,  the  shock  starts  from  corner  0 
shown  in  the  first  sketch.  The  shape  of  the  shock  surface  in 
the  vertical  and  other  planes  is  then  computed  from  the  cal¬ 
culated  "shock  slope." 


slope  = 


Az 


'  ^  -  V”Lcal 


-  1 


shock 


(56) 


Referring  to  the  preceeding  sketch,  the  perturbation  velocities 
associated  with  Vc  ,  Equation  (55) ,  are  used  to  compute  a 
local  Mach  number  for  a  given  traverse.  The  approximate  method 
used  to  accomplish  this  is  based  on  Prandtl-Meyer  compression 
and  is  described  in  Section  4.1.1.  Knowing  the  local  Prandtl- 
lauert  factor,  for  several  traverses  allows  the  calculation 

of  the  shock  shapes  in  a  given  plane  in  accordance  with  the 
following  integration 


x 


s 


rz2 

6  o  dz 

t, average  s 

Z1 


(57) 
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For  the  first  step,  g^  calculated  for  the  first  traverse  is 

used  in  Equation  (57) .  Ihe  shock  shape  integration  is  continued 

up  to  a  specified  number  of  traverses  in  the  planes  (I)  through 

(5)  and  continued  until  6^  equals  approximately  g^.  In  Equation 

(57)  ,  gp  is  obtained  from 

l, average 


p£, average  2 

It  is  clear  that  in  essence  this  shock  determination  method  is 
analogous  to  the  method  for  computing  body  nose  shocks  described 
in  the  preceeding  sections  of  this  report.  Tables  containing 
(x  ,y  , z  )  coordinates  are  calculated  for  all  5  traverse  planes 
shown  in  an  earlier  sketch.  Interpolation  schemes  are  set  up 
to  specify  the  shock  surfaces  in  between  the  traverse  planes. 
Above  the  traverse  plane  the  inlet  effects  are  assumed  to 
propagate  along  free-stream  Mach  cones. 

Now  consider  the  partial  blockage  case  with  the  inlet 
operating  on  design  over  the  upper  portion  of  the  inlet 
opening.  We  will  first  treat  the  shock  shapes  in  the  vertical 
planes  (l),  (3),  and  (5).  Refer  to  the  following  sketch  and  the 
first  sketch  in  this  section.  In  the  inlet  opening,  the  lowest 
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panel  is  blocked  and  has  associated  with  it  the  flow  tangency 
condition,  Equation  (23) .  The  following  simplifying  assumption 
is  made.  The  normal  shock  protrudes  down  from  the  most  forward 
leading  edge  of  the  blocked  panels.  From  its  intersection  with 
the  extension  of  the  duct  wall,  the  coalesced  shock  surfaces  are 
determined  in  the  manner  described  above  for  the  unblocked  or 
full  flow  case.  In  planes  (5)  and  (T),  the  shock  shapes  are 
obtained  from  the  integration  starting  from  the  upper  outboard 
corner  O. 


4 . 3  Wing  Shock  Wave  System 

The  nonlinear  correction  used  to  position  the  wing  shock 
wave  system  is  a  modification  and  extension  to  that  described 
in  Section  3.4  of  Reference  1.  In  that  reference  only  the  wing 
leading-edge  shock  wave  was  considered.  In  the  present  work  a 
correction  is  also  made  to  position  the  expansion  associated 
with  the  wing  trai ling-edge  in  the  correct  location.  The 
assumption  is  made,  as  in  Reference  1,  that  the  shock  patterns 
from  the  leading  and  trailing  edges  depend  on  wing  thickness 
only  and  rotate  with  the  wing  at  angle  of  attack.  This  assump¬ 
tion  is  based  on  experimental  evidence  as  discussed  in  Section 
5.2  of  Reference  1, 

The  wing  coordinate  system  is  shown  in  Figure  5.  For  a 
given  field  point  z^)  under  the  wing  at  which  the  flow 

field  is  to  be  calculated,  an  axial  traverse  is  laid  out  parallel 
to  the  x^j  axis  at  coordinates  yw,z^.  The  nonlinear  correction 
to  linear  theory  to  be  described  determines  the  Mach  number 
behind  the  leading- edge  shock  wave  and  determines  the  shock 
location  along  the  traverse  as  well  as  the  Mach  number  just 
ahead  of  the  trailing  edge  shock  and  the  shock  location  along 
the  traverse.  Let  us  first  consider  the  leading  edge  correction. 

The  perturbation  velocities  due  to  wing  thickness  are 
calculated  using  the  linear  theory  method  described  in  Section 
3.3.4.  For  example,  the  upwash  velocity  due  to  one  source 
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panel  is  given  by  Equation  (32) .  The  velocities  are  calculated 
along  the  traverse  beginning  at  the  root-chord  leading  edge 
and  extending  to  the  point  where  the  wing  effect  is  first  felt 
as  shown  in  the  following  sketch.  Once  this  point  is  found  a 


oo 


careful  search  is  made  downstream  of  this  point  for  the 
maximum  flow  angle  Av  where 


tan  Av 


(59) 


and  the  total  axial  velocity,  U,  is  given  by 


JJ_  =  , 

Vcc  V„ 


(60) 


The  velocities  u^,  v^,  and  ww  are  the  wing  thickness  induced 
velocities  positive  in  the  xw,  yw,  and  zw  directions,  respec¬ 
tively,  shown  in  Figure  5.  With  Av  determined.  Equations  (44), 
(45) ,  and  (46)  are  used  to  calculate  M^.  The  component  of 
in  the  axial  direction,  which  is  associated  with  the  wing 
leading  edge,  is  then 


M1cos  Av 


(61) 
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This  is  the  local  Mach  number  which  will  be  associated  with 
the  flow  immediately  behind  the  leading  edge  shock  wave  at  this 
yw,zw  location. 

Under  the  assumption  that  the  shape  of  the  shock  associated 
with  the  wing  leading  edge  depends  on  wing  thickness  only,  the 
scheme  to  be  described  next  allows  for  the  determination  of  the 
shock  location.  The  assumption  is  justified  by  evidence 
offered  by  comparison  with  experimental  data  in  Section  5.2 
of  Reference  1. 

The  summed  effects  of  the  source  panels  on  the  wing  are 
computed  again  at  points  along  the  same  axial  traverse  dis¬ 
cussed  previously.  However,  the  Mach  number  used  in  this 
calculation  is  M^e  calculated  above  instead  of  the  free-stream 
value.  The  result  is  a  velocity  component  profile  similar 
in  shape  to  the  profile  obtained  on  the  basis  of  the  free-stream 
Mach  number  but  shifted  in  the  upstream  direction.  In  other 
words,  the  point  on  the  xw  traverse  where  the  first  wing  thick¬ 
ness  effects  are  felt,  x^e,  now  lies  ahead  of  the  point  shown  in 
the  previous  sketch. 

Different  local  Mach  numbers,  M^e,  can  be  calculated  for 
additional  traverses  parallel  to  the  xw  traverse  and  located  in 
the  came  yw  plane  using  the  method  just  described.  After  re¬ 
computing  the  flow  field  due  to  wing  thickness  along  the 
traverses  with  the  appropriate  local  Mach  numbers,  M^e,  the 
points  where  the  first  influence  is  felt  describe  a  locus  which 
is  taken  as  the  shock  shape.  A  typical  result  is  indicated 
in  the  sketch  below.  Points  1,  2,  and  3  are  the  locations  where 
wing  thickness  effects  are  first  calculated  on  the  basis  of 
local  Mach  numbers,  M^ei,  M^^,  and  M ie^,  respectively.  The 
asymptote  of  the  locus  determined  this  way  runs  parallel  to 
the  Mach  line  based  on  the  free-stream  Mach  number. 
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The  Mach  number  just  ahead  of  the  wing  trailing  edge 
shock  wave  and  the  shock  location  at  yw,zw  are  determined  in 
a  manner  very  similar  to  that  used  for  the  leading  edge.  The 
thickness  induced  upwash  velocity  distribution  under  a  wing 
at  constant  yw,zw  has  the  characteristic  shape  shown  in  the 
following  sketch.  Just  behind  the  leading  edge  shock  a  down- 


upwash 


*w 


wash  (negative  upwash)  exists  which  reduces  with  downstream 
distance.  As  the  trailing  edge  shock  wave  is  approached, 
upwash  is  produced  which  abruptly  returns  to  zero  after  passing 
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through  the  shock.  To  determine  the  Mach  number  just  ahead 
of  the  shock  the  linear  theory  wing  thickness  velocities  are 
calculated  along  the  xr7  traverse  and  a  search  is  made  for  the 
maximum  flow  angle,  Av,  prior  to  the  sudden  drop  in  upwash 
velocity.  Equation  (59)  and  (60)  are  used  to  calculate  Av. 

This  change  in  flow  angle  is  treated  as  an  expansion  of  the 
free-stream  flow  and  the  Prandtl-Meyer  angle  at  this  point  is 

v,  --  v  +  Av  (62) 

where  v^  is  given  by  Equation  (45)  .  Equation  (46)  is  used  to 
determine  M-^.  The  component  of  M1  in  the  axial  direction  just 
ahead  of  the  wing  trailing-edge  shock  wave  is  taken  to  be 

Mfc6  =  M^cos  Av  (63) 

The  location  of  the  trailing-edge  shock  wave  is  found  by 
repeating  the  above  search  using  Mach  number  Mfce  in  the  thickness 
velocity  calculation.  This  time  the  point  at  which  the  upwash 
velocity  drops,  xtQ,  is  searched  for  rather  than  the  maximum 
flow  angle.  Since  the  flow  is  an  expansion  of  the  free-stream 
flow,  Mfce  >  Mm,  xte  is  downstream  of  the  corresponding  point 
calculated  using  straight  linear  theory.  Comparisons  with  data 
to  be  presented  later  show  the  wing  crailing-edge  influence 
location  to  be  welu  predicted. 

The  complete  wing  trailing-edge  shock  wave  shape  at  a 
particular  yw  can  be  determined  by  repeating  the  above  calcula¬ 
tions  for  a  series  of  zw  positions  as  was  done  for  the  leading- 
edge  shock  wave. 

4.4  Store  Imaging 

In  the  work  of  Reference  1  store-aircraft  interference 
was  accounted  for  by  including  store  effects  in  the  wing- 
fuselage-pylon  u-velocity  panel  boundary  condition.  It  was 
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found  that  to  model  accurately  the  reflection  of  the  store  nose 
shock  wave  an  extremely  large  number  of  panels  had  to  be  dis¬ 
tributed  over  the  affected  areas  of  the  aircraft.  In  addition, 
the  u-velocity  panel  strengths  would  have  to  be  recalculated  at 
each  point  in  the  trajectory  since  the  point  at  which  the  store 
shock  wave  strikes  the  wing-fuselage-pylon  combination  would 
move  and  thus  change  the  boundary  condition. 

To  circumvent  these  situations  an  imaging  scheme  has  been 
adopted  in  which  the  store  shock-wave  reflection  is  accounted 
for  by  imaging  the  store  above  the  wing  and/or  inside  the  fuse¬ 
lage.  The  imaging  schemes  will  now  be  described.  The  methods 
of  imaging  the  circular  store  modeled  by  line  singularities 
with  respect  to  the  wing  and  fuselage  will  first  be  discussed. 

This  will  be  followed  by  a  similar  discussion  for  the  elliptic 
store  which  is  modeled  by  source  panels. 

4.4.1  Circular  Store  Imaging 

A  sketch  of  the  circular  store  wing  imaging  method  is 
shown  in  Figure  11.  This  sketch  shows  the  method  for  zero 
degrees  angle  of  attack.  At  angle  of  attack  the  shock  wave 
shape  is  assumed  to  be  the  same  as  at  zero  degrees  angle  of  attack 
with  the  shock  wave  axis  aligned  with  the  free-stream  ve]ocity 
vector.  The  imaging  scheme  places  the  image  store  in  a  mirror 
image  position  relative  to  the  local  wing  chordal  plane.  The 
first  step  is  to  determine  the  location  at  which  the  calculated 
store  shock  wave  intersects  the  surface  of  the  local  airfoil 
section,  that  is,  the  section  immediately  above  the  store.  To 
do  this,  the  vertical  distance,  z  shown  in  Figure  11,  from  the 

S 

store  nose  to  the  local  wing  chordal  plane  is  determined.  Next, 
the  point  at  which  the  store  shock  wave,  calculated  as  described 
in  Section  4.1.1,  intersects  this  plane  is  found.  If  this  point 
lies  ahead  of  the  local  wing  leading  edge  or  behind  the  trailing 
edge,  the  store  is  not  imaged.  With  this  intersection  point, 
xw  ,  determined,  the  wing  thickness  slope  distribution,  tan  0, 

"S 

is  integrated  to  determine  the  airfoil  section  half  thickness 
at  this  point. 
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tan  6  dxw 


(64) 
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It  is  then  assumed  that  the  shock  strikes  the  local  airfoil  t 

surface  at  zg  =  zg  -  t/2.  The  shock  wave  shape  is  used  to 

determine  the  axial  location  of  the  intersection,  x*  ,  and  the 

" s 

slope  of  the  shock  wave,  6S,  at  this  point. 


drshock 
dxshockJ r=z , 

This  angle  is  shown  in  the  following  sketch. 


6„  =  tan 
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(65) 


The  next  step  in  the  calculation  is  to  determine  the 
incidence  angle  between  the  shock  wave  and  the  normal  to  the 
local  airfoil  section,  0i.  This  angle  is  determined  from  0S 
and  the  slope  of  the  local  thickness  envelope,  tan  0.  The 
thickness  slope  angle  is 

0,  =  tan  1  (tan  0)  *  (66) 

*  xw-xwg 
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The  shock  wave  is  assumed  to  reflect  back  at  angle  9^,  which 
is  equal  to  0^.  The  reflection  angle  relative  to  the  local 
xw  axis  is 


er  =  es  +  26t 


The  point  at  which  the  reflected  shock  wave  strikes  the 
axis  of  the  real  store  is  determined  in  the  following  manner. 
In  the  local  airfoil  coordinate  system  the  store  nose  is 
located  at  xWr,  yWn,  zwn  as  shown  in  the  following  sketch 
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The  storo  center  of  gravity  is  also  located  in  this  system 

assuming  zero  yaw  angle  of  the  store,  that  is  YwCg  =  ywn* 

Then  with  z*  =  t/2,  the  equation  of  the  straight  reflected 
Ws 

shock  wave  is 


XW  )  tan  6r 
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The  equation  of  the  line  representing  the  store  axis  is 


where 


tan  a 


tan  a 


(70) 
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If  0  <  a  there  is  no  intersection  and  therefore  no  shock 

r  ~ 

wave  reflection.  For  0  >  a  the  intersection  of  the  reflected 

r 

shock  wave  and  the  store  body  axis  is  determined  from  Equations 
(69)  and  (70). 


xw  = 

f  1 

tan  a  -  tan  0 
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In  the  store  coordinate  system  the  point  at  which  the  reflected 

shock  wave  strikes  the  store  x  axis  is 

s 

xr  =  ~[xw  “  xWn)cos  a  +  (ZW  ”  ZW  )sin  a  (74) 


To  position  the  first  influence  of  the  image  store,  see 
Figure  11,  at  this  location  requires  that  the  first  source 
and  doublets  on  the  image  store  influence  the  flow  field  at 
on  the  real  store.  This  is  accomplished  by  assigning  the  Mach 
cone  angle  0N  to  the  first  singularities.  As  shown  in  Figure 
11,  BN  =  xr/2i*  Thus,  the  local  Mach  number  which  is  associated 
with  the  nose  is  =  vi^+l.  When  the  store  shoulder  is 

reached  it  is  assumed  that  the  free-stream  Mach  number,  M  , 

00 

characterizes  the  local  flow.  The  shoulder  is  taken  to  be 
the  point  along  the  body  where  the  local  surface  slope  first 
becomes  zero. 
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The  strengths  of  the  sources  and  doublets  on  the  real 

store  are  calculated  on  the  basis  of  linear  theory  using  the 

method  described  in  Section  A. 2  of  Appendix  A.  They  are 

calculated  by  starting  at  the  nose  and  marching  downstream. 

The  boundary  condition  includes  not  only  the  free-stream 

component  but  also  all  of  the  parent  aircraft  induced 

velocities  and,  behind  the  reflected  shock  wave,  those  due  to 

the  image  store.  If  the  store  is  pitching  and  yawing, 

velocities  due  to  this  motion  are  also  included.  The  image 

store  singularity  strengths  are  the  same  as  those  on  the  real 

store  however  the  value  of  $  used  in  the  velocity  calculations 

changes  depending  on  location  on  the  real  store.  Referring  to 

Figure  11,  0N  is  used  for  the  point  at  which  the  reflected 

shock  wave  strikes  the  real  store  axis  and  3  is  used  at  and 

00 

aft  of  the  point  where  the  0^  Mach  cone  originating  at  the 
image  store  shoulder  strikes  the  axis.  Between  these  two 
points,  the  local  Mach  number,  M^,  is  assumed  to  vary 
linearly  with  distance  from  Mjj  =  to  Mro  and  the  3  used 

is  0£  =  Vm|  -  1.  This  provides  a  smooth  transition  from  0N 
to  3^,  as  the  calculation  proceeds  down  the  body. 

A  somewhat  similar  procedure  is  used  for  imaging  the 
store  inside  the  fuselage.  A  sketch  of  the  imaging  method  is 
shown  in  Figure  12,  The  shock-wave  reflection  is  calculated 
in  the  same  way  as  was  done  for  the  wing  imaging  using  the 
fuselage  surface.  With  xr  determined,  two  stores  are  placed 
inside  the  fuselage  using  the  method  of  Reference  15.  The 
image  store  is  placed  at  r  /zB  and  a  center  store  is  placed 
on  the  fuselage  axis.  The  image  store  source  strengths  are 
set  equal  to  those  of  the  real  store  and  the  doublet  strengths 
are  reduced  in  accordance  with  Reference  15.  That  is 
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The  center  store  only  has  sources  whose  strengths  are  the 
negative  of  those  on  the  real  store. 
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Since  the  reflected  shock  wave  strikes  the  real  store  at 
xr,  the  first  singularities  of  both  the  image  store  and  the 
center  store  must  influence  the  flow  field  at  this  point. 
Thus,  the  Mach  cone  angles  for  the  first  singularities  are 
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The  Mach  cone  angle  variation  used  on  the  image  store  is  cal¬ 
culated  in  the  same  manner  as  described  previously  for  wing 
imaging.  The  Mach  cone  angle  variation  used  for  the  sources 
on  the  center  store  is  calculated  such  that  corresponding 
sources  on  the  image  store  and  center  store  influence  the  real 
store  at  the  same  point. 

4.4.2  Elliptic  Store  Imaging 

The  imaging  of  an  elliptic  store,  which  is  modeled  by 
source  panels,  with  respect  to  either  the  wing  or  fuselage 
is  performed  in  a  manner  similar  to  that  described  in  Section 
4.4.1  for  imaging  a  circular  store  with  respect  to  the  wing. 
For  the  elliptic  store,  fuselage  imaging  is  also  done  with 
respect  to  a  reflection  plane.  The  source  panels  account  for 
both  volume  and  lift  effects  so  that  imaging  using  the  circle 
theorem  approach  shown  in  Figure  12  is  not  applicable. 

The  elliptic  store  shock  wave  shape  is  not  axisymmetric 
and  is  calculated  as  a  function  of  meridian  angle  around  the 
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store  as  described  in  Section  4.1.2.  Since  the  store  can  be 
at  arbitrary  attitude  relative  to  the  wina,  interpolation  is 
performed  in  the  calculated  shape  to  define  the  shape  in  a 
plane  defined  by  the  store  nose  and  the  local  streamwise  air¬ 
foil  section  immediately  above  the  nose.  Once  the  shock  shape 
in  this  plane  is  known  the  shock  wave  reflection  is  calculated 
in  the  same  manner  as  described  for  the  circular  store  wing 
imaging. 

The  strengths  of  the  source  panels  on  the  real  elliptic 
store  are  calculated  as  described  in  Section  3.2.6  which 
references  Section  3.2.2.  The  boundary  condition  is  given  by 
Equation  (17).  The  velocity  ^  in  that  equation  is  the 

velocity  normal  to  the  panel  control  point  and  now  includes 
not  only  a  component  of  the  free  stream  but  also  components 
due  to  the  parent  aircraft,  the  image  store,  and  the  store 
yaw,  pitch,  and  roll  angular  motions.  The  latter  are  damping 
velocit  ies. 

A  typical  source  panel  layout  on  an  elliptic  body  is 
shown  in  Figure  3.  The  panel  strengths  are  solved  for  by 
starting  at  the  ring  at  the  nose  and  marching  downstream  ring 
by  ring.  Velocities  induced  at  a  ring  of  control  points  by 
upstream  rings  of  known  panel  strengths  are  included  in  the 
boundary  condition.  The  determination  of  the  3  to  use  in  the 
image  store  velocity  calculation  is  identical  to  that  described 
in  the  preceeding  section  for  the  circular  store. 

As  was  mentioned  earlier,  the  elliptic  store  modeled  by 
source  panels  is  imaged  inside  the  fuselage  with  respect  to  a 
reflection  plane.  To  locate  this  plane  the  first  step  is  to 
determine  the  shock  wave  shape  to  be  used  in  the  calculation. 
This  is  done  in  a  manner  similar  to  that  used  for  the  wing 
imaging.  A  plane  is  constructed  which  contains  the  store  nose 
and  the  fuselage  longitudinal  axis.  Interpolation  in  the 
tabulated  shock  wave  shape  is  performed  to  determine  the  shock 
wave  shape  in  this  plane.  The  intersection  of  the  shock  with 
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the  fuselage  surface  is  calculated.  This  determines  the  loca- 

( 

tion  of  the  reflection  plane  relative  to  the  store  nose.  This 

reflection  plane  is  treated  like  the  wing  chordal  plane  for 

wing  imaging.  At  the  point  where  the  shock  wave  strikes  the 

fuselage  surface,  the  shock  wave  slope  and  the  fuselage  surface  * 

slope  are  determined.  The  location  at  which  the  reflected 

shock  wave  strikes  the  real  store  body  is  calculated  as  for 

wing  imaging  with  the  fuselage  surface  slope  replacing  the 

wing  thickness  envelope  slope.  The  strengths  of  the  source 

panels  are  calculated  in  the  manner  just  described  for  the 

wing  imaged  elliptic  stwre. 
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5.  PARENT  AIRCRAFT  FLOW-FIELD  CALCULATION 

This  section  of  the  report  will  describe  the  procedure 
used  to  calculate  the  parent  aircraft  induced  flow  field  at 
a  field  point  such  as  a  point  on  a  separated  store  during  a 
trajectory  calculation.  For  the  present  discussion  the  parent 
aircraft  is  defined  to  be  all  aircraft  components  except  the 
store  whose  trajectory  is  to  be  calculated.  The  linear  theory 
potential  flow  methods  used  to  model  the  various  aircraft 
components  were  described  in  Section  3.  Nonlinear  corrections 
to  these  methods  which  are  made  to  position  shock  waves 
properly  were  described  in  Section  4. 

The  velocities  will  be  calculated  in  the  fuselage  coordi¬ 
nate  system  shown  in  Figure  5,  the  inertial  coordinate  system 
used  in  the  trajectory  calculation.  For  each  aircraft  com¬ 
ponent  the  field  point  is  located  in  the  coordinate  system 
associated  with  the  component.  The  velocities  are  calculated 
in  that  coordinate  system  and  then  resolved  into  the  fuselage 
system. 

5.1  Circular  Fuselage 

The  first  step  in  calculating  the  velocities  at  a  field 
point  is  to  calculate  the  velocities  due  to  the  fuselage.  If 
the  fuselage  is  circular  and  modeled  by  line  singularities. 
Equation  (16)  gives  the  axial,  radial,  and  tangential  velocities 
at  point  (xn,r_ ,0).  This  coordinate  system  is  shown  in  the 

D  D 

sketch  at  the  end  of  Section  3.2.1  and  is  related  to  the  fuse¬ 
lage  system  by 
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The  velocities  ufi,  v0,  and  wR  in  Equation  (16)  are  trans¬ 
formed  to  velocities  in  the  fuselage  system  by 


*\ 

UF  =  “UB 

vF  =  vBsin  Q  +  wBcos  6 
wF  =  -vficos  9  +  wfisin  6 


(79) 


The  linear  theory  source  and  doublet  constants,  Kn  and 
K,  ,  are  used  in  Equation  (16)  but,  because  of  the  nonlinear 

ci  f  n 

shock-wave  shape  calculated  as  described  in  Section  4.1.1,  the 
Prandtl-Glauert  factor  0  will  not  always  be  that  associated  with 
the  free-stream  Mach  number.  Consider  the  following  sketch 


which  shows  the  plane  formed  by  the  field  point  and  the  fuselage 
longitudinal  axis.  The  nonlinear  shock  wave  in  this  plane  is 
also  shown  and  intersects  the  axial  traverse  containing  the 
field  point  at  xg.  To  cause  the  first  influence  of  the  fuselage 
to  occur  at  this  point,  the  Mach  cone  angle  for  the  first 
singularities  must  be  =  xg/rB.  If  the  field  point  were  at 
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Downstream  of 


this  point,  8N  would  be  used  in  Equation  (16) . 

x„  the  flow  is  assumed  to  return  to  the  free-stream  Mach 
s 

number,  M  .  The  location  on  the  traverse  at  which  this  occurs 

00 

is  taken  to  be  the  point  at  which  a  free-stream  Mach  cone 
originating  at  the  shoulder  location  strikes  the  traverse,  xgh. 
If  the  field  point  is  downstream  of  xg^  then  8^  =  VM£  ~  i 
is  used  in  Equation  (16). 


For  field  points  between  xq  and  xgh  the  8  to  use  in 
Equation  (16)  is  determined  by  linear  interpolation  in  the 
Mach  numbers  at  these  two  points,  MN  =  V 8^  +  1  and  M^,  for 
where 


M*  =  “n  +  (Moo  ‘  *n> 


XB  "  Xs 


xsh  ~  xs 
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and  _ 

8  =  Vm|  -  1 


(81) 


5.2  Ejector  Rack  and  Circular  Stores 

The  velocities  induced  at  a  field  point  by  an  ejector  rack 
or  circular  stores  modeled  by  line  singularities  present  on  the 
parent  aircraft  are  calculated  in  the  same  manner  as  has  just 
been  described  for  the  circular  fuselage.  The  field  point  is 
located  in  the  coordinate  system  of  a  particular  body  and  the 
8  to  be  used  in  Equation  (16)  is  determined.  After  the  axial, 
radial,  and  tangential  velocities  are  calculated,  they  are 
resolved  into  axial,  lateral,  and  vertical  components  in  the 
body  coordinate  system  and  then  transfomred  to  components  in 
the  fuselage  system. 
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5.3  Noncircular  Fuselage 


The  noncircular  fuselage  induced  velocities  are  calculated 
as  described  in  Section  3.5  in  the  paragraph  following 
Equation  (37) .  The  value  of  0  to  be  used  in  the  calculation 
is  determined  in  the  same  way  as  for  the  circular  fuselage. 

The  source  panel  strengths  calculated  on  the  basis  of  linear 
theory  are  used.  The  nonlinear  shock  wave  shape  in  the  plane 
containing  the  field  point  and  the  fuselage  longitudinal  axis 
is  found  by  interpolation  in  the  shock  wave  shapes  which 
were  calculated  in  the  manner  described  in  Section  4.1.2. 

5.4  Engine  Inlets 

If  the  noncircular  fuselage  has  engine  inlets,  the  effect 
of  the  inlet  shock  wave  on  the  flow  field  must  be  accounted 
for.  The  flow  model  for  the  inlets  was  described  in  Section 
3.2.3  and  the  calculation  of  the  inlet  shock  wave  system 
discussed  in  Section  4.2.  At  field  points  ahead  of  the  inlet 
shock  wave  the  flow  field  produced  by  the  noncircular  fuselage 
is  calculated  as  described  in  the  proceeding  paragraph.  When 
the  field  point  lies  on  or  behind  the  inlet  shock  the  procedure 
to  be  described  is  used. 

At  points  in  tht.  region  affected  by  the  fuselage  mounted 
inlets,  the  flow  field  velocities  are  determined  including 
effects  of  the  inlet  shock.  In  many  ways  the  procedure 
resembles  the  one  employed  for  body-nose  shocks  described 
previously  in  this  section.  Consider  the  case  shown  in  the 
sketch  below.  A  coalesced  shock  emanates  from  the  inlet 
opening  corresponding  to  the  partially  blocked  case.  The 
unblocked  or  full-flow  inlet  is  a  degenerate  form  of  the 
situation  shown.  These  two  different  conditions  can  occur 
within  the  framework  of  the  shock  determination  method 
discussed  in  Section  4.2.  We  will  consider  velocity  traverses 
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in  the  vertical  planes  (X),  (3),  and  (5)  shown  in  the  first 
sketch  in  Section  4.2  but  the  method  applies  to  the  ether 
planes  as  well. 

For  a  given  axial  traverse,  the  point  of  intersection  with 
the  inlet  shock  is  determined  using  the  shock  surface  coordi¬ 
nates  calculated  by  the  method  described  earlier.  These  points, 
designated  and  PI^  sketch,  can  lie  ahead  or  aft 

of  the  oblique  shock  extension.  Point  0  is  the  leading  edge, 
in  sideview,  of  the  blocked  panels  on  the  inlet  opening  nearest 
the  cowl  lip.  For  the  unblocked  or  full  flow  case,  point  0 
is  taken  at  the  ramp  leading  edge.  This  point  would  also  be 
used  for  the  fully  blocked  case.  As  such  this  point  is  the 


4 


origin  of  the  shock  system  associated  with  the  inlet  for  both 

cases.  The  procedure  is  to  determine  a  Prandtl-Glauert  factor 

8_  to  position  the  inlet  shock  influence  in  the  correct 
1 1-  s 

location  and  to  return  the  flow  downstream  to  the  free-stream 
value  8^.  For  points  along  the  traverse  between  these  two 
points  a  single  value  of  will  be  used  for  all  fuselage  and 
inlet  panels  until  8^  becomes  equal  to  the  oblique  shock 
Mach  cone  angle.  Downstream  of  this  point  8j  will  be  used 
for  the  inlet  panels  and  8^  for  the  fuselage  panels. 

First  consider  traverse  z^  shown  in  the  sketch.  Here 
the  traverse  intersects  the  shock  ahead  of  the  oblique  shock 
extension.  This  intersection  point,  PlZl,  gives  the  first  or 
starting  value  of  8^,  8j  g.  It  is  obtained  geometrically  by 
connecting  the  intersection  point  PIZl  to  shock  origin  0.  At 
points  on  the  traverse  behind  the  shock,  the  value  of  the 
factor  8^  is  made  to  linearly  vary  from  the  starting  value  to 
8^.  The  length  over  which  this  variation  occurs  ranges  from 
intersection  point  PIZ^  to  an  end  point  fixed  as  follows.  A 
8^  Mach  line  is  drawn  from  a  point  on  the  duct  upper  surface 
located  axially  at  the  cowl  lip  location,  refer  to  the  pre¬ 
vious  sketch.  This  Mach  line  lies  in  the  same  vertical  plane 
as  the  traverse  under  consideration,  and  the  point  of  their 
intersection  with  axial  coordinate  x3  ^  is  the  point  where  8^ 
equals  8^.  The  effects  of  all  source  panels  distributed  on 
the  fuselage  and  inlet  opening  are  then  calculated  by  this 
method  at  points  along  a  traverse  such  as  z^. 

A  slightly  different  approach  is  required  when  the  point 

of  intersection  PI  lies  aft  of  the  oblique  shock  extension, 

z 

as  is  the  case  for  traverse  Z2  shown  in  the  sketch.  The 
Prandtl-Meyer  factor  is  still  determined  by  the  method  just 
described.  However,  since  the  line  connecting  PI^  to  0  now 
lies  aft  of  the  oblique  shock  extension  (or  the  8j  line) ,  the 
source  panels  on  the  inlet  face  will  not  allow  a  solution  with 
8t  a|_  •  To  circumvent  this  problem,  the  axial  coordinates  of 
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points  on  the  traverse  on  and  behind  the  shock  are  mapped  for¬ 
ward  in  accordance  with 


X  =  v 

s,map  s,» 


X-  -  X 
_ Si!L 


XPI  “  xs,« 
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(82) 


where  Xj  is  the  axial  coordinate  of  the  intersection  point  of 
the  traverse  and  the  extension  of  the  oblique  shock  or  Bj  line. 
This  mapping  is  done  in  order  that  the  inlet  shock  will  be 
present  below  the  point  at  which  the  nonlinear  shock  crosses 
the  Bj  line.  The  velocities  calculated  at  these  mapped  points 
are  assumed  to  be  those  at  the  unmapped  points. 


5.5  Elliptic  Stores 

The  velocities  induced  at  a  field  point  by  a  store 
modeled  by  source  panels,  the  elliptic  store  model,  are  cal¬ 
culated  in  the  same  manner  as  described  in  Section  5.3  for 
the  noncircular  fuselage.  The  field  point  is  located  in  the 
store  coordinate  system  and  then  the  plane  containing  this 
point  and  the  store  longitudinal  axis  is  defined.  The  shock 
wave  shape  in  this  plane  is  determined  by  interpolating  in  the 
shock  wave  shapes  calculated  by  the  method  described  in  Section 
4.1.2.  The  linear  theory  panel  strengths  are  used  in  the 
velocity  calculation.  The  value  of  3  to  use  is  determined  in 
the  same  manner  as  was  described  for  the  circular  fuselage 
in  Section  5.1.  After  the  velocities  are  calculated  in  the 
store  coordinate  system  they  are  transformed  to  the  fuselage 
coordinate  system. 

5.6  Wing,  Pylon,  and  Fuselage  Interference  Shell 

The  method  used  to  calculate  the  velocities  due  to  the 
wing,  pylon,  and  fuselage  interference  shell  will  now  be 
described.  The  wing  leading  edge  and  trailing-edge  shock  wave 
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locations  x^e  and  xfce,  and  the  Mach  numbers,  M^e  and  Mtg, 
at  these  locations  are  determined  by  the  method  described  in 
Section  4.3.  The  velocities  due  to  wing  and  pylon  thickness 
and  the  u-velocity  panels  on  the  wing,  pylon,  and  fuselage 
interference  shell  are  calculated  as  described  in  Section  3. 
The  linear  theory  panel  strengths  are  used  and,  for  the  pylon, 
is  used  for  all  field  points. 


For  the  wing  and  fuselage  interference  shell  the  value  of 


to  be  used  is  determined  in  the  following  manner.  If  the 
field  point  is  downstream  of  xtfi,  3^  is  used.  If  the  field 
point  lies  between  x^e  and  xfce  the  value  of  3^  is  determined 
linear  interpolation  in  the  Mach  numbers  at  these  two  points 


by 


Thus,  at  field  point  x, 


(83) 


(84) 


6.  STORE  FORCE  AND  MOMENT  CALCULATION 


The  calculation  of  the  trajectory  of  a  store  separating 
from  an  aircraft  requires  a  knowledge  of  the  forces  and  moments 
acting  on  the  store  at  each  point  in  the  trajectory.  In  order 
to  calculate  these  forces  and  moments,  the  nonuniform  velocity 
field  in  which  the  store  is  immersed,  as  seen  by  the  store, 
must  be  determined  at  each  point  in  time.  This  field  includes 
the  free-stream  velocity,  the  perturbation  velocities  induced 
by  the  parent  aircraft,  image  store  velocities,  and  the 
angular  velocities  due  to  the  store's  pitch,  yaw,  and  roll 
motions.  The  distribution  along  the  body  is  required  in 
order  to  calculate  the  forces  and  moments  acting  on  the  body 
and,  if  fins  are  present,  the  distribution  over  the  fin  sur¬ 
faces  is  required  for  the  empennaae  force  and  moment  calculation. 

This  section  will  first  describe  the  calculation  of  the 
velocity  field  utilizing  the  flow  models  of  Sections  3  and  4. 
Following  this  the  methods  employed  to  calculate  the  body  and 
empennage  forces  and  moments  from  this  nonuniform  flow  field 
will  be  presented,  first  for  the  circular  store  and  then  for 
the  elliptic  store. 

6.1  Velocity  Field  as  Seen  by  Store 

The  x  ,y  ,z  coordinate  system  to  be  used  in  the  force 
s  s  s 

and  moment  calculations  is  shown  in  Figure  13.  It  is  fixed 

in  the  store  with  the  origin  at  the  store  nose  with  the 

velocities  U  ,  V  ,  and  W  positive  in  the  directions  shown, 
s  s  s 

These  are  the  total  velocities  as  seen  by  a  point  on  the  store 
and  each  is  composed  of  a  free-stream  component,  perturbation 
components  due  to  the  parent  aircraft  and  image  store,  and  a 
damping  component  due  to  the  store's  rotational  motion. 

Before  writing  the  expressions  for  these  velocities  let 
us  examine  Figure  14  which  shows  the  coordinate  systems  used 
in  the  trajectory  calculation.  The  equations  of  motion  will 
be  presented  in  Section  7.  Two  coordinate  systems  are  shown 
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in  Figure  14.  The  £,n,G  system  is  an  inertial  system  fixed 
in  the  aircraft  which  is  in  rectilinear  flight  at  uniform 
velocity.  The  positive  £  axis  is  forward  along  the  fuselage 
longitudinal  axis,  the  n  axis  is  positive  laterally  to  the 
right,  and  the  £  axis  is  positive  downward  perpendicular  to  the 
fuselage  longitudinal  axis.  The  origin  of  this  sytem  is  fixed 
in  the  fuselage  nose. 

The  x,y,z  coordinate  system  is  fixed  in  the  store  with 
the  origin  at  the  store  moment  center.  Referring  back  to 
Figure  13,  the  x-axis  lies  along  the  store  longitudinal  axis 
but  is  positive  forward.  The  y-axis  is  parallel  to  the  y 

s 

axis  and  positive  in  the  same  direction  and  the  z  axis  is 
parallel  to  the  z  axis  but  positive  downward.  The  x,y,z 
system  shown  in  Figure  13  is  a  right-hand  system  and  the  store 
angular  velocity  components  about  the  three  axes  are  p,  q, 
and  r,  respectively. 

In  Figure  14  the  origins  of  the  two  coordinate  systems 
have  been  drawn  to  coincide  in  order  to  show  the  angles  used 
to  determine  the  orientation  of  the  store  with  respect  to 
the  inertial  (£,n,0  axes.  The  system  of  angles  chosen 
consists  of  three  rotations  in  the  yaw,  ¥;  pitch,  0;  and 
roll,  4;  sequence  (Reference  16).  The  two  coordinate  systems 
are  then  related  by 


where 
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The  three  velocity  components,  Ug,  Vg,  and  Wg  are 
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The  first  term  on  the  right  side  of  the  equal  sign  in  each 
of  the  above  expressions  is  the  component  of  the  free-stream 
velocity  vector.  If  the  parent  aircraft  is  moving  at  velocity, 
Vm,  and  flying  at  angle  of  attacx,  ap,  then  the  velocity  of 
the  inertial  coordinate  system  relative  to  a  point  in  space  is 


^co  =  VcoCOS  +  Veosin 


(88) 


The  velocity  of  the  store  moment  center  relative  to  the 
moving  inertial  system  is 

=  ee5  +  +  (89 

so  that  the  velocity  of  the  store  relative  to  a  point  in  space 
is 


=  ^  +  Vm  =  (V^cos  ap  +  +  nen  +  (V^sin  oip  + 

(90) 

Since  the  velocity  components  needed  in  Equation  (87)  are  those 
seen  by  a  point  fixed  in  the  store  and  moving  with  the  store, 
Equation  (90)  becomes 

i  =  -  (V^cos  ap  +  |)e?  -  nen  -  (V^sin  ap  +  £)e;  (91) 
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The  components  in  the  x,y,z  coordinate  system  of  Figure  13 
are  then 
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where  [A] ’  is  the  transpose,  or  inverse  since  [A]  is  orthogonal, 
of  the  matrix  given  in  Equation  (86) .  These  velocity  components 
are  related  to  those  in  the  coordinate  system  of  Figure  13  by 
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The  second  term  on  the  right-hand  side  in  the  three 
expressions  of  Equation  (87)  is  the  parent  aircraft  perturba¬ 
tion  velocity  term.  This  is  comprised  of  the  perturbation 
velocities  induced  by  the  fuselage,  wing,  pylon,  rack  and 
other  stores  and  is  calculated  by  the  methods  presented  in 
Section  5.  These  velocities  are  summed  up  in  the  fuselage, 
or  inertial,  coordinate  system.  The  velocity  components  in  the 
x,  y,  and  z  directions  of  Figure  13  are 
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and,  finally,  the  components  in  the  x  ,y  ,z  coordinate  system 

8  p  5 

of  Figure  13  are 


u  »  -u,  v  =  v,  w  =  -w  (95) 

P  P  P 

The  third  terms  on  the  right  side  of  Equation  (87)  are 
velocities  induced  at  a  point  on  the  real  store  by  the  image 
store (s) .  The  calculation  of  the  image  store  induced  velocities 
was  described  in  Section  4.4.  Imaging  is  done  in  a  different 
manner  depending  upon  which  store  modeling  option,  circular 
or  elliptic,  is  used.  The  velocities  induced  by  the  image 
store  at  control  points  on  the  real  store  body  and  fins  are 
calculated. 

The  fourth  terms  in  the  last  two  expressions  in  Equation 
(87),  vs^  and  wg^,  are  included  only  when  the  elliptic  store 
option  is  used  and,  then,  only  when  two  sets  of  fins  are 
present  and  the  option  to  include  trailing-edge  vortex  effects 
is  selected.  These  trailing-edge  vortices  are  shed  by  the 
forward  set  of  fins.  Their  paths  are  calculated  back  to  the 
leading  edge  of  the  rear  fins.  With  the  vortices  fixed  at 
these  positions  the  velocities  induced  at  the  fin  control 
points  are  calculated.  The  trailing-edge  vortex  interference 
will  be  described  in  Section  6.3.3. 

The  remaining  terms  in  Equation  (87)  are  aerodynamic 

damping  terms  due  to  the  rotacional  motion  of  the  store.  The 

rotational  velocities,  p,  q,  and  r,  are  about  the  x  ,  y  , 

s  s 

and  z  axes,  respectively,  and  are  positive  in  the  positive 

s 

sense  of  the  moments  shown  in  Figure  13.  The  values  of  p,  q, 
and  r  at  each  point  in  a  trajectory  are  determined  by  integra¬ 
tion  of  the  equations  of  motion. 

The  store  free-stream  velocity  will  change  during  its 
trajectory  because  of  its  motion  relative  to  the  parent  air¬ 
craft.  The  velocities  to  be  used  in  the  force  and  moment 
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calculation  are  made  dimensionless  by  the  store  free-stream 
velocity.  Thus,  Equation  (87)  becomes 
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where  from  Equation  (90)  or  (91) 


v«  -  [(V^cos  ap  +  £)2  +  n2  +  (Vrosin  ap  +  C)2] 
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(97) 


6.2  Circular  Store  Forces  and  Moments 

This  section  of  the  report  will  describe  the  methods  used 
to  calculate  the  aerodynamic  forces  and  moments  acting  on  the 
separated  store  when  the  circular  store  option  is  used.  In 
this  option  the  store  body  is  modeled  by  line  singularities 
and  one  set  of  fins  can  exist  which  can  be  planar  or  cruciform. 

6.2.1  Body  Forces  and  Moments 

In  the  work  of  Reference  1  s lender-body  theory  was  used 
to  calculate  the  forces  and  moments  acting  on  the  store  body. 
This  method  was  found  to  be  quite  inaccurate  for  a  store 
immersed  in  a  nonuniform  flow  field  at  supersonic  speed.  In  the 
present  work  the  slender  body  method  has  been  replaced  by  one 
in  which  the  store  body  is  modeled  by  line  singularities  dis¬ 
tributed  along  the  body  axis.  Supersonic  line  sources  and  sinks 
are  distributed  along  the  store  centerline  to  account  for 
volume  effects  subject  to  a  nonuniform  axial  velocity  distribu¬ 
tion.  The  nonuniform  upwash  and  sidewash  fields  in  which  the 
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store  is  immersed  are  accounted  for  by  two  distributions  of 
supersonic  line  doublets  on  the  store  longitudinal  axis.  One 
doublet  distribution  cancels  the  effects  of  the  upwash  distri¬ 
bution  and  the  other,  the  sidewash  distribution. 


The  derivation  of  the  equations  used  to  calculate  the 
singularity  strengths,  the  surface  pressure  distribution, 
the  axial  load  distributions,  and  the  overall  forces  and 
moments  is  presented  in  Appendix  A.  The  nonuniform  velocity 
field  in  which  the  store  is  immersed  is  calculated  as  was 
just  described  in  Section  6.1.  As  discussed  in  Appendix  A, 
the  velocities  are  calculated  at  points  on  the  store  axis. 
Equations  (A-13) .  (A-22) ,  and  (A-30)  are  used  to  solve  fox 
the  strengths  of  the  three  singularity  distributions.  The 
velocity  distributions  appearing  in  these  three  equations  are 
related  to  the  velocities  of  Equation  (96)  by 


)  (98) 


The  singularity  strengths  modeling  the  store  in  the 
nonuniform  flow  field  have  now  been  determined.  The  next 
step  in  the  store  force  and  moment  calculation  is  to  calculate 
the  store  surface  pressure  distribution.  Sections  A- 3.1  and 
A-3.2  of  Appendix  A  describe  the  surface  velocity  calculation 
due  to  the  singularities  modeling  the  store.  The  surface 
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pressure  coefficient  calculation  is  discussed  in  Section 
A- 3. 3.  The  Bernoulli  equation  for  compressible,  isentropic 
flow  is  used  and  is  given  by  Equation  (A-35) . 

The  question  arises  as  to  what  velocities  should  be 
included  in  the  surface  pressure  distribution  calculation. 
Those  due  to  the  singularities  modeling  the  store  must  be 
included  as  must  the  damping  velocities.  The  velocities  in 
question  are  those  due  to  the  parent  aircraft  and  the  image 
store (s).  These  latter  velocities  will  be  excluded  based  on 
comparisons  with  experimental  data. 

Figures  15  through  18  show  longitudinal  pressure  dis¬ 
tributions  on  an  ogive-cylinder  store  under  the  left  wing  of 
a  wing-fuselage  combination  in  a  position  at  which  the  store 
nose  shock  wave  reflects  off  the  wing  and  strikes  the  store. 
The  data  are  from  Reference  6.  Two  theoretical  curves  are 
shown  in  the  figures.  The  solid  curve  has  the  velocities 
due  to  the  store  singularities  as  well  as  those  due  to  the 
wing,  fuselage,  and  image  store  included  in  the  surface 
velocity  calculation.  These  velocities  determine  the  pres¬ 
sures.  The  dashed  curve  omits  the  image  store  velocities 
which  start  where  these  curves  begin.  As  can  be  seen,  the 
agreement  between  experiment  and  theory  is  quite  good  when 
the  image  store  velocities  are  excluded. 

Reference  17  gives  some  insight  into  why  the  calculated 
pressures  with  the  image  store  included  are  too  large. 
According  to  the  conclusions  drawn  in  this  reference,  the 
following  phenomena  can  occur;  (1)  when  an  oblique  shock 
impinges  on  an  elongated,  cylindrical  body,  the  main  feature 
is  a  strong  boundary-layer  crossflow  and  separation  on  the 
lee  side  accompanied  by  a  rapid  relieving  of  pressures  on 
the  incident  side,  (2)  details  of  the  pressure  distributions 
depend  on  the  nature  and  thickness  of  the  boundary  layer  and 
the  shock  strength,  and  (3)  on  the  incident  side,  the 
presence  of  the  boundary  layer  also  allows  for  upstream 
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propagation  of  the  pressure  rise  associated  with  the  impinging 
shock.  Thus,  the  pressure  distribution  on  a  body  in  the  region 
enveloped  by  an  oblique  shock  is  very  much  subject  to  viscous 
effects.  Prediction  methods  based  on  nonviscous  theories 
tend  to  overpredict  the  magnitude  of  the  peak  pressures,  and 
therefore,  such  methods  require  suitable  corrections.  One 
improvement  may  be  to  neglect  the  presence  of  the  shock  wave 
in  the  pressure  distribution  calculation  in  the  region  enveloped 
by  the  shock  wave. 

This  has  been  done  in  the  calculations  shown  by  the  dashed 
curves  in  Figures  15  through  18.  In  these  calculations  the 
image  store  contribution  to  the  velocity  field  has  been  elimi¬ 
nated  everywhere.  As  can  be  seen  the  agreement  between  the 
calculation  and  the  experimental  data  is  now  quite  good.  A 
similar  improvement  has  been  observed  for  other  flight  condi¬ 
tions  and  store  positions  relative  to  the  parent  aircraft.  It 
has  also  been  found  that  a  similar  situation  occurs  when  the 
wing  shock  wave  strikes  the  store  in  the  cylindrical  region. 

In  the  example  shown  in  the  figures,  the  wing  shock  strikes 
the  store  near  the  tip  of  the  nose  and  the  effect  is  not  as 
apparent.  In  calculating  the  surface  velocities  to  be  used  in 
the  pressure  distribution  calculation  only  the  velocities  due  to 
the  singularities  modeling  the  separated  store  and  rotational 
velocities  are  included  in  the  present  method. 

The  axial  distributions  of  forces  and  moments  are  cal¬ 
culated  from  the  pressure  distribution  using  Equations  (A- 39) 
and  (A-40)  of  Appendix  A  and  the  overall  forces  and  moments  are 
given  by  Equation  (A-41) . 

6.2.2  Empennage  Forces  and  Moments 

The  circular  store  option  allows  one  set  of  planar  or 
cruciform  fins  to  be  placed  on  the  store  body.  The  method  used 
to  calculate  the  empennage  forces  and  moments  is  that  presented 
in  Section  5.3  and  Appendix  I  of  Reference  18.  The  method 
used  there  is  based  on  a  combination  of  reverse  flow  theorems 


and  slender-body  theory  scaled  for  the  effects  of  aspect  ratio 
as  done  in  wing-body  interference  theory,  Reference  19.  The 
effects  of  interference  between  the  body  and  the  fins  are 
accounted  for. 

The  method  requires  that  the  velocity  field  across  the 
span  of  each  fin  be.  known.  fThe  velocities  are  calculated 
as  described  in  Section  6.1.  With  this  velocity  distribution 
known,  the  method  of  Reference  18  yields  the  normal  and  side 
forces  and  the  pitching,  yawing,  and  rolling  moments  due  to 
the  empennage.  The  total  aerodynamic  forces  and  moments 
acting  on  the  circular  store  are  obtained  by  adding  the  store 
body  forces  and  moments  to  those  for  the  empennage. 

6.3  Elliptic  Store  Forces  and  Moments 

This  section  will  describe  the  methods  used  to  calculate 
the  aerodynamic  forces  and  moments  acting  on  the  separated 
store  when  the  elliptic  store  option  is  used.  The  methods 
used  are  modifications  to  those  described  in  Reference  7. 

The  store  body  is  modeled  by  source  panels.  Up  to  two  sets  of 
fins  can  be  positioned  on  the  body  and  each  set  can  hove  from 
one  to  four  fins.  The  fins  and  the  interferes  u  shell  placed 
on  the  body  to  account  for  fin-body  interference  are  modeled 
by  constant  u-velocity  panels  like  those  used  for  the  wing  and 
fuselage  interference  shell  described  in  Section  3.3.  IT  two 
sets  of  fins  are  present,  interference  of  the  forward  set  on 
the  aft  set  can  be  accounted  for.  Positions  and  strengths  of 
vortices  are  calculated  at  ths  trailing  edge  of  the  forward 
set.  Their  positions  at  the  leading  edge  of  the  aft  set  are 
determined  by  integration  of  their  paths. 

6.3.1  Body  Forces  and  Moments 

The  first  step  in  determining  the  forces  and  moments  act¬ 
ing  on  the  elliptic  store  is  to  determine  the  strengths  of  the 
source  panels  modeling  the  body.  The  method  used  is  the  same 
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as  that  described  m  Section  3.3.2  for  the  noncircular  fuselage 
and  in  Section  3.2.6  for  the  elliptic  store.  In  those  two 
sections  the  boundary  condition  only  included  the  free- stream 
flow.  In  calculating  the  panel  strengths  for  use  in  the  force 
and  moment  calculation,  the  boundary  condition  reflects  the 
velocities,  calculated  as  described  in  Section  6.1.  Thus, 
Equation  (17)  is  replaced  by 
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where,  using  the  velocities  of  Equation  (96) , 
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The  angles  6^  and  ^  are  the  panel  orientation  angles  measured 

in  the  store  body  x  ,y  ,z  coordinate  system.  They  are 
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measured  as  shown  in  the  following  sketch.  The  positive  senses 
are  shown. 
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The  panel  strengths  are  determine  by  using  a  marching 
procedure  over  the  length  of  the  store  body.  This  procedure 
solves  for  the  strengths  ring  by  ring  beginning  at  the  store 
nose.  This  is  done  in  order  that  the  image  store  induced 
velocities  us^,vs^,ws.  in  Equation  (87)  can  be  calculated. 
These  velocities  are  determined  in  the  manner  described  in 
Section  4.4.2. 

The  next  step  in  the  force  and  moment  calculation  is  to 
calculate  the  pressure  coefficient  at  the  control  points  of 
all  of  the  panels.  This  calculation  is  identical  to  that 
describe  in  Section  6.2.1  for  the  circular  store  and  presented 
in  Section  A-3.3  of  Appendix  A.  The  Bernoulli  relationship 
given  by  Equation  (A— 35 )  is  used  with  the  velocity  given  by 
Equation  (A-37) .  As  for  the  circular  store,  only  the  velo¬ 
cities  due  to  the  source  panels  and  store  rotational  motion 
are  included  in  the  pressure  calculation  for  the  same  reasons 
discussed  in  Section  6.2.1.  In  Equation  (A-37) 
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where  u  ,  v  ,  and  w  are  the  source  panel  induced  velocities 
s  s  s 

at  a  control  point  in  the  store  body  coordinate  system  shown 
in  Figure  13.  The  remaining  terms  are  the  velocities  at  the 
control  points  due  to  the  store  rotational  motion,  see  Equation 
(87)  . 
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The  forces  and  moments  acting  on  a  panel  are  calculated 
by  first  determining  the  total  force  acting  at  the  control 
point  in  a  direction  normal  to  the  panel  surface,  Cp^.  This 
force  is 


CF.  =  -CP.VSR 


(102) 
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where  Cp.  is  the  pressure  coefficient  at  the  i  panel  control 
r 

point  and  is  the  area  of  the  i  panel.  This  force  is 
resolved  into  normal-force  and  side-force  components  using  the 
angles  0  and  6  shown  in  the  preceeding  sketch  by 


C„  =  C„  cos  0  cos  6 
N. 
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The  axial- force  component  is  not  used  in  the  present  method 
as  an  axial  force  but  is  used  as  a  contribution  to  the 
pitching  and  yawing  moments.  This  component  is 
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The  moment  contributions  due  to  a  panel,  positive  in  the 
directions  shown  in  Figure  13,  are 
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The  total  forces  and  moments  acting  on  the  store  body  are 
determined  by  summing  Equations  (103)  and  (105)  over  all 
of  the  source  panels  modeling  the  store. 

6.3.2  Empennage  Forces  and  Moments 

The  elliptic  store  option  allows  one  or  two  sets  of  fins 
to  be  present  on  the  store  body.  In  addition,  if  there  are 
two  sets  of  fins  interference  of  the  forward  set  on  the  rear 
set  can  be  accounted  for  by  including  the  velocities  induced 
by  the  trailing  edge  vorticity  associated  with  the  forward 
set  in  the  nonuniform  flow  field.  The  methods  used  model 
the  fins  and  fin-body  interference,  determine  the  trailing- 
edge  vortex  strengths  and  positions,  and  calculate  the  vortex 
paths  along  the  body.  They  are  modifications  of  those  pre¬ 
sented  in  Reference  7.  Some  features  included  in  the  work  of 
that  reference  have  been  deleted  and  others  added.  This 
section  will  describe  the  fin-body  model  and  the  calculation 
of  the  forces  and  moments.  The  trailing-edge  vortex  method 
will  be  discussed  in  the  next  section. 

The  methods  used  to  model  the  fins  and  the  fin-body 
interference  were  extracted  from  program  DEMON 2  described  in 
Reference  7.  Certain  features  were  eliminated  such  as  nose 
vorticity  in  the  boundary  condition,  leading-edge  and  side-edge 
suction  calculations,  and  the  option  of  modeling  a  circular 
store  by  line  singularities.  In  the  present  program  the  store 
can  only  be  modeled  by  line  singularities  by  using  the  circu¬ 
lar  store  option.  Fin  thickness  effects  are  also  not  included. 

Other  features  have  been  added  to  the  present  program. 

Each  set  of  fins  can  be  comprised  of  1,  2,  3,  or  4  fins.  They 
can  be  located  circumferentially  anywhere  on  the  body  and  may 
be  at  any  angle  relative  to  the  tangent  to  the  body  surface 
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in  the  crossflow  plane.  The  program  was  modified  to  include 
parent  aircraft  and  image  store  velocities  in  the  boundary 
condition.  Velocities  at  the  control  points  due  to  store 
angular  motion  are  also  included. 

The  fins  are  modeled  in  the  same  manner  as  described 
in  Section  3.3  for  the  wing  and  fuselage  interference  shell. 

The  fins  and  the  interference  shell  placed  on  the  store  body 
for  fin-body  interference  are  represented  by  constant  u- 
velocity  panels.  A  typical  panel  layout  is  shown  in  Figure 
19.  The  interference  shell  shown  in  the  figure  extends  over 
the  length  of  the  root  chord  of  the  fins  and  the  panel  layout 
is  the  same  as  that  used  for  the  source  panels  modeling  the 
store  body.  The  interference  shell  should  extend  beyond  the 
fin  trailing  edge  to  include  all  of  the  lift  carryover  if 
the  body  extends  beyond  the  trailing  edge.  Use  of  the  u- 
velocity  panel  method  has  shown  that  if  this  is  done 
difficulties  are  encountered  in  the  solution  for  the  panel 
strengths.  To  avoid  this  problem,  the  shell  is  terminated  at 
the  trailing  edge.  Most  of  the  lift  carryover  is  still 
accounted  for. 

The  flow  tangency  condition  at  a  fin  control  point  is 

NWBP  VN  V  Vs 

£  _ —  cos  <j)f  +  - —  sin  <}>f  -  sin  6f  (106) 

n=l  »  »  00  „ 

s  s  s 

v  =  1,2, . . .NPANLS 

In  the  above  notation  there  are  NPANLS  on  the  fins  and 
NWBP  panels  on  the  fins  and  the  interference  shell.  The 
summation  on  the  left  side  of  Equation  (106)  is  the  perturba¬ 
tion  velocity  normal  to  the  vth  panel  induced  by  all  of  the 
u-velocity  panels  on  the  fins  and  interference  shell.  The 
velocities  VSv  and  WSv  are  the  sums  of  the  velocities  given  by 
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Equation  (87)  and  those  induced  at  the  v  control  point  by 
the  source  panels  modeling  the  store  body.  The  angle  is 
measured  as  shown  in  Figure  20  which  is  taken  from  Reference 
3.  The  last  term  on  the  right-hand  side  of  Equation  (106) 
is  the  sine  of  the  fin  deflection  angle.  Positive  fin 
deflection  is  such  that  the  fin  produces  a  force  in  the  direc¬ 
tion  of  Cp  shown  in  Figure  20. 

The  flow  tangency  condition  at  an  interference  sheli 
control  point  is 
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V  =  NPANLS+1 ,NPANLS+2 , . . .NWBP 

These  panels  only  account  for  mutual  interference  between 
the  fins  and  the  body.  Hence,  the  right-hand  side  is  zero. 
Solution  of  the  set  of  equations  given  by  Equations  (106) 
and  (107)  following  the  procedure  described  in  Section  3.3 
yields  the  strengths  of  the  constant  u-velocity  panels  model¬ 
ing  the  fins  and  the  associated  interference  shell. 

The  next  step  in  the  empennage  force  and  moment  calcula¬ 
tion  is  to  calculate  the  pressure  coefficient  at  the  control 
points  of  all  of  the  u-velocity  panels  on  the  fins  and  the 
interference  shell.  The  Bernoulli  relationship  given  by 
Equation  (A-35)  of  Appendix  A  is  used  to  calculate  the 
pressures,  as  before,  with  the  velocity  used  in  the  relation¬ 
ship  given  by  Equation  (A-37).  The  velocity  components  in  this 
equation  are  given  by  Equation  (101).  The  velocities  u  ,  v  , 
and  w  in  Equation  (101)  are  the  sums  of  those  induced  at  the 
control  points  by  the  constant  u-velocity  panels,  the  body 
source  panels,  and  external  vortices  if  they  are  present. 

The  total  force  acting  on  a  panel  in  a  direction  normal 
to  the  panel  surface  is  given  by  Equation  (102) . 
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The  positive  sense  of  this  force  on  fins  in  the  four  quadrants 
is  shown  in  Figure  20.  The  positive  sense  on  a  body  inter¬ 
ference  panel  is  outward  from  the  body.  This  force  is  resolved 
into  normal  and  side  force  components  using  the  angle  <f>£  for 
the  fins  and  8_  for  the  body  interference  panels.  It 

f  H  -A.  ir 

should  be  noted  that  <J>f  for  the  fir.  arrangement  shown  in 
Figure  20  is  positive  in  the  upper  right  and  lower  left 
quadrants,  4>f  =  <|>p,  and  negative  in  the  other  two  quadrants, 

<pf  =  - <j>F .  The  normal-force  and  side-force  components  in  the 
four  quadrants  are 


CN.  =  Cp^cos  <pf 
CYi  “  "CF.Sin  *f 

For  panels  on  the  interference  shell 
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The  angle  8_  nTri  is  defined  in  Figure  20. 

f Dir 

The  moment  contributions  due  to  a  panel  are  calculated 
using  Equation  (105)  with  set  equal  to  zero.  The  total 
forces  and  moments  acting  on  the  fins*  and  the  length  of  body 
covered  by  the  interference  shell  are  calculated  oy  summing 
the  individual  panel  forces  and  moments  over  all  of  the  panels. 

The  forces  and  moments  calculated  in  the  manner  just 
described  are  made  up  of  five  components 

(a)  Fin  alone 

(b)  Body  alone 
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(c)  Those  due  to  body  interference  on  fins 

(d)  Those  due  to  fin  interference  on  body 

(e)  Those  due  to  fin-fin  interference 

The  empennage  forces  and  moments  are  the  increments  due  to  add¬ 
ing  the  fins  to  the  body.  To  determine  these  forces  and  moments, 
the  body  alone  (b)  contribution  must  be  subtracted.  These 
forces  and  moments  were  calculated  as  described  in  Section  6.3.1 
and  are  the  sums  of  those  acting  on  the  body  source  panels  which 
coincide  axially  with  the  interference  shell. 

6.3.3  Trailing  Edge  Vortex  Interference 

For  elliptic  stores  with  two  sets  of  fins,  the  effects  of 
trailing  edge  vorticity  from  the  forward  set  of  fins  can  be 
included  in  the  calculation  of  the  forces  and  moments  acting 
on  the  rear  set  of  fins.  Methods  for  doing  this  are  included 
in  prc  gram  VPATHL  of  Reference  7  and  they  have  been  incorpo¬ 
rated  into  the  present  computer  program.  Only  certain  fin 
arrangements  can  be  handled.  These  cases  are  where  each  fin 
in  the  forward  set  lies  in  either  the  yf  -  0  or  zf  =  0  plane 
shown  in  Figure  19.  For  other  fin  arrangements  interference 
of  the  forward  set  of  fins  on  the  rear  set  cannot  be  accounted 
for  because  of  program  limitation. 

The  steps  in  determining  the  interference  effects  on 
the  forward  set  of  fins  on  the  rear  set  are 

(a)  Calculate  vortex  positions  and  strengths  at  trailing 
edge  of  forward  set  of  fins. 

(b)  Calculate  the  vortex  paths  from  the  trailing  edge 
of  the  forward  set  of  fins  to  the  root-chord 
leading  edge  of  the  aft  set. 

(c)  Calculate  the  vortex  induced  velocities  at  the 
control  points  on  the  rear  fins  and  interference 
shell. 
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The  vortex  strengths  and  positions  at  the  trailing  edge 
of  the  forward  set  of  fins  are  calculated  using  the  method 
described  in  Appendix  B  of  Reference  7.  This  method  requires 
the  spanwise  load  distribution  on  each  fin.  They  can  be 
calculated  from  the  individual  panel  forces  Cp^  determined 
by  the  method  described  in  the  preceeding  section.  The  forces 
on  the  panels  in  £  chordwise  row  are  summed  and  divided  by 
the  panel  width  to  determine  the  force  per  unit  span.  The 
vortex  strengths  and  positions  are  calcualted  from  the  span- 
wise  load  distribution  as  described  in  Appendix  B  of 
Reference  7. 

The  second  step,  (b)  above,  is  to  integrate  the  vortex 
paths  from  these  initial  positions  to  the  root-chord  leading 
edge  of  the  rear  set  of  fins.  The  vortex  strengths  remain 
constant  over  this  region.  The  body  shape  is  elliptical  and 
may  be  expanding  or  contracting.  Also,  the  ratio  of  hori¬ 
zontal  semi-axis  length  to  vertical  semi-axis  length,  a/b, 
need  not  be  constant. 

The  method  used  for  the  vortex  tracking  is  that  described 
for  program  VPATHL  in  Appendix  L  of  Reference  7.  It  makes  use 
of  crossflow  plane  solutions  obtained  from  slender-body  theory. 
Because  of  the  linear  nature  of  the  problem  it  is  possible  to 
superimpose  different  solutions. 

Fundamentally,  at  given  axial  stations,  velocity  components 
are  computed  in  the  crossflow  plane  at  the  points  occupied 
by  a  finite  set  of  external  vortices  in  the  presence  of  the 
elliptic  body.  In  this  calculation,  mutual  interaction  between 
the  vortices,  the  interaction  between  the  vortices  and  the 
body,  the  effects  of  local  angles  of  attack  and  sidewash,  and 
changes  in  body  shape  are  included.  Once  the  lateral  velocity 
components  are  known,  the  vortex  locations  at  the  next  axial 
station  downstream  can  be  determined  by  means  of  an  integra¬ 
tion  scheme.  The  vortices  move  from  one  axial  station  to  the 
next  in  accordance  with  the  flow  angle  calculated  from  the 
lateral  velocity  components. 
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The  crossflow  plane  solution  used  to  account  for  the  local 
angles  of  attack  and  sidewash  is  given  in  Section  6  of 
Appendix  I,  Equations  (1130)  and  (1132).  The  angles  a  and  $ 
appearing  in  the  equations  are  calculated  in  the  present  work 
from  the  velocities  given  by  Equation  (87)  of  this  report 
with  vs^  and  ws^,  set  to  zero. 

Changing  body  shape,  that  is  expansion,  contraction, 
and/or  changing  a/b  ratio,  is  accounted  for  by  using  the  method 
of  Section  8  of  Appendix  I  of  Reference  7,  Equation  (1136). 

One  change  has  been  made  to  allow  for  the  varying  a/b  ratio. 

In  Reference  7,  the  axial  rate  of  change  of  body  cross- 
sectional  area.  S',  is  written  assuming  a  constant  a^b  ratio. 

In  the  present  work  the  following  expression  is  used. 


S' 


ir 


da 
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It  is  obtained  by  differentiation  of  the  following  expression 
for  the  cross-sectional  area  of  an  ellipse. 


S  =  Trab 


(112) 


Finally,  the  velocities  induced  at  an  external  vortex 
by  all  the  other  external  vortices  and  all  of  the  image 
vortices  inside  the  body  are  calculated  using  the  method  of 
Section  5  of  Appendix  I  of  Reference  7,  Equation  (1124). 

The  velocities  described  in  the  preceeding  three  para¬ 
graphs  are  calculated  and  summed  up  at  each  external  vortex 
location.  The  vortices  are  then  moved  to  the  next  axial 
station  using  flow  angles  calculated  from  the  velocities. 

This  process  is  repeated  until  the  vortices  have  been  tracked 
to  the  root-chord  leading  edge  of  the  aft  set  of  fins. 

The  velocities  at  the  control  points  on  the  aft  set  of 
fins  and  the  interference  shell  can  now  be  calculated.  For 
all  control  points  the  vortices  are  assumed  to  be  at  the 


100 


»  -  Hr 


*•«$***  +*f 


'I'*?’  O  m  -i>": _ : _ L_ 


same  y8>z8  locations  as  they  were  at  the  root-chord  leading 
edge.  The  velocities  Vg^  and  Wg^  in  Equation  (87)  due  to  N 
vortices  are  given  by 
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The  notation  used  above  is  that  used  in  Section  5  of  Appendix 
I  of  Reference  7. 


7 .  TRAJECTORY  CALCULATION 


The  complete  derivation  of  the  equations  of  motion  of 
a  rigid  body  with  mass  and  inertia  asymmetries  was  presented 
in  Appendix  II  of  Reference  18.  Since  the  publication  of 
that  report  an  error  in  the  derivation  has  been  found  which 
affects  all  of  the  equations  of  motion.  The  equations  in 
the  appendix  changed  by  correcting  the  error  are  Equations 
(11-15)  through  (11-18),  (11-32),  and  (11-40)  through  (11-43)  .  The 
corrected  derivation  is  presented  in  Appendix  B  of  the  present 
report.  The  following  paragraphs  describe  the  derivation  and 
the  equations. 

Generally,  the  body,  a  store  in  the  present  application, 
will  have  axes  of  geometrical  asymruicry  about  which  the  store 
forces  and  moments  are  determined.  By  mass  asymmetry  it  is 
meant  that  the  store  center  of  mass  does  not  have  to  lie  at 
the  origin  of  this  coordinate  system.  By  inertia  asymmetry 
it  is  meant  that  the  principal  axes  of  inertia  of  the  store 
do  not  coincide  with  the  geometric  axes  of  symmetry  so  that 
the  tensor  of  inertia  possesses  off-diagonal  terms,  the  pro¬ 
ducts  of  inertia. 

In  the  derivation  of  Appendix  B  the  two  coordinate  systems 
discussed  earlier  in  conjunction  with  the  velocity  field  cal¬ 
culation,  Section  6.1,  are  used.  They  are  shown  in  Figure  14. 
The  x,y,z  coordinate  system  is  fixed  in  the  store  and  rotates 
with  the  store.  The  x-axis  lies  along  the  store  longitudinal 
axis  and  is  positive  forward,  the  y-axis  is  positive  to  the 
right,  and  the  z-axis  is  positive  downward.  These  axes  coin¬ 
cide  with  the  geometric  axes  of  symmetry  with  the  origin  at 
the  point  about  which  the  moments  and  products  of  inertia  and 
the  aerodynamic  moments  are  calculated.  The  equations  of 
motion  are  derived  in  this  coordinate  system  since,  by  allowing 
the  coordinate  system  to  rotate  with  the  store,  the  time 
derivatives  of  the  moments  and  products  of  inertia  do  not  appear 
in  the  rotational  equations  of  motion.  The  equations  of  motion 


102 


in  thia  coordinate  system  are  given  by  Equation  (B-15)  and 
(B-40)  of  Appendix  B. 

Since  the  store  coordinate  system  in  which  these  equations 
of  motion  are  written  does  not  allow  the  store  position  and 
orientation  to  be  determined  relative  to  the  parent  aircraft, 
a  coordinate  system  is  introduced  which  is  fixed  in  the  fuse¬ 
lage  nose,  the  system  of  Figure  14.  In  the  present  work 

it  is  assumed  that  the  parent  aircraft  is  flying  at  constant 
velocity,  constant  angle  of  attack,  and  constant  flight  path 
angle  relative  to  the  horizontal  so  that  this  system  is  a  non¬ 
rotating  system.  As  such,  it  can  be  regarded  as  an  inertial 
system  and  the  store  motion  relative  to  this  moving  system 
calculated. 

The  system  is  fixed  in  the  fuselage  nose,  and 

coincides  with  the  xF,yF,zF  system  of  Figure  5,  with  £  forward 
along  the  fuselage  longitudinal  axis,  n  laterally  to  the  right, 
and  c  vertically  downward.  In  order  to  determine  the  orienta¬ 
tion  of  the  store  with  respect  to  these  axes,  a  system  of 
angles  must  be  introduced.  These  angles  are  also  shown  in 
Figure  14  and  consist  of  three  rotations  in  the  yaw, 
pitch,  0;  roll,  sequence.  In  order  to  determine  the  time 
histories  of  these  angles,  differential  equations  expressing 
them  as  a  function  of  the  store  rotational  velocities,  p,  q, 
and  r,  are  required.  These  are  given  by  Equation  (B-l)  of 
Appendix  B.  The  rotational  velocities,  p,  q,  and  r,  are 
shown  in  Figure  14.  The  velocity  about  the  x-axis  is  p  and 
is  positive  right-wing  down.  Those  about  the  y-  and  z-axes 
are  q  and  r,  respectively  and  are  positive  nose  up  and  nose 
to  the  right. 

The  final  differential  equations  which  are  to  be  integrated 
to  determine  the  store  position  and  orientation  relative  to  the 
inertial  system  fixed  in  the  fuselage  are  Equations  (B-l) , 

(B-16)  through  (B-18),  and  (B-41)  through  (B-43) .  The  latter 
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six  equations  are  obtained  from  Equations  (B-15)  and  (B-40) 

by  using  the  direction  cosines  relating  the  two  coordinate 

••  ••  •• 

systems,  Equation  (B-2) .  In  the  equations,  £,  n,  and  c  are 

the  accelerations  of  the  store  moment  center,  the  point  about 

which  the  moments  and  products  of  inertia  are  taken,  in  the 

inertial  coordinate  system.  The  velocities  in  this  system 
•  •  •  »  •  « 
are  £,  n,  and  c  and  the  quantities  xQ,  y  ,  and  zQ  are  these 

velocities  resolved  into  the  store  body  coordinate  system. 

That  is 


(114) 


where  [A] '  is  the  transpose  of  the  matrix  given  by  Equation 
(B-2) .  The  location  of  the  store  center  of  mass  relative 
to  the  store  moment  center  in  store  body  coordinates  is  x,y,zT. 
The  three  rotational  velocities  in  the  store  body  coordinate 
system  are  p,  q,  and  r  and  the  corresponding  accelerations 
are  p,  q,  and  r.  The  quantities  .Tvv,  I„„,  and  I  are  the 

XX  yjr  Z  Z 

moments  of  inertia  and  I  ,  I  ,  and  I  are  the  products  of 

y  z  xz  xy 

inertia.  They  are  defined  by  Equation  (B-36) . 

The  remaining  parameters  in  the  equations  of  motion  are 
the  store  mass,  in,  and  the  forces,  Fx,  Fy,  and  Fz,  and  moments, 
M  ,  M  ,  and  M  ,  acting  on  the  store  which  are  positive  in  the 

X  y  Z 

positive  x,  y,  z,  p,  q,  and  r  directions,  respectively.  The 
three  forces  are 
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In  these  equations  g  ,  g  ,  and  g  are  the  component  of  the 

x  y  * 

gravitational  acceleration  vector  g  in  the  store  body  coordi¬ 
nate  system.  Consider  the  following  sketch  which  shows  the 


Yf  J 


.horizontal 


fuselage 


fuselage  flying  at  velocity  V^,  angle  of  attack  ap,  and  flight 
path  angle  The  gravitational  force  acts  normal  to  the 

local  horizontal  so  that  in  the  inertial  coordinate  system 


g  =  -g  sin  (ap  +  Yf)e^  +  g  cos  (ap  +  Yf)e^ 


(116) 


Then 


-g  sin  (ap  +  yf) 


g  cos  (ap  +  Yf ) 


(117) 


The  second  terms  in  Equation  (115)  are  the  aerodynamic 
forces  acting  on  the  store.  The  axial-force  coefficient,  CA, 
is  specified  and  is  positive  in  the  negative  x  direction.  The 
side-force  and  normal-force  coefficients,  Cy  and  CN,  are 
calculated  by  the  methods  described  in  Section  6  and  are  the 
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total  forces  acting  on  the  store  and  empennage (s) .  If  the 
circular  store  option  is  used,  the  forces  are  calculated  using 
the  methods  of  Section  6.2.  For  the  elliptic  store,  the 
methods  of  Section  6.3  are  used. 

The  third  term  in  the  F  equation  is  a  thrust  force, 

Ffc,  positive  in  the  positive  x  direction.  This  force  acts 
along  the  store  longitudinal  axis  and,  in  the  computer  program, 
is  specified  as  a  function  of  time  from  the  beginning  of 
separation. 

The  last  terms  in  Equation  (115)  are  ejector  force  terms. 
The  ejector  model  in  the  computer  program  allows  the  pylon  or 
rack  to  have  one  or  two  ejector  feet.  The  force  produced  by 
each  foot  is  specified,  in  the  program,  as  a  function  of  time 
or  stroke  length.  The  ejector  model  in  the  computer  program 
allows  the  store  to  move  relative  to  the  ejector  feet  during 
separation  so  that  forces  in  all  three  directions  and  moments 
about  all  three  axes  can  be  produced. 

The  parameters  q«>s  and  SR  are  the  dynamic  pressure  and 
reference  area  used  in  nondimensionalizing  the  aerodynamic 
forces.  In  the  present  work  they  are  taken  as  the  store 
free-stream  dynamic  pressure  and  the  store  maximum  cross- 
sectional  area 


q„  =  ~  p  vf  (118) 

*00  7  roo  oo 

s  s 

Sp  =  irrf  or  -rrab  (119) 

R  max 

where  r _ _  is  the  maximum  radius  of  a  circular  store  and  a  and 

max 

b  are  the  semi-axes  of  an  elliptic  store.  The  density  pro  is 
assumed  constant  at  the  value  for  the  aircraft  flight  altitude 
at  time  t  =  0  and  the  store  free-stream  velocity  is  given  by 
Equation  (97) . 
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The  three  moments  about  the  store  moment  center  are 


Mx  =  +  m(V'  '  V’  +  Me 

S  X 

M  =  q  S_£_C  +  m(g  z  -  g  x)  +  M 

y  R  R  m  ^x  yz  e„ 

3  y 

=  q  SiC  +  m (g  x  -  g  y)  +  M 
z  R  R  n  ^y  e 

s  J  z 

The  first  term  on  the  right-hand  side  of  each  of  the  above 

expressions  is  the  aerodynamic  moment.  The  moment  coefficients, 

C»,  C  .  and  C  are  calculated  by  the  methods  of  Section  6.  The 
t  m  n 

reference  length,  £R,  is  taken  as  the  maximum  store  diameter 
for  a  circular  store  and  the  diameter  of  a  circle  with  the 
same  area  as  the  frontal  area  of  an  elliptical  store.  The 
second  terms  are  moments  produced  by  the  gravitational  force 
when  the  location  of  the  store  center  of  mass  does  not  coin¬ 
cide  with  the  center  of  moments.  The  third  terms  are  the 
moments  produced  by  the  force  exerted  by  the  ejector  feet. 

In  order  to  integrate  the  equations  of  motion  given  by 
Equations  (B-16)  through  (B-18)  and  (B-41)  through  (B-43)  and 
the  equations  determining  the  angular  orientation.  Equation 
(B-l) ,  initial  conditions  of  the  variables  must  be  specified. 

These  are  the  position  of  the  store  center  of  moments  (S,n ,  C), 

•  •  • 

the  translational  velocity  of  this  point  (£,n,C),  the  angular 
velocities  about  the  three  store  axes  (p,q,r),  and  the  initial 
orientation  angles  ('F,0,4>).  With  these  specified  the  equations 
can  be  integrated  with  the  aerodynamic  forces  and  moments 
being  recalculated  at  each  point  in  the  trajectory. 


8.  COMPARISONS  WITH  EXPERIMENTAL  DATA 


This  section  of  the  report  will  present  comparisons 
between  the  results  predicted  using  the  present  computer  pro¬ 
gram  and  experimental  data.  The  comparisons  will  serve  to 
assess  the  accuracy  of  the  methods  described  in  the  pro¬ 
ceeding  sections  for  predicting  flow  fields,  store  loading 
distributions,  store  forces  and  moments,  and  store  traject  xes 
The  data  used  for  the  comparisons  were  obtained  during  wind- 
tunnel  tests  conducted  during  the  present  investigation  and 
during  the  investigation  documented  in  Reference  1.  These 
test  programs  and  the  data  obtained  are  described  in  detail 
in  Reference  6.  Before  presenting  the  comparisons,  the 
wind-tunnel  models  will  be  described. 

8.1  Wind-Tunnel  Model  Description 

The  basic  wing-fuselage  combination  used  in  the  test 
program  is  shown  in  Figure  21.  It  is  made  up  of  nose  N-^ , 
body  B2,  and  wing,  W.  This  configuration  is  designated  N-j^W. 
To  obtain  data  in  the  vicinity  of  a  circular  body  alone,  nose 
was  used  in  conjunction  with  body  B^  shown  in  Figure  22  to 
form  configuration  N^B^.  A  noncircular  body  alone  configura¬ 
tion,  N^B^A^ ,  was  also  tested.  This  was  formed  using  nose  N^ 
shown  in  Figure  23  with  the  noncircular  fuselage  adapter  A^ 
shown  in  Figure  24  attached  to  body  B^ .  Data  comparisons  will 
be  presented  for  all  of  these  basic  configurations. 

Comparisons  will  also  be  presented  for  two  noncircular 
wing-fuselage  combinations.  These  are  designated  N^B^A-j  and 
N3b?wa4.  The  first  configuration  is  built  up  from  components 
described  in  the  preceeding  paragraph.  The  second  configura¬ 
tion  has  A^  replaced  by  duct  assembly  A4  shown  in  Figure  25. 

The  parent  aircraft  model  has  provisions  for  mounting 
pylons  on  the  fuselage  centerline  and  left  wing  1/3  semispan 
or  2/3  semispan  positions.  Figure  26  shows  these  positions. 
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The  centerline  of  the  wing  pylons  is  located  at  the  40  percent 
chord.  The  centerline  of  the  fuselage  pylon  is  located  at 
fuselage  station  (FS)  19.42. 

Two  different  pylon  configurations  were  tested.  The 
double-wedge  pylon  (P2)  shown  in  Figure  27  was  tested  at  the 
three  pylon  positions.  The  model  component  designations  are 

(?2 ^ C '  ^P2^1/3'  an<*  ^P2^2/°*  T^ese  Pylons  are  located  as 
shown  in  Figure  26.  The  swept  pylon  (P^)  was  tested  only  at 
the  1/3  seiuispan  position,  ^3)  1/3*  This  pylon  is  shown  in 
Figure  28. 

Some  tests  were  conducted  with  the  triple  ejector  rack, 

T,  shown  in  Figure  4  attached  to  (P2^i/3*  When  attached,  the 
rack  centerline  shown  in  Figure  4  is  aligned  with  che  pylon 
centerline  shown  in  Figure  27. 

Dummy  ogive-cylinder  stores  were  used  in  conjunction  with 
pylon  (p2^1/3  and  rac,c'  T*  The  store  is  shown  in  Figure  29. 
When  attached  to  (p2^i/3'  the  store  midpoint  is  aligned  with 
the  pylon  centerline  shown  in  Figure  27.  When  used  with  the 
rack  T,  the  stores,  SDQC2  and  SD0C3'  where  attached  at  the 
two  shoulder  locations  with  the  midpoints  aligned  with  the 
rack  centerline  shown  in  Figure  4. 

8 . 2  Flow  Fields 

Comparison  will  now  be  presented  between  flow  fields  pre¬ 
dicted  using  the  present  methods  and  data  from  Reference  6. 

A  number  of  comparisons  will  be  shown  for  parent  aircraft 
configurations  built  up  from  the  circular  fuselage,  wing  com¬ 
bination,  NjB2W.  A  few  will  be  presented  for  the  noncircular 
fuselage  configurations.  The  comparisons  to  be  presented  will 
start  with  the  simplest  parent  aircraft  configuration  and 
increase  in  complexity. 

Figure  30  presents  a  comparison  of  the  calculated  radial 
velocity,  Vg/V^,  about  configuration  with  experimental 


data.  Three  theoretical  curves  are  shown;  the  linear  theory 
solution  using  source  panels,  the  modified  linear  theory  solu¬ 
tion  using  source  panels,  and  the  modified  linear  theory 
solution  using  line  sources.  The  modified  linear  theory  solu¬ 
tions  account  for  the  nonlinear  shock  wave  shape  as  described 
in  Section  4.1.  The  modified  linear  theory  results  predict 
the  location  at  which  the  body  first  influences  the  flow  field 
quite  accurately  but  do  not  show  the  sharp  jump  in  velocity 
as  shown  by  the  experimental  data.  The  theory  shows  a  gradual 
build  up.  This  is  the  nature  of  both  solutions.  The  velocity 
induced  at  the  Mach  cone  location  is  zero  and  increases  with 
distance.  The  source  panel  method  produces  a  saw-toothed 
curve  compared  to  the  line  source  method.  This  behavior  is 
due  to  the  axial  density  of  the  source  panels  as  is  shown  in 
Figure  31. 

In  this  figure  results  are  shown  for  two  different 
axial  spacings  of  source  panels.  The  solid  curve  is  the  result 
of  panels  0.5  inch  long  in  the  axial  direction  and  the  dashed 
curve  is  for  panels  2.0  inches  long.  As  can  be  seen,  increas¬ 
ing  the  axial  density  of  the  panels  smooths  the  calculated 
radial  velocity  distribution.  Thus,  to  obtain  smooth  velocity 
fields  the  axial  length  of  the  source  panels  modeling  the  body 
must  be  as  small  as  possible. 

Figure  32  presents  another  comparison  for  the  radial 
velocity  induced  by  configuration  NjB^  The  radial  distance 
to  the  traverse  location  is  5.86  inches.  Results  for  both 
linear  theory  and  modified  linear  theory  are  shown.  The  non¬ 
linear  method  predicts  the  shock  location  much  better  than 
linear  theory.  The  agreement  between  experiment  and  theory., 
is  as  good  as  that  shown  in  Figure  30  for  rfi  *  2.87  inches. 

Figure  33  presents  a  comparison  of  the  upwash  velocity, 
positive  in  the  negative  zp  direction  of  Figure  5,  5.86  inches 
below  the  centerline  of  N^B^  at  an  angle  of  attack  of  5°. 
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The  agreement  between  experiment  and  modified  lineary  theory 
is  about  the  same  as  it  was  at  0°  angle  of  attack. 

The  next  two  figures  present  comparisons  between  experi¬ 
ment  and  the  modified  source  panel  linear  theory  for  the  flow 
field  about  a  noncircular  body,  in  this  case  In 

Figure  34  the  upwash  velocity  is  shown  for  two  axial  traverses 
below  the  body  centerline.  A  sketch  of  the  cross  section  is 
shown  in  the  figure.  The  traverses  are  below  the  flat  bottom 
of  the  body.  The  agreement  for  the  noncircular  body  is  as 
good  as  was  obtained  for  the  circular  body.  The  roughness 
seen  in  the  calculated  results  can  be  reduced  by  increasing 
the  density  of  the  source  panels  in  the  axial  direction. 

Figure  35  presents  comparisons  for  both  upwash  and  side- 
wash  for  an  axial  traverse  at  yF  =  -2.49  in.  and  zF  =  4.36  in. 
Again  the  comparisons  are  quite  good  except  immediately  behind 
the  nose  shock  wave  influence,  xp  =  -5  inches.  Negative  3ide- 
sh  velocity,  Vg,  is  in  the  negative  yp  direction. 

The  next  nine  figures  will  present  comparisons  of  flow 
field  calculations  with  experimental  data  for  the  circular 
fuselage,  wing  combination  Nj^W.  In  the  calculations  the 
line  singularity  method.  Section  3.2.1,  is  used  to  model  the 
f u  .  age.  The  wing  shock  wave  system  i.‘s  located  in  the  manner 
described  in  Section  4.3. 

Figure  36  presents  results  at  the  wing  one-third  semi¬ 
span  position  for  Mtt  =  1.5  and  =  0°.  The  vertical  location 
relative  to  the  local  wing  chord  is  =  1.37  inches.  This 
position  corresponds  to  the  centerline  of  the  store  shown  in 
Figure  26.  The  xp,yF,Zp  coordinate  system  is  shown  in  Figure  5. 
The  velocities  Ug,  Vg,  and  Wg  are  positive  in  the  -Xp,  +yF, 
and  -Zp  directions.  As  can  be  seen  in  Figure  36,  all  three 
velocity  components  are  well  predicted  by  the  theory.  The 
wing  leading-edge  influence,  xp  *  -17.75  inches,  and  the 
trailing-edge  influence,  xf  =  -26.75  inches  are  positioned 
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properly.  In  the  work  of  Reference  1,  the  wing-fuselage 
velocity  field  calculated  between  the  wing  leading-edge  and 
trailing-edge  influences  was  smoothed.  This  is  no  longer 
being  done  so  that  the  theoretical  curves  shown  are  the  actual 
calculated  results. 

Figures  37  and  38  present  similar  results  for  two  other 

zp  locations/  2.12  and  4,37  inches.  Only  Vs  and  Wg  velocities 

are  shown.  The  agreement  between  data  and  theory  for  U_  is 

s 

similar.  For  both  zp  locations  the  comparisons  between 
experiment  and  theory  are  quite  good.  As  zp  increases,  the 
predicted  location  of  the  wing  leading-edge  effect  moves 
slowly  aft  of  that  shown  by  the  data. 

In  Figure  39  results  are  presented  for  zp  =  1.37  inches 
at  the  wing  two- thirds  semispan  location,  yp  =  -8.0  inches. 

Again,  the  comparisons  are  quite  good. 

Similar  results  are  shown  in  Figure  40  for  the  flow  field 
under  the  fuselage  centerline.  The  zp  location,  2.87  inches, 
is  1.2  inches  below  the  bottom  of  the  fuselage.  The  compari¬ 
sons  between  data  and  theory  are  good.  Some  unexplained 
scatter  occurs  in  the  data.  The  wing  effect  which  occurs  at 
Xp  =  -19  inches  is  well  predicted. 

^  * 

Figures  41  and  42  present  comparisons  between  experiment  * 

*  . 

and  theory  for  =  2,0  and  2.5  The  location  is  at  the  wing 

one-third  semispan,  yp  =  -4  inches,  and  zp  =  1.37  inches.  At  * 

both  Mach  numbers,  the  predicted  wing  leading-edge  influence 

is  downstream  of  that  shown  in  the  data.  At  M  =2.5  the 

00 

magnitude  of  the  jump  is  overpredicted.  Except  in  the  vicinity 
of  the  wing  leading-edge  influence,  the  velocities  are  well 
predicted. 

Figures  43  and  44  present  results  for  Mw  =  1.5  with  the 
wing-fuselage  combination  at  5°  angle  of  attack,  aF  =  5°. 

At  zp  =  1.37  inches,  Figure  43,  the  wing  leading-edge  influence 
is  predicted  downstream  of  that  shown  by  the  data.  The  sidewash 
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velocity  is  overpredicted  in  this  region.  Otherwise,  the 
velocities  are  well  predicted.  The  same  holds  true  at  zp  =  4.37 
inches,  Figure  44,  except  that  the  sidewash  velocity  magnitude 
is  well  predicted  in  the  wing  leading-edge  influence  region. 

The  next  three  figures  will  present  comparisons  of  flow 
field  predictions  with  experimental  data  for  the  circular 
fuselage,  wing  combination  N^B2W  with  the  double  wedge  pylon, 

P2,  shown  in  Figure  27  mounted  at  the  wing  one-third  semispan 
location,  yw  =  -4  inches. 

Figure  45  presents  results  for  Mot  =  1.5  and  ap  =  0°. 

The  vertical  location  relative  to  the  local  wing  chord  is 
zp  =  1.37  inches.  Both  pylon-on  and  pylon-off  data  are  shown 
in  the  figure  to  illustrate  the  large  effect  the  pylon  has 
on  the  flow  field.  Theoretical  curves  are  also  shown  for  both 
cases.  For  all  three  velocity  components,  the  theory  predicts 
the  pylon  effect  quite  well.  Downstream  of  xp  =  -22  inches, 
the  pylon-on  calculation  becomes  quite  rough.  This  is  caused 
by  wing-pylon  mutual  interference.  If  more  constant  u-velocity 
panels  than  the  present  program  is  capable  of  handling  could 
be  used  on  the  wing  and  pylon,  this  roughness  could  be 
minimized.  These  panels  account  for  the  wing-pylon  interference. 

Figure  46  presents  similar  pylon-on  results  at  zp  =  2.87 
inches.  Two  sets  of  data  are  shown  and  give  some  indication 
of  the  data  repeatability.  They  were  taken  during  two  dif¬ 
ferent  test  programs.  The  theory  again  predicts  the  data  quite 
well.  The  pylon  influence  begins  at  x  =  -21.5  inches. 

In  Figure  47  results  are  presented  for  Mm  =  1.5  and 
ct„  =  5°  for  z-  -  1.37  inches.  In  general  the  agreement  is  not 

J?  r 

quite  as  good  as  in  Figure  45  for  ap  =  0°  but  the  trends  are 
predicted.  The  experimental  data  are  not  nearly  as  smooth  as 
at  ap  =  0°.  The  roughness  in  the  calculation  downstream  of 
xp  =  -23  inches  is  again  due  to  wing-pylon  interference. 
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Figures  48  and  49  present  comparisons  between  experiment 
and  theory  for  the  case  of  the  TER  shown  in  Figure  4  added 
to  the  pylon  at  the  wing  one-third  semispan  location.  The 
Mach  number  is  1.5  and  the  angle  of  attack  is  0°.  In  the 
theoretical  calculations,  the  small  rack  pylon  was  modeled 
by  extending  the  wing  pylon  through  the  rack  body  as  shown  in 
the  sketch  at  the  top  of  Figure  48(a).  The  thickness  distri¬ 
bution  for  the  portion  inside  the  rack  body,  the  cross-hatched 
area,  was  set  to  zero.  Constant  u-velocity  panels  are  laid 
out  on  this  portion,  however,  to  cancel  the  sidewash  velocity 
on  this  portion  of  the  rack  body.  Note  that  the  velocity  scale 
in  Figures  48  and  49  has  been  changed  from  that  used  in  the 
previous  figures. 

An  examination  of  the  two  figures  for  the  two  zp  positions, 
2.12  and  3.62  inches,  shows  that  the  general  trends  shown  by 
the  data  are  predicted  quite  well  by  the  theory.  The  theoreti¬ 
cal  curves  are  quite  rough  because  of  the  mutual  interference 
being  accounted  for  by  the  u-velocity  panels.  A  finer  panel 
layout  should  smooth  this  out.  The  experimental  data  also 
exhibit  roughness.  As  the  complexity  of  the  aircraft  configura¬ 
tion  increases,  many  more  shock-wave  interactions  and  reflec¬ 
tions  occur  producing  a  very  complicated  flow  field. 

Figure  50  presents  a  comparison  between  experiment  and 
theory  for  an  aircraft  configuration  with  inlets  and  fuselage 
centerline  pylon.  The  Mach  number  is  1.5  and  the  angle  of 
attack  is  5®.  The  agreement  between  experiment  and  theory 
is  quite  poor.  The  data  show  the  inlet  shock  wave  at  xp  =  -15 
inches  while  the  theory  predicts  it  at  xp  =  -17  inches.  It  is 
apparent  considerably  more  work  has  to  be  done  to  provide  a 
theory  which  will  predict  the  flow  field  for  as  complex  an 
aircraft  configuration  as  this.  There  is  also  experimental 
evidence  that  viscous  effects  are  influencing  the  flow  field. 
This  will  be  shown  in  the  next  figure. 
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Figure  51  shows  four  sets  of  experimental  data  taken  at 
the  same  yF,zF  location  as  the  data  shown  in  Figure  50  but 
at  otp  ■  0°  and  with  the  pylon  removed.  The  circles  are  data 
for  the  circular  fuselage ,  wing  combination  and  the  squares 
are  for  the  noncircular  fuselage,  wing  combination.  As  can 
be  seen  the  fuselage  shape  has  little  effect  on  the  data. 

The  remaining  two  sets  of  data  have  the  inlets  present  on 
the  configuration  and  are  for  two  different  inlet  velocity 
ratios  Vj/V^.  This  ratio  was  determined  during  the  tests  of 
Reference  6  by  measuring  pressures  in  the  duct.  The.  method 
of  determining  and  varying  the  ratio  is  discussed  in  Section 
3.3  of  Reference  6. 

Figure  51  shows  that  by  reducing  the  velocity  ratio  from 
0.8  to  0.4  the  inlet  shock-wave  location  is  moved  forward  from 
Xp  «  -16  inches  to  xp  =  -15  inches.  The  other  significant 
fact  to  be  pointed  out  about  the  two  sets  of  inlet-on  data  is 
that  at  Xp  =  -24  inches  there  are  large  differences  in  the 
velocities.  This  location  is  8  inches  downstream  of  the  inlet 
cowl  lip  (see  the  sketch  at  the  top  of  Figure  50)  where, 
except  for  viscous  effects,  the  two  flows  should  be  nearly 
equal.  Viscosity  appears  to  have  a  large  effect.  Vortices 
from  the  cowl  lip  or  inlet  corner  are  probably  present. 

8.3  Store  Load  Distributions 

This  section  will  present  comparisons  of  store  longi¬ 
tudinal  load  distributions  predicted  using  the  present  theory 
and  measured  during  the  test  programs  described  in  Reference  6. 
The  data  to  be  compared  with  were  obtained  using  the  ogive- 
cylinder  pressure  distribution  model  shown  in  Figure  52. 
Complete  pressure  distributions  were  obtained  by  rolling  the 
model  through  360°  in  10°  increments.  The  forces  per  unit 
length  were  obtained  by  a  circumferential  integration  of  the 
pressures.  The  theoretical  calculations  were  made  using  the 
circular-store  option,  that  is,  the  body  is  modeled  by  line 
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sources,  sinks,  and  doublets.  The  method  is  described  in 

Section  6.2.1.  The  store  coordinate  system  (x  ,y  ,z  )  is 

s  s  s 

shown  in  Figure  13  along  with  the  positive  directions  of  the 
normal-force,  CN,  and  side-force,  C^,  coefficients.  The 
reference  area  used  in  forming  the  coefficients  is  the  store 
frontal  area. 

Figures  53  through  57  present  comparisons  of  calculated 
load  distributions  with  experimental  data  for  the  circular 
fuselage,  wing  combination  J^E2W.  Plotted  in  the  figures 
are  the  force  distributions  per  inch  against  distance  from 
the  store  nose  divided  the  store  length,  6.375  inches.  At 
the  top  of  each  figure  is  a  sketch  showing  the  store  position 
relative  to  the  wing,  wing-pylon,  fuselage,  or  fuselage-pylon 
immediately  above  the  store.  The  store  positions  in  the 
carriage  position  on  the  pylons  are  shown  in  Figure  26.  At  the 
one-third  semispan  location,  yp  =  -4  inches,  the  store 
longitudinal  axis  is  at  zp  =  1.37  inches.  Below  the  fuselage 
centerline,  it  is  at  zp  =  2.83  inches.  Data  could  not  be 
obtained  in  these  positions  since  the  store  had  to  be  rolled. 
The  minimum  zp  positions  where  data  were  taken  are  1.47  and 
2.94  inches. 

Figure  53  presents  results  at  the  wing  one-third  semispan 
position,  yp  =  -4  inches,  for  =  1.5  and  ap  =  0°  at 
zp  =  1.47  inches.  In  all  comparisons  the  store  is  at  the  same 
angle  of  attack  relative  to  the  free  stream  as  the  parent  air¬ 
craft.  Figure  53(a)  shows  the  normal  force  distribution.  Two 
sets  of  data  are  shown.  They  were  taken  during  two  different 
test  programs  and  indicate  the  repeatability  of  the  data. 

From  x  /l  =  0  to  about  0.25  the  down  load  is  due  to  the  flow 
s  s 

field  produced  by  the  wing-fuselage  combination.  The  down  load 
is  underpredicted  by  the  theory.  As  was  seen  in  the  flow-field 
comparisons,  the  wing  shock  wave  location  given  by  the  theory 
is  slightly  downstream  of  that  shown  by  the  data.  This  causes 
the  underprediction. 
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Beginning  at  x  //  =0.25  the  effect  of  the  store  noBe 

5  S 

shock  wave  reflecting  off  the  wing  is  seen.  As  can  be  seen, 
the  theory  predicts  the  effect  quite  well.  Without  the  wing 
image  store  used  in  the  theory  to  account  for  the  reflection, 
the  predicted  result  would  be  approximately  that  shown  by  the 
dashed  curve.  The  theory  predicts  a  down  load  from 
xs/£s  *  0.7  to  0.75  not  shown  by  the  data. 

Figure  53  (b)  shows  the  comparison  for  the  side-force 

distribution.  The  agreement  between  data  and  theory  is  quite 

good.  The  store  shock  wave  reflection  off  the  wing  does  not 

effect  the  side  force  to  any  extent.  The  ref  1  "'ctio:.  off  the 

fuselage  shows  a  small  effect.  The  effect  is  seen  in  the  theory 

beginning  at  x  =  0.6. 

s  s 

Figure  54  presents  comparisons  for  the  same  conditions 

except  the  store  is  at  zp  =  2.87  inches,  two  store  diameters 

below  the  position  shown  in  Figure  26.  The  agreement  is  not 

so  good  as  that  shown  in  the  preceeding  figure.  The  data 

show  the  wing-fuselage  effect  to  begin  at  x_/£  =  0.15  while 

s  s 

the  theory  predicts  it  at  x  /£  =  0.25.  The  location  of  the 

s  s 

shock  wave  reflection  off  the  wing  is  well  predicted  at 
xs/ls  =  °*7, 

A  comparison  for  the  store  under  the  fuselage  centerline 
is  presented  in  Figure  55.  The  theory  predicts  the  wing  and 
reflected  shock  effects  at  about  the  same  location,  x  /£  =  0.3 

b  S 

for  zF  =  2.94  inches  and  xg/£s  =  0.5  for  zp  =  3.59  inches. 

This  is  downstream  of  that  shown  by  the  data.  The  positive 
load  downstream  of  the  down  load  is  not  predicted  by  the 
theory . 

Figures  56  and  57  present  results  for  the  store  at  the 
one- third  semispan  position  for  =  1.5  and  ap  =  5°.  For 
both  zp  locations,  1.47  and  2.87  inches,  the  trends  shown  by 
the  data  are  well  predicted.  There  are  some  regions  where 
the  theory  does  not  predict  details  shown  by  the  data. 
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The  next  four  figures  present  comparisons  of  load  distri¬ 
butions  on  the  store  for  parent  aircraft  configurations  which 
include  a  double  wedge  pylon,  P2  shown  in  Figure  27,  mounted 
at  either  the  wing  one-third  semi span  location  or  under  the 
fuselage  centerline.  Figure  58  shows  results  at  the  one- 
third  semispan  location  for  =  1.5,  otp  =  0°,  and  zp  =  1.47 
inches.  There  is  a  0.1  inch  gap  between  the  top  of  the  store 
and  the  bottom  of  the  pylon.  Data  and  theoretical  results 
for  configurations  with  and  without  a  pylon  are  shown  in  order 
that  the  large  effect  of  the  pylon  on  the  store  load  distri¬ 
butions  can  be  seen.  The  pylon-off  data  and  theory  were 
previously  shown  in  Figure  53.  A  large  effect  of  the  pylon  is 
seen  on  the  normal-force  distribution.  Figure  58(a),  while 
little  effect  is  seen  on  the  side-fcrce  distribution,  Figure 
58(b).  In  general,  the  pylon  effects  are  fairly  well  pre¬ 
dicted  by  the  theory. 

Figure  59  presents  results  for  the  store  one  diameter 
below  the  carriage  position  un  the  pylon,  zp  =  2.12  inches. 

The  agreement  between  data  and  theory  is  quite  good. 

The  effect  of  adding  a  pylon  below  the  fuselage  centerline 
is  shown  in  Figure  60.  Again  both  pylon-on  and  pylon-off  data 
and  theory  are  shown  to  illustrate  the  large  effect  of  the 
pylon  on  the  normal-force  distribution.  The  theory  predicts 
the  trends  shown  by  the  data  quite  well  while  missing  the 
magnitudes  of  the  peaks. 

Figure  61  shows  comparisons  for  the  store  below  the  pylon 
at  the  wing  one-third  semispan  location  for  *  1.5  and  ap  =  5°. 
The  store  is  at  zp  =  1.47  inches.  This  location  is  shown  in 
the  sketch  in  Figure  58.  The  agreement  is  not  as  good  as  that 
shown  in  Figure  58  for  ap  =  0°. 

The  results  shown  in  Figures  53  through  61  are  typical 
of  those  obtained  for  other  flight  conditions  for  similar 
parent  aircraft  configurations.  It  was  shown  in  Section  8.2 
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that  as  the  aircraft  configuration  became  more  complex,  the 
flow  field  predictions  were  not  as  accurate.  This  is 
reflected  in  a  loss  of  accuracy  in  the  load  distribution 
calculations. 

All  of  the  comparisons  presented  have  used  the  line- 

T- 

singularity  method  to  model  the  store.  During  the  test 
program  described  in  Reference  6  pressure  data  were  taken  on 
an  elliptic  store.  Only  eight  pressures  were  measured  around 
■  the  circumference  of  the  store.  This  was  insufficient  to 

determine  accurately  the  load  distributions  in  a  complicated 
nonuniform  flow  field  but  yielded  quite  accurate  results  for 
the  store  in  a  uniform  flow.  Some  comparisons  will  now  be 
presented  in  order  that  the  accuracy  of  the  force  and  moment 
calculation  method  used  for  the  elliptic-store  option  can 
be  examined  for  the  store  in  a  uniform  flow  field.  This  method 
j-  is  described  in  Section  6.3.1.  The  store  is  modeled  by  source 

l  panels  on  the  body  surface  rather  than  by  line  singularities 

|  on  the  body  axis. 

f. 

>  The  elliptic  store  model  on  which  the  pressures  were 

measured  is  the  body  of  the  store  shown  in  Figure  62.  The 
wing  and  tail  fins  were  not  present.  The  comparisons  to  be 
shown  are  for  M  *1,5  with  the  store  at  a  combined  angle  of 
attack  of  10°  and  at  various  roll  angles.  The  combined  angle 
of  attack  is  the  angle  between  the  store  longitudinal  axis 
and  the  free- stream  velocity  vector.  Positive  roll  is  right 
wing  down. 

Figure  63  shows  the  ratio  of  local  static  pressure,  p, 
to  free-stream  static  pressure,  pa,  plotted  against  distance 
from  the  store  nose  divided  by  store  length.  For  this  store 
is  *  6.0  inches.  Data  and  tneory  are  shown  for  the  leeward, 
top,  and  windward,  bottom,  sides  of  the  store.  Except  near 
the  nose,  the  theory  predicts  the  pressures  very  well. 
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Figures  64  through  66  present  comparisons  of  load  distri¬ 
butions  on  the  elliptic  store  at  three  roll  angles,  <f>  =  0°, 

s 

45°,  and  90°.  As  the  store  rolls  the  coordinate  system  shown 
in  Figure  13  and,  therefore,  the  direction  of  positive  forces 
rolls  with  the  store.  At  $  =  0°  the  major  axis  of  the 

ellipse  is  horizontal.  At  all  three  roll  angles  the  theory 
predicts  the  load  distributions  very  well  except  near  the 
base  where  flow  separation  may  be  occurring. 

8.4  Store  Forces  and  Moments 

Comparisons  between  experiment  and  theory  of  the  overall 
forces  and  moments  acting  on  a  store  will  be  presented  in  this 
section  of  the  report.  The  first  comparisons  to  be  shown 
are  for  a  finless  ogive-cylinder  store  in  the  presence  of  the 
circular  fuselage,  wing  combination  with  and  without  a  pylon. 
The  wing-fuselage  combination,  Nj^W,  is  shown  in  Figure  21 
and  the  pylon  is  shown  in  Figure  27.  The  pylon  locations 
when  attached  to  the  aircraft  are  shown  in  Figure  26.  The 
calculated  forces  and  moments  were  obtained  using  the  circular 
store  option  described  in  Section  6.2.1.  Following  this, 
results  will  be  presei  ted  for  two  different  stores  in  uniform 
flow;  that  is,  no  parent  aircraft  is  present.  These  stores 
have  fins  and  the  forces  and  moments  were  calculated  using 
the  elliptic  store  option  described  in  Section  6.3. 

The  store  coordinate  system  (x  ,y  , z  )  is  shown  in 

s  s  s 

Figure  13.  The  positive  directions  of  the  forces  and  moments 
measured  and  calculated  in  that  system  are  also  shown  in  the 
figure.  As  the  store  pitches,  yaws,  and  rolls,  the  coordinate 
system  moves  with  the  store. 

The  next  seven  figures.  Figures  67  through  73,  present 
results  for  the  store  shown  in  the  following  sketch.  It 
has  an  ogive  nose  and  a  cylindrical  afterbody.  The  figures 
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present  the  forces  and  moments  as  a  function  of  Azp.  The 
Azf  =  0  positions  at  the  three  pylon  locations  are  shown  in 
Figure  26  and  correspond  to  the  carriage  positions  on  the 
pylons.  The  store  coordinate  system  (x  , y.z_)  is  kept 
parallel  to  the  fuselage  system  (xF,yF,zF)  shown  in  Figure  5. 
The  reference  area  used  in  forming  the  coefficients  is  the 
store  frontal  area,  0.4418  in  ,  and  the  reference  length  in 
the  moment  coefficients  is  the  maximum  diameter,  0.75  in. 

The  point  about  which  the  moments  are  taken  is  shown  in  the 
sketch. 

Figures  67  through  69  present  results  for  the  store  in 
the  presence  of  the  wing-fuselage  combination  without  a  pylon. 
The  Mach  number  is  1.5  and  the  angle  of  attack  is  0°.  Figure 
67  presents  results  at  the  wing  one-third  semispan  location, 
yF  =  -4  inches.  Results  under  the  fuselage  centerline,  yF  =  0 
inches,  are  shown  in  Figure  68,  and  results  at  the  wing  two- 
thirds  semispan  location,  yF  =  -8  inches,  are  shown  in  Figure 
69.  In  general,  the  trends  shown  by  the  data  are  predicted 
fairly  well.  The  poorest  agreement  between  data  and  theory 
occurs  at  the  two-thirds  semispan  location,  Figure  69.  At 
both  wing  stations,  the  normal  force  is  not  well  predicted 
as  the  carriage  position  is  approached,  AzF  -*■  0. 

12.1 


Figure  70  presents  results  at  the  one-third  semispan 
for  a„  =  5°.  The  store  is  also  at  this  angle  of  attack  rela- 

r 

tive  to  the  free  stream.  Comparing  these  results  with  those 
in  Figure  67  for  ap  =  0°  shows  poorer  agreement  for  ap  =  5°. 

Comparisons  with  a  pylon  present  either  at  the  wing  one- 

third  semispan  location  or  under  the  fuselage  centerline 

are  presented  in  Figures  71  through  73.  The  results  shown  in 

Figure  71  are  for  M  =1.5  and  a„  =  0°  at  the  one-third  semi- 

co  r 

span  location.  The  trends  shown  by  the  data  are  predicted 
quite  well  by  the  theory  except  for  CN  and  Cm  as  Azp  0. 

Figure  72  presents  comparisons  for  the  fuselage  centerline 
location  at  M  =1.5  and  a„  =  0°.  The  trends  shown  by  the 
data  are  well  predicted  by  the  theory.  At  large  Azp  the  theory 
predicts  a  return  to  the  free-stream  values  of  CN  =  Cm  =  0 
sooner  than  seen  in  the  data. 

Comparisons  are  present  in  Figure  73  for  =  1.5  and 
ap  «  5°  for  the  store  under  the  pylon  at  the  wing  one-third 
semispan  location.  The  agreement  between  experiment  and  theory 
is  not  as  good  as  that  shown  in  Figure  71  for  ap  =  0°. 

The  comparisons  presented  in  Figures  67  through  73  used 
the  circular  store  option  to  calculate  the  body  forces  and 
moments.  Section  6.2.1.  It  is  of  interest  to  see  how  the 
elliptic  store  option  predicts  the  store  forces  and  moments. 
This  will  be  done  by  comparing  with  data  obtained  on  two 
different  store  models.  These  data  were  taken  without  a 
parent  aircraft  present;  that  is,  the  store  is  in  the  free 
stream.  The  methods  used  to  calculate  the  forces  and  moments 
are  described  in  Section  6.3. 

The  first  comparisons  are  for  the  store  shown  in 
Figure  74.  It  has  an  ogive  nose,  cylindrical  body,  and  tri¬ 
form  fins.  The  forces  and  moments  measured  on  this  store  at 
=  1.5  are  shown  in  Figure  75  along  with  those  predicted  by 
the  theory.  Data  are  shown  for  two  store  roll  angles,  <j>_  =  0° 
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and  90°.  The  fin  orientations  at  these  two  angles  and  the 
orientations  of  the  store  coordinate  system  are  shown  in  the 
following  sketch.  Positive  directions  of  CN  and  are  also 


shown  in  the  sketch.  Also  presented  in  Figure  75  are  results 
for  <|>g  =  0°  at  negative  ac.  This  corresponds  to  <j>g  =  180°  at 
positive  ac.  This  fin  orientation  is  as  shown  below. 


The  agreement  between  experiment  and  theory  shown  in  Figure  75 

is  quite  good.  Comparisons  of  C^,  Cn,  and  are  presented 

only  for  <j>  =  90°.  These  coefficients  are  zero  for  the  4>  =  0° 

s  s 

and  180°  roll  orientations. 


j  'W.  • 
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Figure  76  presents  comparisons  between  experiment  and 
theory  for  the  elliptic  store  shown  in  Figure  62.  This  store 
has  an  elliptic  body  with  a  monoplane  wing  and  an  interdigitated 
tail.  Comparisons  are  shown  in  Figure  76  for  three  roll  angles, 
<j)  =  o°,  50°,  and  90°.  It  should  be  remembered  that  as  the 

S 

store  rolls  the  coordinate  system  in  which  the  forces  and 
moments  are  measured  rolls  with  the  store  as  shown  in  the  pro¬ 
ceeding  sketches.  Thus,  C.T  and  c  are  zero  for  <j>  =  90°,  Cv 

N  m  s  x 

and  C  are  zero  for  4>  =  0°,  and  C»  is  equal  to  zero  for  <j>  =  0° 

n  s  -c  s 

and  90°.  It  can  be  seen  the  method  of  Section  6.3  predicts 
the  forces  and  moments  quite  well.  Because  of  the  close 
coupling  of  the  wing  and  tail  fins,  the  trailing-edge  vortex 
interference  described  in  Section  6.3,3  has  not  been  included 
in  the  calculations.  In  this  short  of  a  distance,  the  vortices 
are  not  fully  developed. 


8.5  Store  Trajectories 

Comparison  between  theory  and  data  will  now  be  presented 
for  two  store  separation  trajectories.  The  data  are  from 
Reference  6.  The  trajectory  calculation  method  was  described 
in  Section  7  of  this  report.  One  trajectory  uses  a  circular 
fuselage,  wing,  and  pylon  combination  for  the  parent  aircraft. 

The  second  configuration  is  made  up  of  a  noncircular  fuselage, 
engine  inlets,  wing,  and  pylon.  In  order  to  simulate  full- 
scale  flight  conditions  the  wind-tunnel  model  components  shown 
in  Figures  21  through  28  have  been  scaled  up  by  a  factor  of 
twenty. 

For  both  trajectories  the  simulated  altitude  is  40,000  feet, 
the  Mach  number  is  1.5,  and  the  aircraft  angle  of  attack  is  5°. 
The  store  angle  of  attack  at  time,  t,  equal  zero  is  the  same 
as  the  aircraft.  In  both  cases  the  aircraft  flight  path  angle, 
Y^,  is  zero  degrees.  That  is,  the  aircraft  is  in  horizontal 
flight,  not  climbing  or  diving.  At  t  =  0  seconds  the  store  is 
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in  the  carriage  position  on  the  pylon  and  is  not  moving 
relative  to  the  parent  aircraft.  For  both  trajectories  aero¬ 
dynamic  damping  was  set  to  zero.  The  store  thrust  is  zero. 

In  both  trajectories  the  following  ejector  force,  time 
history  is  used  for  0  <  t  <  0.075  seconds. 

F  =  1.41448  x  105t  +  2.68614  x  106t2  -  2.220?  x  108t3 
e 

+  3.52968  x  109t4  -  1.70086  x  1010t5,  lbs 

At  t  =  0  seconds  Fg  =  0  pounds.  The  force  reaches  a  maximum 
of  Fe  2  3,000  pounds  at  t  s  0.03  seconds  and  then  decreases 
to  Fg  s  2,000  pounds  at  t  =  0.075  seconds.  This  force  is 
shown  graphically  in  Figure  24  of  Reference  6.  When 
t  >  0.075  seconds,  Fe  =  0  pounds.  One  ejector  foot  is  used  in 
both  trajectories,  however,  the  location  is  different. 

The  store  center  of  mass  in  both  cases  is  taken  to 
coincide  with  the  moment  centers  so  that  in  the  equations  of 
mot,  ~it\ 


x 


z  = 


0 


and ,  therefore , 


I 

xy 


xz 


0 


The  parent  aircraft  configuration  used  for  the  first 
trajectory  is  shown  in  Figure  77.  The  dimensions  shown  are 
full  scale  in  feet.  The  configuration  consists  of  the  circular 
fuselage,  wing  combi  ion  shown  in  Figure  21  with  pylon 

P3  shown  in  Figure  2  a  attached  at  the  wing  one-third  semispan 
position.  An  ogive-cylinder  store  with  cruciform  tail  fins 
is  located  as  shown.  The  tail  fins  are  oriented  45°  from  the 
vertical  and  horizontal.  Details  of  the  store  in  model  scale 
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are  shown  in  Figure  78.  For  this  trajectory,  the  one  ejector 
foot  strikes  the  store  1.5  feet  full  scale  behind  the  moment 
reference  location  shown  in  Figure  78.  In  both  experiment  and 
theory  the  store  rolling  moment  is  set  to  zero.  The  store 
mass  and  inertia  characteristics  are 

m  =  40  slugs 
Ixx  =  20  slug-ft2 
V  -  -  ’00  slug-ft2 

The  trajectory  results  for  this  case  are  shown  in  Figure 
79.  Figures  79(a),  (b) ,  and  (c)  show  the  location  of  the  store 
center  of  gravity  (moment  reference  location)  relative  to  its 
location  in  the  carriage  position  on  the  pylon.  The  £,n»C 
coordinate  system  coincides  with  the  xF,yp,zp  system  shown 
in  Figure  5.  The  store  moves  aft,  outboard,  and  down  relative 
to  the  aircraft.  The  axial  and  vertical  motion  are  well 
predicted  by  the  theory  while  the  small  outboard  lateral  motion 
is  underpredicted.  Figures  79(d)  and  (e)  show  the  pitch  and 
yaw  motions  of  the  store  relative  to  the  carriage  attitude, 
0=5°,  f  =  0°.  The  theory  underpredicts  both  motions. 

The  parent  aircraft  used  for  the  second  trajectory  com¬ 
parison  is  shown  in  Figure  80.  It  is  configuration  (P2) c* 

These  model  components  are  shown  in  Figures  21,  23,  25,  and  27. 
The  dimensions  shown  in  Figure  80  are  full  scale.  The  elliptic 
store  in  model  scale  is  shown  in  Figure  62.  When  the  store  is 
mounted  on  the  pylon,  the  reference  mark  shown  in  Figure  62 
is  at  the  pylon  50  percent  chord  location.  For  this  trajectory, 
the  one  ejector  foot  strikes  the  store  1.0  foot  behind  the 
store  moment  center,  the  assumed  center  of  gravity.  The  store 
mass  and  inertia  characteristics  are 


m  *  20  slugs 


I  =  20  slug-ft2 

V " '  250  slu9-ft2 

The  trajectory  results  for  this  case  are  shown  in  Figure 
81.  The  axial  motion  is  well  predicted  and  there  is  no  lateral 
motion  since  the  store  is  under  the  fuselage  centerline  and 
the  configuration  is  symmetric.  For  the  first  0.2  second  the 
vertical  motion  is  well  predicted  by  the  theory  but  for 
t  >  0.2  second  the  agreement  is  poor.  This  can  be  explained 
in  the  following  way.  The  initial  vertical  motion  is  due  to 
the  ejector  stroke,  0  <  t  <  0.75  second.  At  the  end  of  this 
time  the  aerodynamic  forces  are  the  only  forces  other  than  the 
gravitational  force  acting  on  the  store.  If  the  pitch  motion 
shown  in  Figure  81(d)  is  examined,  it  is  seen  the  theory 
predicts  the  store  to  pitch  nose  up  while  the  data  show  the 
store  pitching  nose  down.  The  store  angle  of  attack  relative 
to  the  free  stream  is  the  sum  of  aF  and  A0  where,  in  this 
case,  <xF  =  5®.  The  theory  predicts  the  store  to  always  be 
at  positive  angle  of  attack,  A0  >  -5®,  which  causes  it  to 
want  to  lift  toward  the  aircraft.  The  data  show  the  store 
pitching  nose  down  to  negative  angle  of  attack,  A0  <  -5°, 
which  would  fly  the  store  away  from  the  aircraft.  The  poor 
prediction  of  the  initial  pitch  motion  is  not  surprising  in 
view  of  the  poor  velocity  field  prediction  shown  in  Figure  50. 
Those  are  the  velocities  along  the  store  centerline  at  t  =  0 
seconds. 
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9 .  CONCLUDING  REMARKS 


This  report  presents  the  results  of  an  investigation  which 
has  been  conducted  with  the  objective  of  extending  and  improving 
the  supersonic  store  separation  prediction  method  of  References 
1  and  2.  In  this  earlier  work  the  fuselage  and  single  store 
v/ere  limited  to  bodies  of  revolution.  The  store  could  have  one 
set  of  planar  or  cruciform  fins.  Wing  dihedral  and  incidence 
relative  to  the  fuselage  were  not  accounted  for.  An  ejector 
rack  could  not  be  modeled.  The  forces  and  moments  acting  on  the 
separating  store  in  the  nonuniform  flow  field  were  calculated 
using  slender-body  theory  and  wing-fuselage- store  interference 
effects  were  not  accurately  modeled.  All  of  these  shortcomings 
have  been  eliminated  in  the  present  work. 

In  the  present  work  linear,  potential  flow  methods  are  used 
to  model  the  various  aircraft  components.  Nonlinear  corrections 
are  made  to  these  models  to  position  shock  waves  more  accurately. 
Two  fuselage  models  are  available.  A  circular  fuselage  can  still 
be  modeled  by  a  distribution  of  line  sources,  sinks,  and  doublets 
on  the  fuselage  axis.  A  noncircular  fuselage,  including  engine 
inlets,  is  modeled  by  a  distribution  of  source  panels  on  the 
body  surface  which  account  for  both  volume  and  angle  of  attack 
effects. 

A  layout  of  constant  u-velocity  panels  on  the  wing  and  pylon 
(if  present)  model  the  lifting  surfaces.  The  wing  may  have 
dihedral  which  can  vary  across  the  span.  Wing  incidence  rela¬ 
tive  to  the  fuselage  is  included  in  the  wing  boundary  condition. 
Wing  and  pylon  thickness  effects  are  modeled  by  an  additional 
layout  of  souce  panels.  Wing  and  pylon  mutual  interference  is 
fully  accounted  for.  Mutual  interference  between  the  wing  and 
pylon  and  the  fuselage  is  handled  by  constant  u-velocity  panels 
placed  on  an  interference  shell  which  extends  over  part  of  the 
fuselage. 
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An  ejector  rack  model  is  included.  The  rack  body  is 
modeled  as  a  body  of  revolution  by  a  distribution  of  line 
sources,  sinks,  and  doublets  on  the  body  axis. 

Up  to  seven  stores  may  be  included.  They  may  be  modeled 
as  bodies  of  revolution  by  line  sources,  sinks,  and  doublets 
on  the  body  axis  or  as  an  elliptical  body  by  a  distribution 
of  source  panels  on  the  body  surface  to  account  for  both 
volume  and  lift  effects. 

In  conjunction  with  the  analytical  work,  an  extensive 
and  systematic  wind-tunnel  test  program  was  carried  out  to 
provide  flow-field  data  for  checking  and  improving  the 
accuracy  of  the  flow  models  and  the  flow-field  calculation 
methods.  In  addition,  store  axial  load  distributions,  store 
forces  and  moments,  and  store  separation  trajectories  were 
measured  to  provide  data  to  be  used  in  evaluating  the  methods 
used  in  calculating  the  store  forces  and  moments  and 
trajectories.  The  wind-tunnel  test  program  is  described  in 
detail  in  Reference  6. 

Flow-field  comparisons  presented  in  this  report  show 
that  the  flow  fields  produced  by  a  circular  fuselage,  a  non¬ 
circular  fuselage,  a  circular  fuselage  with  wing,  and  a 
circular  fuselage  with  wing  and  pylon  are  quite  accurately 
predicted  by  the  theory.  The  agreement  deteriorates  with 
increasing  Mach  number  and  configuration  complexity.  Compari¬ 
sons  are  also  presented  for  a  configuration  which  has  a  triple 
ejector  rack  attached  to  a  wing  mounted  pylon.  The  flow- fie Id 
prediction  is  fair.  Improved  predictions  could  be  obtained 
if  more  singularities  could  be  used  in  modeling  the  aircraft. 

A  comparison  presented  for  a  complex  aircraft  configuration 
including  engine  inlets  shows  poor  agreement  between  data  and 
theory.  Additional  work  is  required  to  improve  the  inlet 
shock  wave  model.  The  wind-tunnel  data  also  show  significant 
viscous  effects. 


129 


i|i,nnipifW—i|iWH;i  Inu  m  miliir*  >W tt iHiV  r~ '  in;  *»' . 


.  5  * 


i 

i 


/'■vSBi* 


In  the  present  work,  mutual  interference  between  the 
separating  store  and  the  wing-fuselage  combination  is 
accounted  for  by  an  imaging  scheme.  Image  stores  are  placed 
above  the  wing  and  inside  the  fuselage  to  account  for  reflec¬ 
tions  of  the  store  nose  shock  wave  off  of  these  surfaces. 
Comparisons  with  experimental  load  distributions  show  this 
imaging  scheme  to  work  quite  well. 

Two  options  exist  in  the  computer  program  for  calculating 
the  store  forces  and  moments.  The  circular  store  option 
uses  line  sources  and  sinks  and  two  sets  of  line  doublets 
placed  on  the  store  axis  to  model  the  store  body  in  the 
nonuniform  flow  field.  This  method  replaces  the  slender-body 
theory  method  used  in  Reference  1.  The  method  used  to  calcu¬ 
late  empennage  forces  is  unchanged.  It  is  based  on  slender- 
body  theory  with  a  correction  for  aspect  ratio  and  includes 
an  application  of  reverse-flow  theorems.  Comparisons  with 
data  for  store-body  load  distributions  and  total  forces  and 
moments  show  that  the  line-singularity  method  gives  quite 
accurate  results  for  parent  aircraft  configurations  consisting 
of  a  circular  fuselage  with  wing  or  wing  and  pylon. 

The  second  option,  the  elliptic  store  option,  models 
the  store  body  in  the  nonuniform  flow  field  by  a  distribu¬ 
tion  of  source  panels  placed  on  the  body  surface.  The  store 
can  have  two  sets  of  fins.  They  are  modeled  by  distributions 
of  constant  u-velocity  panels  placed  on  the  fins  and  on  inter¬ 
ference  shells  on  the  body  to  account  for  fin-body  interference 
Comparisons  between  data  and  theory  are  presented  for  two 
stores  in  a  uniform  stream,  an  ogive-cylinder  store  with 
triform  fins  and  an  elliptic  store  with  monoplane  wing  and 
interdigitated  tails.  The  theory  predicts  the  forces  and 
moments  for  both  configurations  quite  well. 

Comparisons  between  theory  and  experiment  are  presented 
for  two  store  separation  trajectories.  The  first  is  for  an 
ogive-cylinder  store  separating  from  a  pylon  on  the  wing  of 
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a  circular  fuselage,  wing  combination.  The  theory  predicts 
the  separation  trajectory  well  enough  that  it  could  be  used 
for  preliminary  design  purposes.  The  parent  aircraft  used  in 
the  second  comparison  is  comprised  of  a  noncircular  fuselage 
with  engine  inlets,  wing,  and  fuselage  centerline  pylon.  An 
elliptic  store  with  monoplane  wing  and  interdigitated  tail 
is  separated  from  the  pylon.  The  agreement  between  theory 
and  experiment  is  poor.  The  reason  for  this  is  thought  to 
be  that  the  flow  field  induced  by  the  inlets  is  not  modeled 
accurately.  Additional  work  must  be  done  to  improve  the 
inlet  model.  Additional  comparisons  should  also  be  made  to 
validate  the  method  further. 


>&•'  f?  (b)  Radial  velocity 


Figure  2.-  Flow  field  below  a  circular  arc  body  (L/D=5) 
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Figure  4.-  Triple  ejection  rack,  TER. 
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Figure  12.-  Circular  store  fuselage  imaging. 
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Figure  13.-  Coordinate  systems  fixed  in  separated  store 
and  used  in  force  and  moment  calculation. 
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Figure  14. ~  Coordinate  Sj-^tem  used  in  trajectory 

calculation . 


Figure  16.-  Pressure  distribution  on  inboard 
store  in  presence  of  wing- fuselage ;  Hw  = 
a  =  0°,  yF  =  -4  in.,  z  =  1.47  in. 
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♦F  -  fin  dihedral  angle,  PHIDIH 

♦  -  dihedral  angle  of  right  upper 

1  °r  fin,  PHXPR-PHIDXH 

♦  ,  *>  dihedral  angle  of  left  lower 

fin,  PHIFh-PHIDIH 

♦  “  dihedral  angle  of  right 

3  °r  U  lower  fin,  PHIFU-  -PHIDIH 

* .  _  ■  dihedral  angle  of  left  upper 

4°rD  fin,  PHIFD- -PHIDIH 

6  -  fin  location  polar  angle, 
THETIT 

6.  «  polar  angle  of  right  upper 

1  °r  R  fin,  THETR-THETIT 

0.  *  polar  angle  of  left  lower 

2  °r  L  fin,  TOETL-THETIT 

0,  u  »  polar  angle  of  right  lower 
°r  fin,  THETI-  -THETIT 

#.  .  *  polar  angle  of  left  upper 

fin,  THETD-  -THETIT 


Body  interference  panels: 

^iP.j 


&tf  "  sf,i  "  *f,i-l 
Ayf  “  yf,i  ~  yf,i-l 


sin  8 


2.BIP 


cos  0 


2,BIP 


Vas*  +  Ay2 


Va*j  +  Ayj  J 


these  functions  are  used 
in  the  transformation 
(rotation)  of  the  body 
interference  panel  to  the 
reference  (xs,ys,zs)  or 
<x.f»yf.rf>  system 


Figure  20,-  Geometrical  angles  associated  with  case  involving 
fins  on  body  with  elliptical  cross  section  and 
force  coefficients  associated  with  fins  and 
body  interference  panel. 


igure  23.-  Noncircular  nose 


Figure  24.-  Noncircular  fuselage  adapter. 


Throat  r~  To  match  wing  contour 


Figure  25.-  Duct  assembly  for 


All  dimension* 
in  inches 


Figure  26.-  Ogive-cylinder  force  and  moment  store  in  the 
carriage  position  at  the  fuselage  centerline,  1/3 
semispan,  and  2/3  semispan  pylon  stations. 


All  dimensions  in  inches 


Figure  27.-  Double-wedge  pylon,  P2. 
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Figure  28.-  Swept  pylon,  P, 


Figure  34.-  Upwash  velocity  below  the  centerline  of  configuration  N3B1A3 
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(a)  Axial  velocity. 

Figure  36.-  Flow  field  under  the  wing  at  the  one-third 
semispan  location;  Mffl  =  1.5, 
ap  =  0°,  zp  =  1.37  in. 
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(c)  Upwash  velocity. 
Figure  36.-  Concluded. 
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(a)  Sidewash  velocity 
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(a)  Axial  velocity. 


Figure  40.-  Flow  field  under  the  fuselage  centerline 
M  =  1.5,  a„  =  0°,  z„  =  2.87  in. 
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(b)  Upwash  velocity 


Figure  40.-  Concluded. 
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(a)  Sidewash  velocity. 


Figure  41.-  Flow  field  under  the  wing  at  the  one- third 
semispan  location;  M  =  2.0, 
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(b)  Sidewash  velocity. 
Figure  45.-  Continued. 
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(b)  Sidewash  velocity. 
Figure  46.-  Continued. 
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(b)  Sidewash  velocity. 


Figure  48.-  Continued. 
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(b)  Upwash  velocity 


Figure  51.-  Concluded. 
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(a)  Normal-force  distribution. 

Figure  53.-  Store  load  distributions  under  the  winq  at  the  one- 
third  semispan  location;  Mm  =  1.5, 
ap  -  0°,  Zp  =  1.47  in. 


<b)  Side-force  distribution 


Figure  53.-  Concluded. 
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ZF  “  1.47  in. 


(a)  Normal-force  distribution. 


Figure  58.-  Effect  of  pylon  on  load  distributions  on  store  at 
one-third  semispan  location;  *  1.5, 
ciy  =  0°,  zF  =  1.47  in. 
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(b)  Side-force  distribution 


Figure  58.-  Concluded 


(b)  Side-force  distribution 


Figure  59.-  Concluded. 
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(a)  Normal-force  distribution. 


Figure  61.-  Effect  of  pylon  on  load  distribution  on  store 
at  one-third  semispan  location 
=  1.5,  aF  =  5  ° ,  z-,  =  1 . 4  7  1  n . 


sssure  distributions  on  the  top  and  bottom 
the  elliptic  store?  M  =1.5, 
ctc  =  10°,  <|>s  =  0°. 


229 


Theory 


J _ 1 _ L 

.2  .4  .6 


(b)  Side-force  distribution. 
Figure  65.-  Concluded. 


232 


;%■ 


233 


(a)  Normal-force  coefficient 


AzF,  in. 

(b)  Pitching-moment  coefficient. 


Figure  68.-  Forces  and  moments  on  the  ogive-cylinder  store  under 
the  fuselage  centerline;  =  1.5,  ap  =  0°. 
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Figure  69.-  Forces 
the  two- third 


and  moments  on  the  ogive-cylinder 
semi  span  location;  =  1.5,  otp  -• 


store  at 
0°. 
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(c)  Side-force  coefficient 
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(a)  Normal-force  coefficient. 


Azp,  in. 

(b)  Pitching-moment  coefficient. 


72.-  Forces  <•  .id  moments  on  the  ogive-cylinder  store 
under  the  pylon  at  the  fuselage  centerline; 

Moo  =  !-5'  «F  =  0°. 
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(a)  Normal-force  coefficient. 


forces  and  moments  on  the  ogive-cylinder  store  under 
the  pylon  at  the  one-third  semispan 
location;  =  1.5,  ap  =  5°. 


(a)  Normal-force  coefficient. 
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(c)  Side-force  coefficient. 
Figure  75.-  Continued. 
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(e)  Rolling -moment  coefficient 


Figure  75,-  Concluded. 
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Figure  7c.-*  Forces  and  moments  on  the  elliptic  store  with  wing 
and  tail  in  the  free  stream,  *  1.5. 


(c)  Side-force  coefficient. 


4.0  8.0  12.0  16.0  20.0 

ac,  deg. 

(d)  Yawing-moment  coefficient. 
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(e)  Rolling-moment  coefficient. 


Figure  76.-  Concluded. 


Figure  77.-  configuration  for  trajectory  of  ogive-cylinder 

store  with  cruciform  fins. 
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APPENDIX  A 


NORMAL-  AND  SIDE-FORCE  DISTRIBUTIONS  ACTING  ON  AN 
AXISYMMETRIC  STORE  IN  SUPERSONIC, 

NONUNIFORM  FLOW 

A-l  INTRODUCTION 

A  three-dimensional  method  for  calculating  load  distribu¬ 
tions  and  overall  forces  and  moments  acting  on  a  store  in 
supersonic,  nonuniform  flow  in  the  vicinity  of  the  parent  aircraft 
is  presented.  The  procedure  is  summarized  as  follows. 

The  parent  aircraft  is  flow  modeled  first  and  the  strengths 
of  the  singularities  used  to  model  the  fuselage,  wing,  pylon, 
ejector  rack,  and  other  stores  are  known.  Since  the  separated 
store- to-parent  aircraft  interference  is  accounted  for  by  means 
of  an  imaging  scheme  described  in  Section  4.4,  the  strengths  of 
the  aforementioned  singularities  are  not  affected  by  the 
separated  store.  If  this  interference  were  to  be  included  in 
the  boundary  condition  applied  at  points  on  the  aircraft,  the 
singularity  strengths  would  be  affected  by  store  position.  The 
store  is  immersed  in  the  nonuniform  flow  field  induced  by  the 
parent  aircraft  and  the  store  image (s).  The  store  is  flow 
modeled  in  this  nonuniform  flow  field  and  the  pressures  on  the 
store  are  determined. 

The  next  step  is  to  flow  model  the  store  in  the  nonuniform 
flow  field.  In  order  to  account  for  volume  effects,  supersonic 
line  sources  and  sinks  are  disturbed  along  the  store  centerline. 
Appendix  I  of  Reference  1  describes  a  solution  representing 
volume  effects  cf  an  axisymmetric  body  in  supersonic,  uniform 
flow.  The  solution  can  be  extended  to  account  for  nonuniform 
axial  flow.  The  constants  in  the  linearly  varying  source  or  sink 
strengths  expressions  are  obtained  from  the  flow  tangency 
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condition  for  the  store  aligned  with  the  local  flow  streamlines. 

As  a  consequence,  the  constants  are  related  to  the  product  of 
the  body  slope  measured  in  a  meridional  plane  and  the  local 
axial  velocity  at  a  control  point  on  the  body  surface. 

In  order  to  cancel  effects  due  to  nonuniform  upwash  and 
sidewash,  supersonic  line  doublets  are  distributed  along  the 
store  axis.  The  solution  for  the  cross  flow  problem  associated 
with  an  axisymmetric  body  in  supersonic  uniform  flow  is  also 
given  in  Appendix  I  of  Reference  1.  This  solution  can  be  extended 
to  account  for  varying  upwash  and  sidewash  at  points  along  the 
store  centerline  induced  by  the  parent  aircraft  and  store 
image (s).  The  constants  in  the  linearly  varying  doublet  strength 
equations  are  determined  from  an  application  of  the  flow 
tangency  condition  at  control  points  on  the  surface  of  the  store 
body. 

The  velocities  induced  by  the  parent  aircraft  and  store 
image (s)  are  computed  at  points  along  the  centerline.  Nonlinear 
corrections  are  applied  in  this  process  as  described  in  Section  4 
of  this  report.  These  velocities  are  transferred  to  the  body 
surface  for  use  in  the  flow  tangency  conditions  associated  with 
the  axial  and  crossflow  problemc.  This  constitutes  an  approxi¬ 
mation  to  the  flow  tangency  condition  in  that  the  velocities  to 
be  cancelled  at  the  surface  point  are  taken  equal  to  the  velocities 
on  the  centerline.  Thus  as  far  as  the  singuxarity  strengths 
are  concerned,  no  account  is  taken  of  the  lateral  variation  of 
axial  velocity,  upwash  and  sidewash  in  the  region  occupied  by 
the  store.  The  error  introduced  by  this  approximation  is  a 
function  of  the  extent  of  nonuniformity  of  the  flow  coupled 
with  the  slenderness  of  the  store.  The  procedure  used  to  cal¬ 
culate  velocities  induced  by  the  store  image (s)  is  described  in 
Section  4.4.  The  image  of  the  store  is  represented  by  the  same 
singularities  as  the  real  store.  On  account  of  the  step-by-step 
solution  procedure  to  be  described  below,  the  effects  from  the 
image  store  required  in  the  boundary  condition  of  the  real  store 
are  induced  by  singularities  with  known  strengths. 
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Once  the  strengths  of  the  singularities  representing  the 
store  are  known,  velocities  induced  by  the  singularities  are 
computed  at  the  control  points  on  the  store  surface.  Pressures 
are  computed  from  the  velocities  including  free  stream  components 
using  the  Bernoulli  relationship.  Normal-  and  side-force 
distributions  are  obtained  by  integrating  the  pressures  along 
the  circumference  of  a  store  body  cross  section. 

The  store  modeling  technique  will  be  described  in  detail. 
This  is  followed  by  a  discussion  of  the  pressure  calculation 
method  and  the  procedures  used  to  calculate  the  force  distribu¬ 
tions  and  overall  forces  and  moments  acting  on  the  store. 


* 
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A 


A-2  FLOW  MODEL  FOR  AXISYMMETRIC  BODY  IN 
NONUNIFORM  FLOW 

In  accordance  with  the  earlier  work  (References  1  and  20)  in 
supersonic  store  separation  and  missile  aerodynamics,  bodies  of 
revolution  with  pointed  noses  in  uniform,  supersonic  flow  are 
represented  oy  the  following  sets  of  singularities.  Supersonic 
line  sources  account  for  volume  effects  and  supersonic  line 
doublets  model  effects  of  upwash  and  sidewash.  In  this  report 
the  solutions  to  the  source  and  doublet  problems  are  modified 
to  include  effects  of  nonunifo'-m  flow.  The  modifications  affect 
the  flow  tangency  conditions  as  described  in  the  following 
section. 


A-2.1  Flow  Tangency  Boundary  Condition  at 
the  Body  Surface 

The  sketch  below  shows  a  body  of  revolution  in  a  flow 
field  characterized  by  velocity  components  U,V,W  which  depend  on 
the  axial  coordinate  xB  only.  In  the  application  of  the  method 
to  an  ejected  store,  these  velocity  components  include  effects 
due  to  store  pitch  and  yaw  damping,  transnational  velocities  of 
the  store  c.g.  and  store  angle  of  pitch  and  sideslip  in  addition 
to  the  perturbation  velocity  components  induced  by  the  parent 
aircraft  and  store  image.  The  rectangular  coordinates  :<B,yB,zfi 
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with  origin  at  the  store  nose  and  a  cylindrical  coordinate 
system  xB'rB»9  are  indicated  in  the  sketch.  Also  shown  are  the 
nonuniform  free-stream  velocity  components  U (xfi) ,V (xB) ,W(xB) 
in  the  xB,yB,zB  directions,  respectively.  Perturbation 
velocities  uB,vfi,wB  are  in  the  xg,  rB,  and  6  directions.  In  a 
meridional  plane,  the  body  contour  is  given  by  r_  =  R(x_) .  The 
flow  tangency  condition  given  by  Equation  (4)  in  Reference  1  is 
now  restated  with  the  inclusion  of  variable  velocities  U,  V,  and 
W.  Note  that  the  free-stream  velocity  is  the  resultant  flow 
velocity,  with  components  Us  ,  Vs  ,  and  Ws  ,  as  seen  by  the 

00  CD  00 

store.  In  the  application  to  an  ejected  store,  the  resultant 

velocity  at  infinity  as  seen  by  the  store  includes  effects  due 

to  store  angle  of  pitch  and  sideslip  and  translational  velocities 

of  the  store  c.g.  The  resultant  velocity  has  components  Ue  ,  Vs 

°00  °00 

Ws  .  M  is  related  to  U„  .  The  boundary  condition  states  that 

O.  °00 

the  resultant  flow  vector  at  the  body  is  aligned  with  the  local 
body  slope. 


B  W  V 
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The  equation  governing  the  flow  potential,  <{>,  is  expressed  ir. 
the  cylindrical  coordinate  system  shown  in  the  previous  sketch. 
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Because  Equation  (A-2)  is  linear  in  the  solution  can  be 
composed  as  the  sum  of  three  potentials 


*(VrB'9>  "  *a(VV  +  *o„(xB'rB'e)  +  *c„(xB'rB'0)  (A_3) 

W  V 
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which  satisfies  Equation  (A-2)  .  Potential  <j»  is  the  axial 
potential.  Potentials  <pCw  and  <J)Cv  are  crossflow  potentials 
associated  with  upwash  W  and  siewash  V,  respectively.  In  terms 
of  cylindrical  coordinates,  the  perturbation  velocities  are 
related  to  the  two  potentials  in  accordance  with  the  following 
expressions. 
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Substituting  for  the  perturbation  velocities  in  Equation  (A-l) 
gives  the  following  result. 
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This  expression  for  the  flow  tangency  condition  can  be  split  into 
three  parts  for  bodies  of  revolution. 
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Equation  (A-6a)  is  the  boundary  condition  for  the  axial  flow 

problem  with  the  axial  flow  velocity  U  a  function  of  Xg. 

Similarly,  Equations  (A-6b)  and  (A-6c)  involve  only  crossflow 

velocities  associated  with  upwash  W  and  sidewash  V,  respectively. 

The  sum  of  the  three  equations  yields  the  complete  flow  tangency 

condition,  Equation  (A-5).  Thus,  the  first  term  on  the  right- 

hand  side  of  Equation  (A-3)  together  with  Equation  (A-6a) 

constitute  the  axial  problem.  According  t-  Reference  21,  the 

reference  velocity  for  the  axial  problem  is  the  velocity  far 

from  the  body  parallel  to  the  body  axis  Uro  and  the  associated 

Mach  number  is  M  .  The  second  term  on  the  right-hand  side  of 

«a 

Equation  (A-3)  and  Equation  (A-6b)  represent  the  upwash  problem. 
The  third  term  on  the  right-hand  side  of  Equation  (A-3)  and 
Equation  (A-6c)  are  the  sidewash  problem.  For  these  two  problems, 
the  reference  velocity  refers  to  the  full  free-stream  velocity,  Vw 

and  the  Mach  number  is  M  . 

00 

The  three  problems  will  now  be  discussed  in  detail.  Where 
convenient,  reference  will  be  made  to  the  pertinent  work  in 
Reference  1. 

A- 2. 2  Axial  Flow  Problem 


The  technique  for  modeling  the  volume  effects  due  to  a 
body  of  revolution  in  uniform  flow  is  described  in  Appendix  I, 
Section  1-2.1,  of  Reference  1.  The  analysis  presented  there 
needs  to  be  modified  to  include  the  effect  of  the  variation  of 
the  longitudinal  velocity  component  U  with  coordinate  xB,  the 
distance  along  the  body  centerline. 

The  axial  flow  potential  <b  is  governed  by  Equation  (A-2) 

a 

arter  omitting  the  angular  dependence  term.  The  result  is 


(A-7) 
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In  accordance  with  Section  9.13  in  Reference  21,  the  reference 
flow  associated  with  the  axial  flow  problem  corresponds  to  the 

flow  parallel  to  the  body  axis  at  infinity,  Us  ,  as  seen  by  the 

00 

store.  The  boundary  condition  given  by  Equation  (A-6a)  must  be 
satisfied  at  points  on  the  body  surface.  It  is  rewritten  below. 
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From  Equation  (1-10)  in  Appendix  I  of  Reference  1, 
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Substituting  Equation  (A-9)  into  (A-8)  gives  the  following  general 
result  for  constant  K 
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For  body  slope  dR/dxfi  and  axial  flow  velocity  U  specified  at 
some  xD  station,  Equations  (A-9)  and  (A-10)  describe  the  flow 
field  due  to  one  line  source  resulting  in  a  conical  :flow  field. 
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The  velocities  are  functions  of  parameter  xn/8  r_,  only  and  thus 

O  Ct  D 

are  invariant  along  lines  for  which  x^/r^  are  constant. 

Following  the  analysis  described  in  Reference  21,  the 
axial  flow  around  an  axisymmetric  body  aligned  with  the  flow 
can  be  obtained  from  a  superposition  of  conical  solutions  of 
the  type  discussed  above.  The  composite  source  distribution  for 
the  axial  flow  solution  is  given  by  a  superposition  of  linearly 
varying  line  sources.  In  fact,  the  distribution  corresponds  to 
a  superposition  of  basic  cone  solutions 

f (?)  +  Kc?  +  K-^?  -  ?x)  +  K2(?  -  ?2)...  (A-ll) 

In  this  specification,  constants  and  K2  are  negative  for  a 
typical  nose  shape  and  =  0 ,  =  Xb^  -  3a*Bdi'  ^2  =  XBd2  " 

6arBd2'  etc*  are  the  origins  of  the  line  singularities  expressed 
in  terms  of  body  definition  points  (xB(i,  rBd)  to  be  discussed 
next.  The  constants  are  evaluated  in  a  step-by-step  manner  as 
described  in  Appendix  I  of  Reference  1.  Note  that  the  numerator 
in  Equation  (A-10)  can  be  viewed  as  a  modified  body  slope. 

Thus,  the  axial  flow  analysis  for  a  body  of  revolution  in 
uniform  flow  can  be  directly  adapted  to  the  case  involving  non- 
uniform  flow. 

Let  NXBODY  be  the  number  of  line  sources  with  origins 
equally  spaced  on  the  body  centerline.  The  number  of  control 
points  also  equals  NXBODY.  The  xfi  coordinates  of  the  control 
points  are  half-way  in  between  the  xBd  coordinates  assigned  to 
uhe  body  definition  points.  Line  source  origins  ?Q, ?2» . . ,  are 
related  to  the  body  definition  point  coordinates  as  mentioned 
above.  Geometrically,  the  line- source  origins  are  the  inter¬ 
sections  of  the  Mach  lines  drawn  from  the  body  definition  points 
with  the  body  or  xg  axis.  The  relative  positions  of  the  body 
definition  points,  line- source  origins,  control  points,  and 
Mach  lines  are  shown  in  the  following  sketch. 
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In  accordance  with  Equation  (A-10) ,  the  first  constant  KQ  is 
calculated  from 


(A-12) 

where  xBl  and  rBl  are  the  coordinates  of  the  first  control  point. 
The  n^h  constant  associated  with  the  (n+l)^  line  source,  Kn, 
is  determined  from  the  flow  tangenoy  condition  applied  at  the 
(n+l)th  control  point  (XBn+1*  rBn+1> •  The  boundary  condition 
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includes  the  effects  of  the  line  sources  with  known  constants 
Kk-1  an<^  or*9ins  positioned  ahead  of  the  origin  5n  of  the 

line  source  under  consideration  (1  <  k  £  n) .  These  effects  are 
added  to  the  axial  flow  term  in  the  numerator  of  Equation  (A-10) . 
Thus,  constant  is  given  by 


K  =  =■  f  n  =  1,2,...  (NXBODY-1) 
n  d 
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The  terms  in  the  above  equation  are  specified  as  follows: 
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=  Mach  number  associated  with  flow  velocity 
a  Us  parallel  to  body  axis  at  infinity 
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[Equation  (A-14)  continued] 
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Ufx_  1  =  axial  flow  velocity  along  the  body  centerline 
v  n+1-1  at  xn 
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After  completion  of  the  step-by-step  procedure  for  calculating 
constants  Kq  and  Kn  where  n  =  1 , 2 , . . . (NXBODY-1) ,  the  axial  and 
radial  perturbation  velocities  induced  at  a  field  point  with 
coordinates  (xg,  rfi)  by  the  line  sources  associated  with  the 
axial  flow  problem  are 
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(A-15 ) 


The  upwash  and  sidewash  crossflow  problems  and  the  associated 
line-doublet  solutions  will  now  be  described. 

A-2.3  Crossflow  '  ro:  lems  Due  to  Upwash  and  Sidewash 

The  analysis  presented  here  is  an  extension  to  the  cross- 
flow  analysis  in  Appendix  T  of  Reference  1.  The  added  features 
are : 
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(1)  Separate  sets  of  doublets  for  upwash  and  sidewash 

(2)  Allowance  for  variation  with  distance  along  the  body 
centerline  of  the  upwash  and  sidewash 

A- 2. 3.1  Solution  for  Upwash  Line  Doublets 

The  crossflow  potential  4>Cw  f°r  upwash  must  satisfy  Equation 
(A- 2)  subject  to  the  boundary  condition  given  by  Equation  (A-6b) . 
It  is  shown  in  Reference  1  that  the  solution  can  be  constructed 
from  the  general  (not  the  cone  solution  discussed  previously) 
axisymmetric  solution.  The  result  is  a  linearly  varying  doublet 
solution  for  upwash  with  the  doublet  axis  parallel  to  the  zg 
axis  shown  in  the  sketch  of  Section  A- 2.1. 
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The  reference  flow  Mach  number  is  M  associated  with  the  free 

00 

stream  velocity  Vro  seen  by  the  store.  Factor  3  is  therefore 
specified  as 
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(A-17 ) 


The  crossflow  around  an  arbitrary  body  of  revolution  can 
be  obtained  by  a  superposition  of  doublet  solutions  in  the  same 
manner  used  for  the  axial  flow  problem.  Thus,  the  composite 
doublet  distribution  for  the  upwash  crossflow  solution  is 
obtained  by  generating  a  superposition  of  linearly  varying  line 
doublets 
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dw(E)  =  Kd] 


(5  -  EX>  +  K 


(A-18) 


'Ill 


in  terms  of  constants  K<jw  Q,  K<jw  1 ,  Kjjw  2  and  origins  =  0, 

£>1  ~  “  0rBdi'  C2  “  xBd2  "  0rBd2  of  the  line  singularities. 

These  constants  are  also  evaluated  in  the  same  step-by-step 
manner  employed  in  the  axial  flow  problem  but  with  the  boundary 
condition  specified  by  Equation  (A-6b) .  This  flow  tangency 
condition  must  be  satisfied  at  the  control  points  on  the  surface 
of  the  body  shown  in  the  sketch  of  Section  A-2,2  concerned  with 
the  axial  flow  problem.  Equation  (A-6b)  is  rewritten  as 
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Direct  dif ferentation  of  the  crossflow  potential,  Equation 
(A-16) ,  gives  the  perturbation  velocities  associated  with  the 
line  doublet  for  upwash. 
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Note  that  the  partial  differentials  (representing  perturba¬ 
tion  velocities)  in  Equation  (A-19)  are  given  by  the  first  two  of 
Equation  (A-20)  and  contain  the  factor  cos  9.  Consequently, 
the  cos  0  factor  is  cancelled  in  Equation  (A-19) .  Referring 

again  to  the  sketch  in  Section  A- 2. 2,  the  first  constant  K^.. 

W  f  o 

is  obtained  from  the  boundary  condition  applied  to  the  first 
control  point  (xjj^,  z^)  with  a  substitution  from  Equation  (A-20) 
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Let  NXBODY  be  the  number  of  line  doublets  for  upwash  with 
origins  equally  spaced  on  the  body  centerline.  In  general,  the 
nth  constant  is  explicitly  determined  from  the  boundary  condition 
applied  at  the  (n+l)th  control  point  as  follows 


K-  =  ^  ;  n  =  1,2, . . . (NXBODY-1)  (A-22) 
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[Equation  (A-23)  continued] 
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The  step-by-step  procedure  indicated  by  Equation  (A-22) 

together  with  Equation  (A-21)  provides  the  solution  for  constants 

Kdw  Kdu  where  n  =  1,2,... (NXBODY  -  1).  The  axial,  radial, 
n  #  O  vi  f  n 

and  tangential  perturbation  velocities  induced  by  the  line 
doublets  associated  with  the  upwash  problem  at  a  field-point 
(xB»rB,0)  are 
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A- 2. 3.2  Solution  for  Sidewash  Line  Doublets 


for  this  problem,  the  doublet  axis  is  parallel  to  the  direc¬ 
tion  of  the  axis  shown  in  the  first  sketch  of  Section  A- 2.1. 
The  solution  can  be  related  to  the  solution  for  the  upwash  by 
a  rotation  of  coordinate  axes  yfi  and  zQ  as  shown  in  the  following 
sketch.  It  is  seen  that  the  relation  of  the  y^,z^  or  r ^,0* 
system  to  the  V  vector  is  the  same  as  the  relation  of  the  yB»*E 
or  rB,0  system  to  the  vector  W.  In  fact,  the  former  can  be 
derived  from  the  latter  by  means  of  a  clockwise  rotation  in 


angle  6  through  90°.  As  a  consequence.  Equations  (A-20)  and 
(A-24)  for  perturbation  velocities  due  to  upwash,  replace  cos  0 
with  sin  6,  and  sin  8  with  -cos  6.  The  resulting  expressions 
represent  the  perturbation  velocities  induced  by  sidewash.  The 
values  for  constants  Kd„  „  through  K*.,  are  obtained  from 
Equations  (A-21)  and  (A-22)  replacing  W(xb1)/Voo  and  Wtxg^^J/V^ 
with  VUb^/V^  and  V(xDn+1)/Vro. 

Alternatively,  it  can  be  shown  that  the  cross-flow  potential 
for  sidewash  can  be  related  to  the  general  axisymmetric  solution 
as  follows 


4>c  (xB,rB,9)  *  sin  0 


(A- 25) 


For  the  upwash  problem  the  sin  6  factor  is  replaced  with  cos  8. 

The  potential  for  the  sidewash  problem  nan  then  be  derived  in  the 
manner  used  for  the  upwash  problem  as  summarized  in  Appendix  I 
of  Reference  1. 

The  composite  doublet  distribution  for  the  sidewash  cross- 
flow  solution  is  given  by  a  superposition  of  linearly  varying 
line  doublets 

dv(5)  *  K e  +  K .  <5  -  e.)  +  K.  U  -  5o)  +  ...  (A-26) 

V  °V,o  ^,1  A  v#  2 

in  the  same  manner  used  for  the  upwash  case,  Equation  (A-18) . 

The  first  constant  K<av  is  obtained  from  the  boundary  condition 
given  by  Equation  (A-6c)  which  is  rewritten  below 
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In  accordance  with  the  transformations  mentioned  above,  the 
perturbation  velocities  associated  with  one  line  doublet  for 
sidewash  can  be  obtained  from  Equation  (A-20) .  The  result  is 
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[Equation  (A—  28 )  continued] 
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In  Equation  (A-27) ,  the  partial  differentials  representing 
perturbation  velocities  are  given  by  the  first  two  of  Equation 
(A-28)  and  contain  the  factor  sin  0.  Thus,  the  factor  sin  0 
is  canceled  in  Equation  (A-27) .  The  first  constant  Kdv  is 
obtained  from  the  boundary  condition  applied  to  the  first  control 
point  (xg^rij.)  shown  in  the  sketch  of  Section  A- 2. 2 


cosh 


-1 


B, 


fir 


B, 


where 


XB,  "  ixB/2 


280 


<  XS**  v 


Let  NXBODY  be  the  number  of  line  doublets  for  sidewash  with 
origins  equally  spaced  on  the  body  centerline,.  In  general,  the 
nfc^  constant  is  explicitly  determined  from  the  boundary  condi¬ 
tion  applied  at  the  (n+l)th  control  point  as  follows. 


[Equation  (A-31)  continued] 
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The  step-by-step  procedure  indicated  by  Equation  (A-30) , 

together  with  Equation  (A-29)  provides  the  solution  for  constants 

Kd„  Kd„  „  where  n  -  1,2, . . . (NXBODY-1) .  The  axial,  radial, 
v  f  o  v  t  n 

and  tangential  perturbation  velocities  induced  by  the  line 
doublets  associated  with  the  side-wash  problem  at  a  fieldpoint 
(xB,rB,0)  can  be  determined  from  Equation  (A-24)  with  the 
transformations  discussed  earlier.  The  results  are 
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[Equation  (A-32)  continued] 
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Knowing  the  strengths  of  the  line  sources,  sinks,  and  doublets, 
the  pressure  can  now  be  computed  at  points  on  the  surface  of 
the  store  body.  The  pressures  are  then  integrated  to  yield  the 
force  distribution  on  the  store.  The  technique  used  will  now 
be  described. 


A-3  CALCULATION  OF  LOADS  ON  THE  STORE  IN 
SUPERSONIC ,  NONUNIFORM  FLOW 


This  section  is  concerned  with  the  methods  used  to  compute 
the  pressure  distributions,  load  distributions,  and  forces  and 
moments  acting  on  an  axisymmetric  store  in  nonuniform,  supersonic 
flow.  The  strengths  of  the  singularities  used  to  represent  the 
store  in  such  a  flow  field  have  been  determined  by  the  methods 
described  in  the  previous  section.  Velocities  induced  by  the 
store  singularities  at  points  on  the  store  body  surface  are 
calculated.  The  Bernoulli  velocity-pressure  relationship  is  then 
used  to  compute  pressure  coefficients  at  points  on  the  body  sur¬ 
face.  Integration  around  the  body  circumfernece  yields  the  axial 
distributions  of  normal  and  side  forces  and  the  pitching  and 
yawing  moments.  Finally,  the  axial  distributions  are  integrated 
to  obtain  overall  forces  and  moments. 

In  the  calculation  of  pressures  on  the  store  body  summarized 
above,  velocity  components  induced  by  the  store  image  and  parent 
aircraft  are  not  included.  In  Section  6.2.1  it  was  shown  that 
this  assumption  improves  agreement  between  the  theory  and  data 
for  the  store  pressure  distribution.  Over  the  region  of  the  store 
where  shock  waves  impinge,  this  assumption  is  rationalized  on  the 
basis  of  pressure  bleed-off  by  the  strong  boundary- layer  cross- 
flow.  A  description  of  these  viscous  effects  is  given  in 
Section  5.4  of  Reference  1.  Note  that  the  strengths  of  the 
singularities  modeling  the  store  in  the  nonuniform  flow  field 
are  affected  by  effects  from  the  image  store  (s)  and  parent  air¬ 
craft.  Away  from  the  shock  impingement  regions,  their  effects 
should  probably  be  included  in  the  pressure  calculations.  At 
the  present  time,  only  the  sources  or  sinks  and  the  doublets 
associated  with  the  store  itself  induce  velocities  at  the  points 
on  the  body  surfaces  where  pressures  are  calculated. 
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A- 3.1  Effects  Due  to  Store  Volume 


f  * 


Store  volume  is  represented  by  a  distribution  of  linearly 
varying  line  sources  (supersonic  cone  solutions)  along  the  store 
body  axis.  The  strengths  of  the  line  sources  are  related  to  the 
store  body  surface  slopes  measured  in  a  meridional  plane.  The 
reference  velocity  used  in  the  1  me  source  strength  calculation 
is  the  component,  USoo  measured  along  the  store  axis  of  the  free- 
stream  velocity  as  seen  by  the  store.  This  free-stream  velocity, 
designated  V  ,  includes  effects  of  translation  of  the  store 
relative  to  the  parent  aircraft. 

At  a  point  (x0,yB,zB)  on  the  surface  of  the  store,  the 
perturbation  velocities  induced  by  the  line  sources  representing 
store  volume  are  calculated.  Thus,  at  a  point  (xB,yB,zB)  on 
the  surface  of  the  store,  the  perturbation  velocity  components 
are  expressed  as  follows. 
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(A-33) 


Note  that  these  perturbation  velocities  are  expressed  in  terms  of 
the  axial,  radial,  and  tangential  components  aligned  with  the 
cylindrical  coordinate  system  (xB,rB,0)  indicated  in  the  sketch 
of  Section  A- 2.1.  Also,  the  reference  velocity  is  the  one 
actually  seen  by  the  store,  Vw,  mentioned  earlier.  Terms  sub¬ 
scripted  with  "a"  pertain  to  the  volume  effects  of  the  store 
(axial  solution)  and  are  given  by  Equation  (A-15) .  In  that 
equation,  the  reference  velocity,  U__,  must  be  scaled  up  to  V 

w  g  w 

as  shown  in  Equation  (A-33) . 
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At  *-Ms  stage,  the  perturb'  on  velocities  specified  in 
Equation  (A-33)  include  effe  ...  of  the  store  volume.  However, 
at  the  point  (xB>yB/ zr^  '  the  b°dy,  a  net  flow  normal  to  the 
surface  exists  which  must  be  counteracted  by  additional  singu¬ 
larities  distributed  along  the  store  body  axis.  Their 
contributions  will  be  discussed  in  the  next  section. 

A-3.2  Effects  Due  to  Crossflow 

The  line  sources  distributed  along  the  store  axis  satisfy 
the  condition  of  no  flow  through  the  body  surface  provided  the 
store  is  aligned  with  the  flow.  This  case  is  dealt  with  in  the 
previous  section.  In  general,  the  store  is  not  aligned  and 
therefore  subjected  to  crossflow  as  well.  In  order  to  cancel 
flow  normal  to  its  surface,  line  doublets  are  distributed  along 
the  store  axis  in  addition  to  the  line  sources. 

The  crossflow  is  associated  with  the  nonuniform,  supersonic 
flow  field  induced  by  the  parent  aircraft  in  which  the  store  is 
immersed.  It  includes  effects  of  store  pitch  and  yaw  damping. 

In  addition,  store  translational  velocities  are  accounted  for  in 
the  specification  of  the  crossflow.  All  these  contributions 
are  combined  into  a  crossflow  velocity  vector  whose  variation  with 
distance  xfi  along  the  store  axis  can  be  specified.  This  vector 
can  be  resolved  into  components  along  the  z0  and  yfi  axes  of  the 
store  resulting  in  upwash  component,  W(xfi),  and  sidewash  component, 
V(xfi) .  Subject  to  these  crossflow  velocity  components,  the 
strengths  of  the  line  doublets  for  upwash  and  sidewash  are  then 
determined  using  the  methods  described  in  detail  in  Section  A. 2. 3. 

At  a  point  (xB,yB,zfi)  on  the  surface  of  the  store,  the  per¬ 
turbation  velocities  induced  by  the  upwash  and  sidewash  doublets 
are  calculated  using  Equations  (A-24)  for  uB/dw»  vB,dw'  and 
WB,dw  and  Equation  (A-32)  for  UB,dv»  vB,dV'  and  wB,dv*  They  are 
added  to  the  store  volume  perturbation  velocities  to  give 
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Components  (uB,vfi,wB)  are  the  axial,  radial,  and  tangential 
velocity  components  in  the  cylindrical  coordinate  system  shown  in 
the  sketch  of  Section  A— 2.1.  For  use  in  the  Bernoulli  pressure 
equation,  to  be  described  next,  these  perturbation  velocity 
components  will  be  transformed  into  the  components  aligned  with 
the  rectangular  coordinate  system  (xB,yB,£B).  The  reference 
velocity,  V^,  is  the  free-stream  velocity  as  seen  by  the  store. 

As  mentioned  before  in  connection  with  Equation  (A-33) ,  the 
perturbation  velocities  induced  by  the  line  sources  and  sinks 
associated  with  the  store  volume  are  given  by  Equation  (A-15) , 
but  the  reference  velocity  must  be  scaled  up  to  Vffi  as 
indicated  in  Equation  (A-33) . 


A-3.3  Pressure  Coefficient  Calculation 


At  a  point  (X3»yBfzg)  on  the  surface  of  the  store,  the 
pressure  coefficient  is  calculated  using  the  Bernoulli  relation¬ 
ship*  For  compressible,  isentropic  flow,  the  Bernoulli  equation 
is  given  by 


r 


<  1  + 


v. 


1 


y/y-i 


*\ 

1 


(A-35 ) 


The  ratio  of  specific  heats  for  air  is  y  and  equals  1.4. 
Reference  velocity,  V^,  and  Mach  number,  Mro,  are  as  seen  by  the 
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store  at  infinity.  The  magnitude  of  the  resultant  velocity,  V, 
at  a  point  on  the  body  surface  can  be  obtained  in  the  following 
manner . 

The  perturbation  velocities  are  now  known  at  the  point 
(X3'^B'ZB^  on  surface  of  the  store  through  Equation  (A-34) 

in  terms  of  components  aligned  with  the  cylindrical  coordinate 
system  (xB,rB,0)  attached  to  the  store  as  shown  in  the  sketch  of 
Section  A-2.1  and  repeated  below  for  convenience.  A  transforma¬ 
tion  converts  the  perturbation  velocity  components  of  Equation 
(A-34)  into  components  aligned  with  rectangular  coordinate  system 


(x^,yQ,Z-,)  with  origin  at  the  store  nose  and  directions  along 

id  £5  JtS 

the  major  axes  of  the  store  also  shown  in  the  sketch  above. 


Longitudinal  component  uB  remains  as  is. 


Let  vyfi  and  wZfi 


denote  the  perturbation  velocity  components  in  the  yB  and 
directions,  respectively.  As  such,  all  three  perturbation 


- — t— — — — 

"y  •  *  a-..'  ' 


I 


velocity  components  are  specified  in  a  body  axis  coordinate 
system  (as  opposed  to  the  wind  axis  system  aligned  with  the 
free-stream  direction) .  Components  VyB  and  wZb  are  related  to 


and  wB  as  follows. 


XB  UB 


Yr  vr  w_ 

— —  =  —  sin  0  +  — -  cos  0 

CO  v  00  v  CO 


(A-36) 


ZR  vr  WR 

— cos  0  -  sin  0 

Voo  Vco  Vco 


The  perturbation  velocities  specified  in  the  above  equation 
are  now  specified  in  the  cartesian,  body  coordinate  system. 
The  magnitude  of  the  resultant  velocity,  V,  required  in 
Equation  (A-35)  can  be  expressed  in  terms  of  the  perturbation 
velocities  in  the  body  axis  system  by  the  following  equation. 


cos  a 


■B  . 


sin  3, 


2  2  2 

2w  u  +  v  +  w 

ZR  XR  Yr  ZR 

+  -rr-5-  sin  a  +  — 2 - - B 

Voo  S  v2 

00 


(A-37) 


Angles  a  ,  S  ,  and  a  are  the  included  angle  of  attack,  side- 
C  f  s  s  s 

slip  angle,  and  pitch  angle,  respectively,  associated  with  the 
store.  In  accordance  with  page  5  of  Reference  19,  the  angles  can 
be  related  to  the  components  along  the  body  coordinate  system  of 
the  free-stream  velocity,  V  ,  as  seen  by  the  store.  , 
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sin  as 

sin  BS 


V 


oo 


2  2  2 
sin  a  „  =  sin  a„  +  sin  8 
c ,  s  s  s 


J 


(A-38 ) 


Special  care  must  be  taken  with  the  term  inside  the  square 
brackets  in  Equation  (A-35)  and  raised  to  the  power  y/(Y  ~  D • 

It  is  possible  that  numerical  inaccuracies  can  occur  in  the 
determination  of  the  perturbation  velocities.  In  certain  cases, 
the  inaccuracies  can  cause  the  bracketed  term  to  be  negative 
resulting  in  a  negative  pressure  in  the  flow.  This  is  physically 
impossible  and  the  bracketed  term  must  therefore  be  limited  to 
values  which  are  equal  to  or  larger  than  zero. 


A-3.4  Loading  Distributions  Along  the  Store  Axis 


At  a  given  axial  location,  the  pressure  coefficients  are 
calculated  at  points  on  the  circumference  of  a  cross  section  of 
the  store  body  using  the  method  described  above.  In  terms  of 
length  xn  along  the  store  axis  from  the  nose,  and  the  store 

■D 

frontal  area  as  reference  area,  expressions  for  the  rate  of  change 
of  normal  and  side  force  with  the  axial  length  can  be  written. 

The  normal  and  side  forces  act  along  the  zfi  and  yfi  axes, 
respectively. 


dCN  . 

r(xB) 

dXE 

Sref 

dCy 

„  r(xB) 

dxB 

Sref 

Cpcos  0  d0 


CpSin  0  d0 


(A- 39) 
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In  these  expressions,  angle  6  (radians)  is  as  shown  in  the 
previous  sketch  and  is  positive  clockwise  from  the  positive 


zQ  axis. 


The  distributions  of  pitching  and  yawing  moments  are  given 


x_  -  x  dC„ 
B  m  N 


XB  -  xm  dCY 


(A-40) 


In  Equation  (A-40) ,  d  is  the  maximum  diameter  of  the  store  and 

x  is  the  location  on  the  body  axis  of  the  moment  center, 
m 

Due  to  the  nonlinear  dependence  of  the  pressure  coefficient 
on  the  resultant  velocity,  V,  analytical  integration  of 
Equation  (A-39)  is  laborious  but  not  necessarily  impossible. 

It  is  convenient  to  integrate  Equation  (A-39)  numerically.  At 
a  finite  number  of  discrete  points  on  the  circumference  of  a 
body  cross  section,  the  pressure  coefficient  is  calculated  from 
the  perturbation  velocities  by  means  of  Equations  (A-35)  and 
(A-37).  The  products  Cpcos  0  and  Cpsin  0  are  formed  at  the 
discrete  points.  An  integration  scheme  based  on  a  Simpson's 
rule  takes  in  the  products  and  the  A0  interval  between  points  and 
provides  the  numerical  answers. 

A-3.5  Overall  Forces  and  Moments  Acting  on  the  Store 

The  overall  force  and  moment  coefficients  can  be  determined 
by  multipling  Equations  (A-39)  and  (A-40)  by  Axfi  and  summing  over 
the  length  of  the  body.  The  interval  Axfi  is  given  by  the  spacing 
between  the  axial  stations  at  which  the  loading  distributions 
are  calculated.  Thus,  the  normal-force,  side-force,  pitching- 
moment  and  yawing-moment  coefficients  are  then  given  by  the 
following  expressions,  respectively. 


r:,s 


where  N  is  the  number  of  axial  stations  on  the  body  centerline 
at  which  the  circumferential  pressure  distributions  are  cal¬ 
culated. 


APPENDIX  B 


EQUATIONS  OF  MOTION  OF  A  RIGID  BODY  WITH 
MASS  AND  INERTIA  ASYMMETRIES 


B-l  INTRODUCTION 

The  purpose  of  this  appendix  is  to  present  the  derivation 
of  the  equations  of  motion  of  a  rigid  body,  the  store,  with 
mass  and  inertia  asymmetries.  By  mass  asymmetry  we  mean  that 
the  center  of  mass  is  not  at  the  origin  of  the  coordinate  system 
fixed  in  the  body  and  used  in  the  equations  of  motion.  Generally 
the  external  store  will  have  axes  of  geometrical  symmetry  which 
will  be  used  in  determining  the  forces  and  moments  acting  on 
it,  but  the  center  of  mass  in  certain  applications  will  not  lie 
on  these  axes  of  symmetry.  By  inertia  asymmetry,  we  mean  that 
the  principal  axes  of  inertia  do  no  coincide  with  the  geometric 
axes  of  symmetry  of  the  store.  For  these  geometric  axes,  the 
tensor  of  inertia  will  thus  possess  off-diagonal  terms,  the  so- 
called  products  of  inertia. 
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B-2  AXIS  SYSTEMS  AND  EULERIAN  ANGLES 

In  the  derivation  of  the  equations  of  motion  consider  first 
an  inertial  set  of  axes,  <£,  n>  C) ,  fixed  in  the  aircraft  which 
is  in  rectilinear  flight  at  uniform  velocity.  Then  a  system  of 
axes  (x,  y,  z)  fixed  in  the  store  is  defined.  The  x  axis  is 
forward  along  the  store  longitudinal  axis,  the  y  axis  is  posi¬ 
tive  to  the  right  looking  forward,  and  the  2  axis  is  vertical 
downward  so  that  a  right-handed  system  exists.  The  £,  n,  c 
system  is  an  airplane  fixed  system  with  £  measured  forward 
along  the  fuselage  longitudinal  axis,  n  laterally  to  the  right 
and,  5  vertically  downward.  The  origin  is  at  the  fuselage  nose. 

In  order  to  determine  the  orientation  of  the  store  with 
respect  to  the  £,  n,  t  axes,  a  system  of  angles  will  be  intro¬ 
duced.  The  system  of  angles  to  be  used  is  that  shown  in  Figure 
14.  The  origins  of  the  inertial  and  store  axis  systems  are 
taken  at  the  same  point  to  show  the  angular  orientation.  The 
system  consists  of  three  rotations  in  the  yaw,  pitch,  roll 
sequence.  This  order  is  important.  The  three  angles  will  be 
designated 


y  yaw 
0  pitch 
#  roll 
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B-3  DIFFERENTIAL  EQUATIONS  FOR  THE  ORIENTATION  ANGLES 


In  specifying  the  position  and  orientation  of  the  store 
relative  to  the  aircraft,  we  will  determine  the  coordinates 
£,  n»  and  c  of  the  origin  of  the  x,  y,  z  coordinate  system  and 
the  angular  quantities  V,  0,  and  4.  The  equations  of  motion 
will,  however,  be  written  in  terms  of  the  following  quantities 
and  their  time  derivatives:  £,  n,  t,  P,  q,  r.  Here  p,  q, 
and  r  are  the  components  of  the  store  angular  velocity  about  the 
x,  y,  and  z  axes,  respectively.  A  set  of  differential  equations 
is  required  for  Y,  0,  and  4  in  terms  of  p,  q,  and  r  which  are 
to  be  integrated  along  with  the  equations  of  motion  to  keep 
track  of  the  angular  orientation  of  the  store.  These  equations 
are  from  Reference  16. 


'i'  =  (q  sin  4  +  r  cos  4) /cos  0 

&  =  q  cos  4  -  r  sin  4 

4  *  p  +  q  sin  4  tan  0  +  r  cos  4  tan  0 


(B-l) 


With  initial  values  of  the  orientation  parameters  0,  and  4, 
the  foregoing  equation  can  be  integrated  to  obtain  the  variations 
of  these  quantities  with  time. 

From  the  values  of  0,  and  4  it  is  possible  to  determine 
the  orientation  of  the  x,  y,  z  axes  with  respect  to  the  £,  n»  C 
in  the  usual  direction  cosine  way.  That  is 


K 

X 

n 

=  [A] 

y 

c 

z 
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v/here  (Reference  1C) 


cos  0  cos  ¥  sin  $  sin  0  cos  ¥  cos  $  sin  0  cos  ¥ 

-  cos  $  sin  ¥  +  sin  $  sin  ¥ 

[A]  *  cos  0  sin  ¥  sin  $  sin  0  sin  ¥  cos  $  sin  0  sin  ¥ 

+  cos  $  cos  ¥  -  sin  $  cos  ¥ 

-sin  0  sin  4>  cos  0  cos  $  cos  0 

_  — 

(B-2) 

These  data  are  also  necessary  in  converting  forces  cal¬ 
culated  in  the  x,  y,  z  system  to  the  £,  n/  £  system. 


B-4  TRANSLATIONAL  EQUATIONS  OF  MOTION 

We  will  now  write  the  equations  of  motion  for  transla¬ 
tion  of  the  origin  of  the  x,  y,  z  body  axes.  The  origin  is 
not  assumed  to  be  at  the  body  center  of  mass.  Consider  the 
following  sketch. 


Let  a  be  the  vector  from  the  origin  of  the  inertial  system 
point  P,  to  the  origin  of  the  body  system,  point  0.  Let  r' 
be  the  vector  from  the  origin  of  the  body  axes  to  the  body 
center  of  mass. 


—  — 


r  =  ex  x  +  ey  y  +  cz  2 


(B-3) 


1  i 


Let  r  be  the  vector  from  the  origin  of  the  inertial  system  to 
the  body  center  of  mass.  Then 


r  =  a  +  r ' 


(B-4 ) 


If  ?  is  the  vector  force  acting  on  the  body  the  equation  of 


translational  motion  is 


d2? 


(B-5) 


This  equation  will  now  be  expressed  in  terms  of  the  force  com¬ 
ponents  (Fx'Fy'Fz)  along  x,  y,  and  z;  the  position  coordinates 


.  ■«  ■;  #*»**■■  ’  -  -  '*•«' ' 


which  from  Equation  (B-5)  is  equal  to 


y-  =  -  f"  F  e  +  F  e  +  F  e  1 
m  L  x  x  y  y  z  zj 


(B-12) 


The  acceleration  of  0  relative  to  P  can  be  expressed  in 
terms  of  the  moving  body  coordinate  system  as  follows: 


d2£ 


=  x_e 


da 


dt 


2  '"ox  ■* o  y 


+  ye  +  z  e  +  w  x  -3-rr 
Jo  y  o  z  dt 


(B-13) 


da 

where  ^  is  the  velocity  of  0  relative  to  P  and  may  be  expressed 
in  terms  of  components  along  the  moving  body  coordinate  system 
x  ,  yQ,  and  zQ.  The  body  angular  velocity  w  is  given  by  Equation 
(B-6 ) .  Therefore 


,2 

"dt2”  =  ^  +  Y°e*  +  V.  +  (q2°  ‘  Y°r)ex 


+  (xQr  -  pzo)ey  +  (pyo  -  xoq)ez 


(B-14 ) 


Combining  Equations  (B-l] )  ,  (B-12),  and  (B-14)  gives 
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The  accelerations  along  the  inertial  axes  are  then 


XQ  +  (qzQ  "  ryQ) 


H  =  [A]  yQ  +  (rxQ  -  pzQ) 


20  +  (pyQ  -  qxQ) 


where  the  transformation  matrix  is  given  by  Equation  (B-2) . 

The  final  equations  of  translational  motion  in  the  inertial 
coordinate  system  are  now  written  using  the  following  direction 
cosine  notation  for  the  [A]  matrix  of  Equation  (B-2) 


cos  ($,x)  cos  (£,y)  cos  (£,z) 
[a]  -  cos  (n/x)  cos  (n»y)  cos  (n#z) 
cos  (C,x)  cos  (c,y)  cos  (C#z) 


't  +  [y  cos  (£  ,z) 
+  [x  cos  ( E,,y ) 


-  z  cos  (4,y)]p  +  [z  cos  (£,x)  -  X  cos  (£,z)Jq 

-  y  cos  (?,x)]  r 

F  1 

=  +  x(q  +  r  )  -  ypq  -  zpr  cos  (5,x) 

F  —  —  2  2  — 

+  ~  xpq  +  y(p  +  r  )  -  zqr  cos  (£,y) 

'Fz  -  -  -  2  2  ’ 

+  — —  xpr  -  yqr  +  z(p  +  q  )  cos  (5,z) 


(B-16 
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L y  cos  (n»z) 
[x  cos  (n,y) 


£  y  cos  U,z) 
[x  cos  (K,y) 


z  cos  (n,y)Jp  +  [z  cos  'r'»x'  "  x  008  (r'»z)J<3 
y  cos  (n#x)]  r 

F„  _  «  2  _  _ 

—■  +  x(q  +  r  )  -  ypq  -  zpr  cos  (n»x) 

-  . 

F 

+  -  xpq  +  y(p2  +  r2)  -  zqr  cos  (n,y) 

m 

'F  _  2  2  I 

+  -  xpr  -  yqr  +  z  (p  +  q  )  cos  (n>z) 


z  cos  U,y)jp  +  [z  cos  (r,,x)  -  x  cos  U»z)J 
y  cos  (C  #x)  j  r 

'F  _  0  2  -  -  1 

+  x(qz  +  r  )  -  ypq  -  zpr  cos  (c,x) 

+  -  xpq  +  yCp^  +  r  )  -  zqr  cos  U,y) 

"*  F 

+  —  -  xpr  +  yqr  +  ¥(p2  +  q2)  cos  (C,z) 

m 


B-5  EQUATIONS  OF  ROTATIONAL  MOTION 


The  equations  of  rotational  motion  will  be  written  in 
terms  of  body  axes  because  the  moments  and  products  of  inertia 
are  not  functions  of  time  in  these  axes.  Also,  the  moments 
due  to  aerodynamic  forces  are  calculated  in  this  axis  system 
which  is  an  additional  convenience.  Consider  a  general  body 
with  axes  x,  y,  z  having  its  origin  at  point  0  and  an  inertial 
system  £,  n ,  £  having  its  origin  at  point  P.  Let  p  be  the 
position  vector  of  an  element  of  mass  in  the  inertial  system 
and  p'  be  the  position  vector  in  the  body  system. 


Let  tp  be  the  torque  of  the  forces  acting  on  the  body  with 
respect  to  point  P.  Then  the  law  of  the  conservation  of  angular 
momentum  in  inertial  coordinates  is 


* 

Lp  =  St 


(B-19) 


where  the  angular  momentum  vector  is  given  by 


&  =  P  x 

P  dt 


(B-20) 


the  integral  being  taken  over  the  body  mass. 

For  the  store  the  origin  0,  of  the  body  system,  will  generally 
be  offset  from  the  body  center  of  gravity  by  position  vector 


2 r  Sr*- 


r’,  and  the  tensor  of  inertia  for  the  x,  y,  z  axes  will  generally 
contain  products  of  inertia.  The  torque  will  be  known  with 
respect  to  the  origin  0,  and  the  angular  velocity  is  given  by 
Equation  (B-6)  in  terms  of  components  p,  q,  and  r  in  body  axes. 

We  will  now  transform  Equation  (B-19.'  to  body  axes  in  terms  of 
the  foregoing  quantities.  The  angular  momentum  becomes 

Sp  =  j  (a  +  p • )  X  A.  (S  +  p ’ )dm  (B-21) 


Let 


so  that 


sP  =  t> 


f|  dm  +  |  s*  X  ||  dm  +  ac 


(B-22) 


and 


d. 

"dt 


-ax 


1 0—it 

+  l 
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*|dm 

at 


_ o 

dt 


H  *!?■ 


dm  x 


d2a 

dt2 


(B- 


Certain  simplifications  arise  in  the  foregoing  expression 
from  the  propr  rty  of  the  center  of  mass 

(B-24 ) 

(B-25) 

(B-26) 
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|  ( p 1  -  r 1 ) dm  =  0 

so  that 

* 

p ’ dm  =  r ’ m 

and 


Tl’TWfPPIffWf 


Also 


Equation  (B-29)  for  the  conservation  of  angular  momentum  with 
respect  to  origin  O  has  the  same  form  as  that  with  respect  to 
origin  P  plus  an  additional  term. 

With  regard  to  the  three  terms  of  Equation  (B-29) ,  the 
torque  will  have  three  components  along  the  x,  y,  *.  axes 


t,  =  e  M  +  e  M  +  e  M, 
o  xx  y  y  z  z 


(B-31) 


The  acceleration  vector  d2a/dt2  has  components  ..id  t 
which  can  be  related  to  the  components  of  d  a/dt*  given  by 
Equation  (B-14)  as  follows 


+  *  ryJ  *  Ceos  (x,£)  +  ricos  (x,n)  +  £cos  (x,5) 


Vo  +  (rxo  “  Pzo) 


=  Coos  (y,£;)  +  ncos  (y,n)  +  coos  (y,C)  >  (B-32) 


+  (py  -  qx  )  =  £cos  (z,£)  +  ficos  (z,n)  +  ccos  (z,C) 


The  direction  cosines  in  the  above  equation  are  obtained  from 
the  transpose  of  [A],  Equation  (B-2) . 

With  the  following  notation 


r'  *  exx  +  eyy  4-  e2z 
and  Equation  (B-14)  we  have 

m?'  x  fit =  [y(zo +  pyQ  “  “  i(y0  +  rxo -  pzo}] 

dt 


(B-33) 


+  mey  [z(5t0  +  qiQ  -  ryQ)  -  x<»0  +  pyQ  -  qkQ)] 
+  mez  [x(yQ  +  rxQ  -  piQ)  -  y(xQ  +  q'zQ  -  ryQ)J 


(B-34) 


With  regard  to  the  remaining  term  in  Equation  (B-29) ,  we 
have  from  Equation  (B-22) 


d»o 

d 

r  *.  x  sii 

dt 

"  dt 

p  dt 

We  first  evaluate  and  then  differentiate 


=  J  p  *  x  dm  =  Jp'  x  (w  x  p  •  )dm  =  j  ju>(p'  *  p' )  -  p '  (p '  •  u)  j  dm 

=  |  u>(x2  +  y2  +  z2)dm  -  j  (exx  +  e^y  +  ezz)  (px  +  qy  +  rz)dm 
=  ex  j^p  (y2  +z2)dm  -  q  j  xy  dm  -  r  J  xz  dm 


-p  |  xy  dm  +  q  j  (x2  +  z2)  dm-  r  j  yz  dm 


-P 

. 


xz  dm  -  q  [  yz  dm  +  r  [  (x2  +  y2)dm 


]  (B-: 


The  moments  and  products  of  inertia  are  defined  as  follows: 
*xx  =  {  <y2  +  Ixy  =  {  xy  dm 

Iyy  =  j  (*2  +  z2)dra  Ixz  =  |  xz  dm  * 

Izz  =  {  (*2  +  y2)dm  Iyz  -  j  y2  dm 


» 

=  (X 

i 


2  2 

*  +  z  j  dm  I 


=  f  (x2  +  y2) dm  I 

* 


(B-36) 


so  that  Equation  (B-35)  becomes 

3°  -  K  Oxx  -  ** 


xy  "  rIxz. 


+  ey  l"  pIxy  +  qIyy  ~  rIyz] 
+  K  ['  -  <3lyz  +  rIzz] 


(B-37) 


In  body  coordinates  the  momenta  and  products  of  inertia  are  not 
functions  of  time  for  a  rigid  body.  With  the  help  of  the  rela¬ 
tionships 


we  find 


"dtT  =  ex  L^>Ixx  "  qIxy  "  ^Ixz  +  rq(Izz  ~  Iyy) 

-  (q2  -  r2)Iyz  -  p(qlxz  -  rlxy>] 

+  ®y  ['  PJxy  +  ^yy  -  ^yz  +  rP(Ixx  '  Izz> 

-  (r  -  p  ) Ixz  -  q  <rIxy  “  P1yz^] 

+  K  [-  Kz  -  Vyz  +  ixzz  +  M'V  -  Jxx> 

-  (P2  -  q2Uxy  -  r(plyz  -  qlxx>]  (B-39) 


Equating  the  components  of  Equation  (B— 31 )  to  the  sum  of 
Equations  (B-34)  and  (B-39)  yields  the  three  desired  equations 
of  rotational  motion 


M„  -  m  [y(Sc  +  pyD  -  qxQ)  -  z  <yQ  +  -  p !„>]  +  plxx  -  qlxy  -  il, 

+  V  '  (q2-r2)1yz'plqIxz-rIxy> 

M  -m  [7(x0+qi0-ri0)  -;(z0  +  py0-qi0)]  -  plx  +ql  -rl, 


+  ‘P*1**-1..’  -  (r2  *  P2)Ixz-'5(rIxy-PIyz) 


Mx  -  m  [  x  (y0  +  rxQ  -  piQ)  -  y  (*e  ♦  qi0  -  rx0)  ]  -  Pl„  -  qlyit  + 


r\  r\ 

+  pq  (I  —  I  )  —  (p^  —  q  '  I  -  r  (pi  -  ql  ) 
*4V  YY  xx  ^  4  xy  ^  yz  4  xz' 


(B-40) 


The  rotational  equations  of  motion  in  the  inertial  coordinate 
system,  using  Equation  (B-32) ,  are 

m  £  y  cos  (z,^)  -z  cos  (y,  £)]  i  +  m  [  y  cos  (z,n)  -z  cos  (y,n)]ri 

+  m[ycos  (2,0  -  z  cos  (y,C)]  H  +  !XXP  “  JXy4  “  rXBr 

-  Mx-rq(Izz-V  +  (q2-r2,Iy«  +  p(qIx.-rIxy> 

(B-41) 

m  [z  cos  (x,£)  -x  cos  (z,5)]|j  +  ra^z-  cos  (x,n)  -x  cos  (z,n)]fi 
+  m  [z  cos  (x,c)  -  x  cos  (z , ;)]  5  -  Ixyp  +  Iyyq  -  Iyzr 

-  \  -  lIxx  -Zzz>  +  fr2  -  P2)  *xz  +  q  (rIxy  *  V 

(B-42) 


>*»•  *  ***-.<i4  1 


.wmc ■w*Kcer;sr:-w^  '*Vr 


... 


COS 


„[*  cos  <y,U  -y  cos  U.5)]E  +  >"[X  cos  (y,n)  'V 

+  m  ^  x  cos  (y»  t)  -  Y  cos  &  “  V  *y2^  zz* 


(x,n)]n 


2  -2'i....  +  r(pi„!,-qixa> 

(B-43) 


=  Mz  -  Pq  dyy  -  Ixx)  +  <P  "  <*  >  ■ Axy  T  ‘  '*‘y*  ”'** 


% 
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B-6  OVERALL  SOLUTION  OF  EQUATIONS  OF  MOTION 


Certain  input  data  and  initial  conditions  are  to  be 
specified  before  any  integration  of  the  equations  of  motion 
can  be  performed.  First  the  input  data  must  include  the  body 
mass,  m,  and  the  moments  and  products  of  inertia.  Secondly, 
the  offset  of  the  origin  of  the  body  coordinate .system  from 
the  center  of  mass  must  be  specified.  We  also  need  the  initial 
body  orientation  as  specified  by  the  three  initial  angles. 

With  regard  to  initial  conditions,  the  initial  position 
(£,  C)  of  the  origin  of  the  body  coordinate  system  and  its 
velocity  (£,  n,  £)  must  be  specified  together  with  the  initial 
angular  velocity  components  p,  q,  and  r.  We  may  then  integrate 
Equations  (B-16)  through  (B-18)  and  Equations  (B-41)  through 
(B-43)  to  determine  £,  n,  S,  p,  q,  and  r  as  functions  of  time. 
Simultaneously  we  must  integrate  Equation  (B-l)  to  keep  track 
of  the  body  orientation.  The  aerodynamic  force  and  moment 
components  are  calculated  step-by-^tep  during  the  integration 
of  the  trajectory. 
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